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This  work  deals  with  the  analytical  and  numerical  modeling  of  two 
different  types  of  solar  cells:  solar  cells  made  of  crystalline 
silicon  and  thin-film  single- junction  amorphous  silicon  solar  cells. 
The  modeling  work  on  crystalline  silicon  solar  cells  is  done  by 
deriving  optimum  design  principles  for  high  conversion  efficiency. 
Optimization  of  the  performance  through  minimization  of  the 
recombination  losses  is  emphasized.  The  effect  of  certain  three- 
dimensional  geometries  on  the  operation  of  crystalline  silicon  solar 
cells,  and  other  bipolar  devices,  is  analyzed.  A reciprocity  theorem 
is  presented  that  relates  the  current  of  the  device  under  illumination 
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with  the  normalized  minority-carrier  density  under  bias  in  the  dark. 

In  the  same  context,  an  experimental  technique  is  described  and 
demonstrated  that  allows  measuring  the  minority-carrier  lifetime  and 
diffusivity  in  the  base  of  a solar  cell. 

The  modeling  of  amorphous  silicon  solar  cells  emphasizes  the 
interpretation  of  numerical  solutions  for  the  device  operating  in  the 
steady  state  under  optical  and  terminal  excitation.  These  solutions 
are  obtained  through  a computer  program  that  solves  the  Poisson  and 
the  current  continuity  equations.  Based  on  this  interpretation  we 
propose  an  equivalent-circuit  model  for  a-Si:H  solar  cells.  This  model 
clarifies  the  concept  of  the  limiting  carrier  and  enables  derivation 
of  analytical  formulas  for  the  cell  photocurrent  under  strongly  and 
weakly  absorbed  light. 

Finally,  we  compare  the  physics  and  the  modeling  of  crystalline 
and  amorphous  silicon  solar  cells  focusing  on  transient  experiments 
and  on  steady-state  recombination. 
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CHAPTER  ONE 
INTRODUCTION 

This  dissertation  mainly  deals  with  the  analytical  and  numerical 
modeling  of  crystalline  and  amorphous  silicon  solar  cells,  although 
some  results  are  general  enough  for  application  to  bipolar 
transistors,  optical  detectors  and  other  allied  devices.  It  consists 
of  two  parts.  The  first  part,  which  treats  crystalline  silicon, 
develops  optimized  modeling  and  consequent  design  principles  for  high 
efficiency  solar  cells  and  emphasizes  minimization  of  recombination 
losses.  The  first  part  also  provides  solutions  for  several  three- 
dimensional  transport  problems  in  crystalline  silicon  that  have 
implications  for  device  characterization  and  that  enable  extension  of 
the  optimum  design  principles  to  three  dimensional-solar  cells.  In 
addition,  an  experimental  technique  is  presented  and  demonstrated  that 
allows  characterization  of  silicon  solar  cells  by  measuring  the 
minority-carrier  lifetime  and  diffusion  coefficient  in  the  principal 
region  of  the  device. 

The  second  part  is  on  the  modeling  of  thin  film  amorphous  silicon 
(a-Si:H)  p/i/n  solar  cells.  To  enable  this  modeling,  a computer 
program  is  developed  that  solves  numerically  the  complete  system  of 
Shockley  equations.  These  solutions  are  treated  as  computer 
experiments  designed  to  reveal  the  physical  mechanisms  controlling 


1 


2 


charge  transport  and  quantum  efficiency.  Based  on  the  physical 
interpretation  of  the  numerical  solutions  we  propose  analytical  models 
for  the  terminal  current-voltage  and  photocurrent-voltage  relations. 
Uniformly  absorbed,  strongly  absorbed,  and  AMI . 5 light  are  studied. 

The  interdependencies  among  such  internal  variables  as  electric  field, 
recombination  rate,  and  carrier  densities  are  discussed.  A major  point 
in  the  analysis  of  a-Si:H  solar  cells  is  the  separation  of  currents 
into  photogenerated  and  back- inj ected  components.  This  separation 
clarifies  the  concept  of  the  limiting  carrier  and  helps  in 
understanding  the  physical  mechanisms  constraining  the  quantum 
efficiency.  It  also  enables  the  derivation  of  an  equivalent-circuit 
model  that  qualitatively  accounts  for  the  dependence  of  internal 
variables  on  optical  and  terminal  excitation.  Finally,  we  compare  the 
physics  and  modeling  of  single  crystal  and  amorphous  silicon  solar 
cells.  We  discuss  the  differences  in  the  physical  mechanisms  that 
control  charge  transport  in  the  two  devices  and  the  applicability  to 
amorphous  silicon  solar  cells  of  concepts,  like  minority-carrier  and 
low-level  injection,  borrowed  from  silicon  devices. 

The  work  that  follows  is  presented  in  five  chapters.  Chapter  two 
contains  a reciprocity  theorem  and  topics  in  three-dimensional 
transport  in  silicon  devices.  Chapter  three  is  on  the  optimized 
modeling  of  silicon  solar  cells.  In  Chapter  four  a new  experimental 
method  is  demonstrated  that  allows  measuring  minority-carrier  lifetime 
and  the  diffusion  coefficient  in  the  base  of  a silicon  solar  cell.  In 
Chapter  five  numerical  and  analytical  results  on  a-Si:H  solar  cells 
are  presented.  In  Chapter  six  we  compare  the  physics  and  the  modeling 
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of  crystalline  silicon  and  amorphous  silicon  solar  cells.  Chapter 
seven  summarizes  the  conclusions  reached  in  this  dissertation. 
Appendices  treat  particular  problems  arising  in  the  development  of  the 
various  chapters.  Each  chapter  provides  references  to  the  literature 
of  previous  related  work. 


CHAPTER  TWO 

RECIPROCITY  THEOREM  AND  SOLUTION  OF  THREE-DIMENSIONAL  TRANSPORT  PROBLEMS 

USING  ONE -DIMENSIONAL  MODELS 

2 . 1 Introduction 

This  chapter  shows  that  in  certain  three-dimensional  devices  a 
number  of  output  parameters  can  be  determined  or  related  to  the  input 
without  solving  the  actual  boundary  value  problem.  These  results  are 
derived  by  exploiting  the  symmetry  of  the  structures  involved  as  well 
as  existing  solutions  for  related  one -dimensional  problems.  Several 
conclusions  follow,  regarding  the  modeling  and  the  characterization  of 
devices,  which  are  three-dimensional  either  by  construction  (bipolar 
transistors,  diodes)  or  by  intention  to  improve  performance  (solar 
cells,  detectors). 

The  work  is  presented  in  two  sections.  Section  2.2  presents  a 
reciprocity  theorem  that  allows  calculation  of  the  internal  quantum 
efficiency  of  an  illuminated  p-n  junction  by  using  the  solution  of  the 
transport  equations  of  the  same  device  in  the  dark.  This  theorem  is 
used  in  Section  2.3  which  treats  transport  in  planar  devices  having 
mixed  boundary  conditions  at  one  surface.  The  transit  time  and 
transport  factor  of  shallow  planar  transistors  and  diodes  are  found  to 
be  the  same  as  in  their  one -dimensional  counterparts.  With  respect  to 
solar  cells  and  photodiodes  we  conclude  that  texturizing  the  cell 
surface  to  increase  the  effective  absorbance  leaves  the  base  bulk 
recombination  unchanged.  Finally,  partitioning  the  emitter  to  reduce 
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its  area  and  recombination  losses  results  in  the  same  open-circuit 
voltage  as  of  the  one -dimensional  base-limited  solar  cell  having  the 
same  base  and  a continuous  emitter. 

2 • 2 Generalized  Reciprocity  Theorem  For  Semiconductor  Devices 

Here  we  present  a reciprocity  theorem  that  relates  the  device 
short-circuit  current  induced  by  a carrier  generation  source  to  the 
minority- carrier  quasi -Fermi  level  in  the  dark.  This  theorem  is  an 
extension  of  Donolato's  [1]  reciprocity  relation  for  the 
photogenerated  charge  collection  by  a p-n  junction.  His  proof  was 
restricted  to  regions  where  only  the  lifetime  was  allowed  to  vary  with 
position.  We  extend  this  proof  to  p-n  junction  devices  having 
three-dimensional  space  dependencies  of  all  material  parameters.  In 
addition,  we  prove  that  reciprocity  holds  also  for  transients  and 
state  another  form  of  reciprocity  relating  quasi-Fermi  levels  and 
point  generation  in  a quasineutral  region  under  low  injection 
conditions . 

The  conventional  reciprocity  theorem  relating  input  voltages  and 
output  currents  in  the  base  of  a transistor  was  proven  by  W.  Shockley, 
M.  Sparks  and  G.  K.  Teal  [2]  in  one  of  the  early  papers  on  p-n 
junction  transistors.  Here  we  use  Shockley's  procedure  and  properly 
generalize  in  order  to  prove  reciprocity  for  charge  collection. 

The  principal  quasineutral  region  of  the  diode,  where  carrier 
generation  and  recombination  takes  place,  is  bounded  by  the  junction 
Sj  and  the  rest  of  the  surface  SL  (lateral  and  back  surface) 
characterized  by  a position  dependent  surface  recombination  velocity 
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s'.  The  region  is  assumed  to  be  p-type  in  low  injection.  Extension  to 
n-type  is  straightforward.  The  space  charge  region  edge  Sj  is  assumed 
bias  independent.  Later  in  this  section  we  will  elaborate  on  the 
validity  of  this  assumption.  Two  different  excitation  cases  are 
considered.  In  case  1 the  junction  is  forward  biased  by  a voltage  V 
in  the  dark.  In  case  2 the  junction  is  short  circuited  and  a point 
generation  source  of  strength  G is  applied  at  point  r' . If  J,  n and 
J',n'  are  the  minority  (electron)  current  density  and  excess  minority 
carrier  concentration  for  cases  1 and  2 respectively,  then  the 
boundary  conditions  for  each  case  are 

n<rJ> 

~"(rj)  " exp(eV/kT)  -1  , n'(rj)  = 0,  (2.1) 


where  rj  is  a point  on  the  junction  edge  Sj  and  n0  is  the  equilibrium 
minority-carrier  concentration.  The  continuity  equation  for  cases  1 
and  2 is 


V- J 


y-  + GS(r-r'),  (2.2) 


where  S( r - r')  represents  a unit  point  carrier  generation  source 
applied  at  r = r' , and  r is  the  lifetime.  Because  of  low  injection, 
the  electric  field  E,  the  mobility  /i  and  the  lifetime  r are 
independent  of  the  excitation.  The  system  is  linear. 

Let  vector  A be  defined  by  the  equation 


A 


n' 


n. 


J. 


(2.3) 
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The  surface  integral 


A - ds 


, nn'  , 
-e(— — s' 


n„ 


n’n 

ru 


s ' ) ds 


(2.4) 


From  (2.1)  this  integral  on  Sj  becomes 


[exp(eV/kT) -1] 


J'  - ds 
JSj 


[exp(eV/kT)-l]  I', 


(2.5) 


where  I'  is  the  current  due  to  the  point  source  ecxiting  the  p-type 
region  through  the  junction.  Therefore,  if  S is  the  entire  surface 
enclosing  the  quasineutral  region,  then  from  (2.4)  and  (2.5)  the 
total  surface  integral  of  A becomes 


A-ds  = A-ds  + A-ds  = [exp(eV/kT) -1]  I'.  (2.6) 

Js  JsL  JSj 

By  using  (2.2),  the  Einstein  relationship  for  non- degenerate 
minority  carriers,  D =*  (kT/e)/x,  and  the  expression  for  the  electric 
field  experienced  by  the  minority  carriers  E =-  -(kT/e)V  In  nQ , the 
divergence  of  A can  be  written  as 


V-A  =■  V-(—  J')  - V • (— — J)  - 
no  no 

~ [Vn  - nV  In  nQ]-J'+  -JJ-  V-J'-  — [Vn'  - n'V  In  nQ]-J  - — V-J  = 
o no  no  no 

( Vn  + -pr  nE)  • (e/in' E + eDVn')  - e-^-  (-  + GS(r-r')) 

llo  no  r 

- -J-  (Vn'  + -J-  n'  E)  • (e/mE  + eDVn)  + e—  (-  — ) = 
no  no  r 
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~jr  ^kT(  Vn  + -r™  nE)-(-J-  n'E  + Vn')  - e — - e—  GS(r-r')) 

no  kl  kT  nQr  nQ  v ' ' 


— MkT(Vn'  + n'E)-(-4?nE  + Vn)  + e-^- 


n 


kT 


- e — GS (r  - r') • 


(2.7) 


By  using  Gauss's  theorem  and  from  Eqs.  (2.6)  and  (2.7) 


J A ds  ‘ I 


A • ds  = J V • A dv , 

[exp(eV/kT)-l]I  '=  - e f — G5(r-r')dv  = -e  G.  (2  8) 

J no  n0 (r  ) 1 

In  (2.8)  the  ratio  n(r')/nQ(r')  is  the  normalized  excess  minority- 
carrier  density  at  r - r',  Because  of  the  linearity  of  the  system  this 
density  is  linearly  proportional  to  excess  carrier  density  at  the 
junction  edge  of  the  region:  n(r')/n0(r')  = ms(r')  [exp(eV/kT) -1] , 
where  ms(r')  is  the  saturation  value  of  the  normalized  excess 
minority-carrier  density  at  r=r',  in  the  dark.  According  to  (2.8), 
therefore,  the  collection  probability,  Cp,  of  a carrier  generated  at 
point  r'  is 


CD(r')  = ms(r') . 


(2.9) 


Equation  (2.9)  is  similar  to  Donolato's  final  expression.  It  is  more 
general  because  we  include  position  dependencies  of  n,  D and  n0 . 

In  a similar  manner,  one  could  prove  another  form  of  the 
reciprocity  theorem  applying  to  a quasineutral  region  in  low 
injection:  the  minority-carrier  quasi-Fermi  level  at  point  1 due  to  a 
point  source  applied  at  point  2 is  the  same  as  the  minority-carrier 
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quasi-Fermi  level  at  point  2 due  to  the  same  source  applied  at  point 
1.  These  theorems  also  hold  for  transients.  This  is  proved  by- 
replacing  n,I  and  n',I'  in  (2.2)  by  their  Laplace  transforms  and  1/r 
by  s + 1/r  for  zero  initial  conditions.  Here  s is  the  Laplace 
transform  variable.  The  result  is  that  (2.9)  for  transients  becomes 


Equation  (2.9)  could  be  useful  in  solar  cell  and  photodetector 
theory  since  it  allows  the  calculation  of  short-circuit  current  from 
the  analysis  of  the  device  under  dark  conditions: 


where  1,2  represent  the  quasineutral  emitter  and  base  of  the  cell. 
Equation  (2.11)  ignores  generation  in  the  space  charge  region  and 
holds  under  the  same  assumptions  underlying  the  proof  of  the 
reciprocity  theorem.  Equation  (2.11)  was  derived  on  the  assumption  of 
a bias  independent  boundary  Sj . In  a realistic  device  Sj  moves  with 
bias.  This  displacement,  however,  is  of  the  order  of  the  Debye  length 
and  in  this  transition  area  the  minority  carrier  Fermi  level  is  flat 
and  Cp(r)  - ms(r)  = 1.  Therefore,  the  assumption  of  a bias  independent 
Sj  does  not  restrict  the  applications  of  the  theorem  unless  the  quasi 
neutral  region  thickness  is  comparable  to  the  Debye  length. 

The  major  assumption  in  deriving  (2.8)  was  low  level  injection 
so  that  the  system  could  be  treated  as  linear.  Low  injection 
guarantees  linearity  of  the  system  even  in  the  presence  of  the  effects 
of  heavy  doping  or  graded  chemical  composition.  In  this  case,  in 


Cp(r' ,s)  - ms (r ' ,s) . 


(2.10) 


(2.11) 
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equation  -Vlnn0  — eE/kT,  E will  be  the  minority  carrier  field  which 
includes  the  quasi-fields  arising  from  the  position  dependent  energy 
gap,  electron  affinity,  and  density  of  states. 

Therefore,  the  original  Shockley's  proof  of  reciprocity  can  be 
extended  to  heavily  doped  three-dimensional  base  regions  to  obtain 
relations  between  input  voltages  and  output  currents.  In  this  sense, 
Shockley's  proof  is  more  general  than  the  recently  published  [3]  proof 
of  reciprocity  in  one -dimensional  transistor  base  regions  which  takes 
into  account  heavy  doping  effects. 

The  practical  significance  of  reciprocity  is  that  it  can  reduce 
the  number  of  unknowns  in  a problem  and  reduce  the  amount  of 
calculations  needed.  As  an  example,  in  the  recent  paper  by  del  Alamo 
and  Swanson  [4] , the  quantum  efficiency  of  an  emitter  could  have  been 
found  by  using  (2.11),  since  the  minority-carrier  density  as  a 
function  of  position  was  already  numerically  calculated  under  dark 
conditions.  In  addition,  their  forward  and  reverse  transport  factors 
and  saturation  currents  are  related  by  the  equation  afJof  = arJor. 
Thus,  the  present  reciprocity  theorem  can  find  use  in  the  measurement 
of  such  device  parameters  as  the  surface  recombination  velocity  at  a 
polycrystalline  silicon  contact. 

Although  not  discussed  here,  reciprocity  relations  enter  the 
physics  of  many  different  phenomena,  such  as  heat  conduction  or 
thermoelectric  effect.  These  relations  are  based  on  fundamental 
considerations  relating  to  detailed  balance  under  thermal  equilibrium 
and  Onsager's  principle  of  microscopic  reversibility. 
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2 • 3 Exact  Solution  of  Three-Dimensional  Transport  Problems  Using 

One-Dimensional  Models 

This  section  shows  that  in  certain  three-dimensional  devices  a 
number  of  output  parameters  can  be  determined  by  exploiting  existing 
solutions  for  related  one -dimensional  problems.  The  devices  considered 
are  mainly  shallow  planar  structures,  their  three-dimensionality 
arising  from  mixed  boundary  conditions  or  non-uniform  excitation.  The 
devices  either  extend  laterally  to  infinity  (Fig.  2.1(a))  or  they  end 
at  a passivated  lateral  surfaces  of  arbitrary  shape  forming  thereby 
the  general  orthogonal  cylindrical  structure  shown  in  Fig.  2.1(b).  By 
passivated,  we  mean  that  the  surface  recombination  velocity,  S,  is 
zero.  The  devices  have  one -dimensional  counterparts  having  the  same 
volume  properties  but  uniform  boundary  conditions  at  the  front  and 
back  planar  surfaces.  Each  device  and  its  one -dimensional  counterpart 
are  described  by  a common  transport  equation  and  by  the  same  values  of 
certain  parameters,  such  as,  the  decay  rate  time  constant,  the  transit 
time,  the  transport  factor  and  the  open- circuit  photovoltage.  Existing 
closed-form  expressions  for  these  parameters  for  the  one -dimensional 
case,  therefore,  immediately  apply  to  three-dimensional  structures. 

In  addition,  we  describe  analytically  certain  measures  of  current  and 
charge  for  devices  having  non-planar  front  and  back  surfaces.  All 
these  results  are  based  on  a second  order  linear  differential  equation 
which  is  derived  from  the  continuity  equation  and  is  presented  in  the 


next  subsection. 
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Fig.  2.1.  Three-dimensional  geometries  under  study. (a)  An  infinite 
slab  having  two  surfaces  at  z = 0 and  z = W and  an 
internal  plane  at  z=Z.  (b)  A cylindrical  structure  bounded 
by  a lateral  suface  orthogonal  to  front  and  back  surfaces. 
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2.3.1  Solution  of  the  Continuity  Equation  in  Planar  Symmetry  Cases 
Now  we  consider  that  the  sample  shown  in  Fig.  2.1(a)  is  excited 
by  a lazer  beam.  The  laser  beam  is  meant  as  an  example  of  excitation. 
The  following  analysis  is,  in  general,  valid  for  any  other  type  of 
excitation,  for  example,  thermal  generation,  and  minority-carrier 
injection  over  a potential  barrier.  The  sample  is  n-type  (extension 
to  p-type  is  straightforward). 

The  continuity  equation  for  the  hole  density  is 

3P 

(x.y.z, t)  = -V • Jp(x,y , z , t)+G(x,y , z , t)  - U(x,y,z,t),  (2.12) 

where  Jp  is  the  hole  flux  density  and  G and  U are  generation  and 
recombination  rates  in  excess  of  their  thermal  equilibrium  values. 
Integration  of  (2.12)  over  the  z - Z'  (0  < Z'  < W)  plane  yields 


3Q 

3t 


(Z\t)  = - 


dxdy  + 


f fcdxdy 


Udxdy 


rc° 

^px  3Jpy 
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dxdy  + 


f fo^y 

-ao  -oo 


f Udxdy. 


(2.13(a)) 


-oo  * -oo 


Here  Q(Z',t)  is  the  number  of  holes  per  unit  distance  on  the  z = Z' 
plane,  and  Jpx,  Jpy,  Jpz  are  the  components  of  Jp=Jpxx  + Jpyy  + Jpzz 
(x.y>z  are  the  unit  vectors  of  the  cartesian  coordinate  system).  The 
above  integrals  are  finite  for  realistic  generation  profiles. 
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By  Gauss's  theorem, 


r°° 


J - co 


dJ 


px 


3J 


py. 


3x 


dy 


dxdy  - J(Jpxx  + Jpyy ) ‘ ds , 

Poo 


(2.13(b)) 


where  T*,  is  the  infinite  circle  on  the  z=Z'  plane,  and  ds  is  the 
infinitesimal  vector  normal  to  I"*,  . The  integral  /roo  (Jpxx  + Jpyy)ds 
is  the  hole  flux  crossing  the  infinite  circle, per  unit  distance.  This 
integral  is  zero  for  any  realistic  generation  function,  since  the 
total  flux  crossing  the  entire  surface  at  the  infinite  is  bounded  by 
the  total  generation  rate.  Therefore,  (2.13(a))  becomes 


3Q 

dt 


3J 


_pz 

3z 


dxdy  + JJg  dxdy  - Jju  dxdy, 


(2.14) 


where  from  now  on  we  omit  the  (infinite)  limits  of  integration  for 
brevity.  Equation  (2.14)  is  a general  equation  and  relates  the  total 
charge  and  current  crossing  a plane  to  the  generation- recombination 
rate  on  the  same  plane.  We  focus  now  our  attention  to  planar  devices 
in  low  injection  where  internal  parameters  like  electric  field,  E, 
diffusivity,  D,  mobility,  n,  and  lifetime,  r,  are  functions  only  of  z 
and  are  independent  of  the  minority- carrier  density.  Expressing  the 
minority  flux  density  in  terms  of  drift  and  diffusion  of  the  excess 
holes,  P'  = P-P0,  and  the  recombination  rate  as  U = P'/r,  we  can 
rewrite  (2.14)  as 


3Q' 

3t 


„ *121  + . mEz)  30 : . (Siiiia  + 

3z  dz  3z  dz 


— ) Q'  + 


IJ 


G dxdy, 


(2.15) 


where  Q'  = //  P'dxdy  is  the  number,  per  unit  distance,  of  excess  holes 


15 


on  the  z=Z'  plane.  Equation  (2.15)  is  linear  provided  that  Ez , D,  /z 
and  r are  independent  of  the  Q'  because  of  low  injection.  The  same 
equation  also  holds  for  the  device  shown  in  Fig.  2.1(b)  since  the 
current  crossing  the  passivated  lateral  surface  is  zero. 

The  recombination  velocities  S]_,  S2  of  the  front  (z=0)  and  back 
(z=W)  surface  of  the  sample  enter  to  determine  the  boundary 
conditions  of  Eq.  (2.15): 

SX  Q'(t,0)  - D |^-  (t , 0)  , S2  Q'(t,  W)  =-  -D  . (2.16) 

Under  such  boundary  conditions  the  solution  of  Eq.  (2.15)  depends 
only  on  the  total  generation  rate  on  the  z = Z'  plane  and  not  on  G 
as  a function  of  x and  y.  A similar  conlusion  has  been  reached  in  [5], 
Knowledge  of  Q'(z,t)  does  not,  of  course,  provide  P'(x,y,z,t). 

However,  it  does  provide  the  total  minority  current  in  the  z 
direction,  Iz  =-eDdQ'/dz  + e/zEzQ'  , and  the  total  number  of  excess 
holes  on  the  z=Z'  plane. 

Therefore,  when  one  seeks  to  determine  either  the  total  minority 
current  or  the  number  of  excess  holes,  (2.15)  suffices,  provided  only 
that  the  transport  parameters  of  the  semiconductor  do  not  depend  on 
the  injection  level  and  are  functions  only  of  z.  Equation  (2.15)  is 
much  more  convenient  to  use  in  the  determination  of  (measurable) 
current  than  is  the  standard  continuity  equation.  In  the  following 
subsections  we  will  apply  (2.15)  to  derive  solution  for  three 
dimensional  transport  problems  in  planar  devices. 
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2.3.2  Photo luminescence  Decay  Response 

Photoluminescence  decay  is  one  of  the  methods  used  in  determining 
the  minority-carrier  lifetime  in  heavily  doped  silicon  [6].  The 
method  consists  of  exciting  a silicon  sample  by  a laser  pulse  and  thus 
generating  excess  electron-hole  pairs.  The  lifetime  is  derived  from 
the  decay  rate  of  the  excess  charge.  This  decay  as  a function  of  time 
is  obtained  by  monitoring  the  radiation  emitted  by  the  few 
electron-hole  pairs  that  recombine  through  radiative  recombination. 

In  principle,  this  method  is  more  reliable  than  any  other  methods 
because  it  involves  no  assumptions  regarding  the  band  structure  and 
the  transport  parameters  in  heavily  doped  silicon.  On  the  other  hand, 
the  method  has  the  drawback  that  surface  recombination  prohibits  the 
decay  from  being  purely  exponential.  Moreover,  the  effect  on  the 
decay  rate  of  the  size  and  shape  of  the  laser  beam  was  not  clarified 
previously.  Recent  papers  [7,8]  have  investigated  the  time  dependence 
of  the  luminescence  decay,  and  in  particular  the  influence  of  surface 
recombination  on  the  decay  rate  for  special  types  of  excitation. 

Ioannou  and  Gledhill  [7]  solved  the  continuity  equation  analytically 
for  the  total  excess  charge  in  the  case  of  a finely  focused  beam, 
while  Tyagi  et  al.  [8]  simplified  the  problem  to  one  dimension  by 
considering  uniform  excitation.  In  both  cases  it  is  found  that  the 
effect  of  the  surface  recombination  can  be  minimized  by  waiting 
several  lifetimes,  after  the  laser  is  turned  off,  before  measuring  the 
decay  rate.  Here  we  extend  the  analysis  for  the  general  case  of 
arkitrary  excitation,  and  prove  that  the  beam  size,  shape  and 
uniformity  does  not  affect  the  decay  rate.  An  upper  bound  for  the 
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error  due  to  surface  recombination  is  derived.  In  agreement  with 
references  6,7,8  the  effect  of  surface  recombination  is  found  to  be 
significant  only  in  the  initial  part  of  the  transient. 

We  now  apply  (2.15)  to  the  study  of  the  decay  of  the  excess 
charge  induced  by  a laser  pulse  in  a semi- inf inite  (W  = <»)  sample 
having  position  independent  doping  concentration.  For  low  injection 
minority  carriers  flow  by  diffusion  and  for  D independent  of  z (2.15) 
simplifies  to 


— “ D * — + Gt(t)(l-R)a  e‘az , (2.17(a)) 

3t  dz  r 

where  Gt(t)  is  the  total  photon  flux  of  the  laser  beam,  R is  the 
reflection  coefficient,  and  a is  the  absorption  coefficient. 

The  boundary  conditions  are 

Q'(0,z)  - 0,  Q'(t,«)  = 0,  Q'(t,0)S  = D 3Q^t,0)  , (2.17(b)) 

where  S is  the  surface  recombination  velocity. 

If  Gt(t)  is  the  unit  impulse  function,  Gt(t)  = <S(t),  then  Laplace 
transformation  of  (2.17(a))  with  time  yields 

D3-Q,^Z)  - (s  + — ) Q'  = - (l-R)a  e"QZ , (2.18) 

dz  T 

where  s is  the  Laplace  transform  operator.  The  boundary  conditions 
for  (2.18)  are  the  Laplace  transformation  of  (2.17(b)).  The  solution 


is 
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. aexp( -qz)  (1-R)  /T,  . rr 

Q (s  , z) ^ t—  - (r+a)ayD(l-R) 


-Da  +s+l/r 


x 


(7s+l/r  - yDDexp(-y(s+l/r)/D  z 
( -Da2+s+l/r ) ( -Dr2+s+l/T) 


(2.19) 


where  r — S/D. 

If  Q't(t)  is  the  total  number  of  excess  carriers  and  Q't(s)  its 
Laplace  transformation,  then 


Qt'(t) 


Q'(t,z) 


dz 


z)  dz 


_ 1-R (i-R)(r-ta)aD<yi+lA  - ypr)  (2  2Q) 

-Da2+s+l/r  /(s+l/r)  ( -Da2+s+l/r ) ( -DrVs+l/r ) 

Taking  the  inverse  Laplace  transform  of  (2.20),  we  obtain 

Qt'Ct)  ” Ta(l-R) (T-a) *lexp(-t/r) [a‘lexp(L2a2t/r)erfc(La/t7r) - 
r'lexp(L2r2t/r)erfc(LT/t/r) ] , (2.21) 

where  L — /Dr  is  the  hole  diffusion  length  and  erfc  is  the 
complementary  error  function.  This  result  holds  independently  of  the 
beam  shape  and  is  the  same  as  Equation  (8)  of  Ioannou  and  Gledhill  [7] 
derived  in  the  case  of  a finely  focused  beam. 

If  the  laser  pulse  has  a flux  N(t)  that  lasts  for  a finite  time 


T,  then  the  response  Q''t(t)  will  be  the  convolution  of  N(t)  with 

Q't(t): 
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Q"t(t)  - N(t')Q't(t-t')dt'  . (2.22) 

J0 

If  N(t)  is  a square  pulse  of  amplitude  N and  width  T,  then  (2.22) 
yields  for  t > T 

Q"t(t)  = ra(l-R)  (T-a)  ’1Nrexp(-t/r) 

x { - (L2a2  -1)  "1a'1  F(La,t)  + (rV-l)  'V1  F(rL,t) 

+ LfaV-l)'1  - (L2T2  -1) _1]  F(l,t)  },  (2.23(a)) 

where 

F(y,t)  - exp(T/r)exp(y2(t-T)/r)erfc(yV(t-T)/r)  - 

exp (y2  t/r ) erf c (yft/r) . 

If  S - T — 0,  (2.23(a))  reduces  to 

Q"t(t)  -Nr  (l-R)[exp(T/r)  -1]  exp(-t/r),  t > T (2.23(b)) 

which  is  a purely  exponential  function  with  time  constant  t. 

As  a more  general  statement,  if  S - 0 the  decay  of  the  excess 
charge  will  be  exponential  independently  of  the  shape  of  the  surface 
at  which  the  sample  terminates  and  independently  of  the  presence  of 
internal  fields.  This  can  be  proved  by  integrating  the  continuity 
equation 


3P' 

at 


-V-  Jr 


P' 


r 


(2.23(c)) 


over  the  volume  of  the  sample  and  applying  the  divergence  theorem: 
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HsT  - { Vds  - 7~  ■ (2.23(d)) 

J s 

Here  Q'  is  the  excess  minority-carrier  charge  and  the  integral  above 
is  the  hole  flux  at  the  surface.  This  vanishes  for  S = 0.  Therefore, 

dQ'  Q' 

d? r and  Q'(t)  - Q'(0)  exp(-t/r).  (2.23(e)) 


When  S > 0,  (2.23(a))  implies  that  Q"t(t)  is  not  exponential.  A 
measure  of  the  deviation  of  Q"t  from  the  exponential  function 
(Constant  exp(-t/r))  follows  from  considering  the  function: 

-d[  lnQ'  ' t(t)  ] dQ"t/dt  1 db(t)/dt 

dt  Q"t  ~ + b(t)  • (2-24) 

where  b(t)  is  the  negative  of  the  expression  in  the  curly  bracket  of 
(2.23(a)).  The  smaller  that  [db(t)/dt)/b(t) ] is  the  more  accurate  the 
estimation  of  t by  the  slope  of  lnQ"t  is. 

One  can  prove  that  [db(t)/dt]/b(t)  is  an  increasing  function  of 
S and  1/T.  As  a worst  case,  we  consider  infinite  recombination 
velocity  and  impulse  excitation.  Then  Q"t  becomes 

Q"t  =N'(1-R)  exp(L2a2t/r)  erfc(La/t/r)  exp(-t/r),  (2.25(a)) 

where  N'  is  the  strength  of  the  pulse.  In  such  a case 

db(t)/dt  _ 1/Jn  - La/t/r  exp  (L2q2  t/r ) erfc  (La/t/r ) L a/t/r  _ 

, , s n , I • (2.25(b)) 

b(t)  exp(L  a t/r)  erfc/(La(t/r ) ) t 


This  ratio  is  an  increasing  function  of  La  and  when  La  -*■  then 
[db(t)/dt]/b(t)  -*•  l/2t . Therefore,  in  the  worst  case,  by  measuring 
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the  slope  lnQ"^  at  time  t,  the  time  constant  derived  will  be 
r/[l+(r/2t)]  instead  of  r.  Thus  r/2t  is  an  upper  bound  for  the 
error  in  determining  r.  This  bound  is  much  lower  than  the  one, 
2.47r/t,  considered  by  Tyagi  et  al  [8].  Consequently,  the  maximum 
error  of  Dziewior  and  Schmid  [6]  due  to  surface  recombination  is 
approximately  r/[2(5r)]  = 10%,  since  from  their  data  it  appears  that 
they  waited  five  lifetimes  before  measuring  the  slope.  Thus  their 
results  are  reliable  in  so  far  as  surface  recombination  is  concerned. 

If  an  accuracy  of  20%  in  the  recombination  lifetime  is  required, 
then  0.8  - l/[l+(r/2t)]  which  implies  t=2t.  Now  the  question  arises, 
how  can  one  know  when  t is  equal  or  greater  than  2 r?  If  La(t/r)1/2>2 
(2.25(a))  can  be  simplified: 

Q"t(t)  * Q"t(0)  — gyP<-Vr)  . (2.26) 

(7rt/r)-*-/2La 

This  implies  that  Q"t(2r)  * 0.05  Q"t(0)  1/La.  If  L is  not  known  but 
one  can  estimate  its  range , then  a can  be  chosen  so  that  La  < 5 , 
implying  Q"t(2r)  >0.01  Q"t(0).  Therefore,  by  waiting  until  the 
response  drops  to  1%  of  its  initial  value,  t will  be  greater  than  2r 
and  the  error  less  than  20%. 

2.3.3  Transit  Time  and  Transport  Factor  in  Planar  Devices 

Here  we  will  use  (2.15)  to  show  that  shallow  planar  devices  have 
transit  times  and  transport  factors  that  are  equal  or  related  to  the 
transit  times  and  transport  factors  of  corresponding  one  dimensional 
devices.  The  devices  considered  are  made  on  substrates  which  are 
infinite  or  terminate  at  a passivated  (S“0)  cylindrical  surface 
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orthogonal  to  x-y  plane,  as  shown  in  Fig.  2.1.  Infinite  in  this 
context  means  that  the  lateral  surface  is  several  diffusion  lengths 
away  from  the  injecting  boundary.  The  device,  shown  in  Fig.  2.2,  is 
homogeneous  across  planes  parallel  to  x-y  planes . The  z component  of 
the  electric  field,  Ez,  minority-carrier  lifetime,  r,  and  mobility,  n, 
are  functions  only  of  z . The  back  surface  of  the  diode  is 
characterized  by  a uniform  recombination  velocity  S2  while  the  part  of 
the  emitter  not  occupied  by  the  junction  has  S-S^.  Initially  we  set 
equal  to  zero.  Later  we  discuss  the  effect  of  finite  on  the  transit 
time.  The  emitter  is  assumed  to  have  zero  thickness.  To  the  device 
shown  in  Fig.  2.2  corresponds  a one -dimensional  device  bounded  by  the 
same  planes  and  having  the  same  mobility,  lifetime,  and  field  strength 
as  functions  of  z,  as  the  device  in  Fig.  2.2.  The  two  diodes  differ 
only  in  the  boundary  conditions  at  the  top  surface:  the  diode  in 
Fig.  2.2  has  mixed  boundary  conditions  while  its  one -dimensional 
counterpart  has  uniform  boundary  conditions.  The  similarity  between 
the  two  problems  will  be  used  in  the  following  cases  to  show  that 
certain  three-dimensional  transport  problems  can  be  solved  in  terms  of 
one -dimensional  models  and  in  some  cases  can  have  closed  form 
solutions . 

The  base  transit  time  of  the  diode  of  Fig.  2.2.  is  given  by 


Q ( z ) dz 


-D(0)dQ(0)/dz  + /i(0)Ez(0)Q(0)  ’ (2,27) 


where  Q^.  is  the  total  number  of  excess  holes  (for  convenience  holes 
are  taken  to  be  the  minority  carriers)  stored  in  the  base  when  the 
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C MW  Ut*a 


'Back  Contact 


Fig.  2.2  A shallow  diode  formed  on  a semiconductor  slab.  The  dashed 
lines  connecting  the  emitter  with  the  back  contact  are  the 
minority-carrier  flow  lines,  while  the  dashed  line  that  is 
parallel  to  the  front  and  back  surface  defines  an  internal 
plane . 
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base  minority  current  is  ID.  In  (2.27),  ID  = I2(0)  = -D(0) dQ(0)/dz  + 
/i(0)Ez (0)Q(0)  because  there  is  no  surface  recombination  since  S]_  - 0. 
Here  Iz(z)  is  the  total  hole  flux  in  the  z direction. 

In  steady  state  and  under  no  volume  excitation  (2.15)  becomes 

2 

D(z)  + A(z)^2_  + b(z)Q  = 0,  (2.28) 

dz  dz 

where  A(z)  — dD/dz  - pEz  and  B(z)  = -d(/iEz)/dz  -1/r.  The  boundary 
condition  at  z - W is 

M( W)Ez(W)Q(W)  - D(W)dQ(W)/dz  - S2Q(W).  (2.29) 

The  other  boundary  condition  could  involve  the  total  charge  on 
the  front  surface,  Q(0) . Knowledge  of  Q(0)  would  require  a complete 
solution  of  the  three-dimensional  boundary  value  problem.  This  is  not 
necessary,  however,  if  only  rt  is  to  be  calculated.  Consider  the 
one -dimensional  diode  having  the  same  n(z) , Ez(z),  D(z) , t(z),  W and 
S2.  Here,  of  course,  the  junction  covers  the  entire  front  surface.  We 
will  prove  that  the  two  diodes  have  the  same  transit  time.  Let  Q' (z) 
be  the  excess  hole  density  and  J'(z)  the  hole  current  density  of  the 
one  dimensional  device.  The  quantity  Q' (z)  will  obey  the  same 
differential  equation  (2.28)  and  the  same  boundary  condition  (2.29), 
as  Q(z)  does. 

From  (2.29) 


dQ/dz  dQ'/dz  /i(W)Ez(W)-S2 

q “q7  m 


(2.30) 


Thus,  the  Wronskian  of  Q and  Q'  is  zero  at  z=W  and,  therefore,  is  zero 
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everywhere.  Consequently,  for  any  qiven  bias  V the  ratio  Q(z)/Q'(z)  is 
independent  of  z and  depends  only  on  the  geometry  of  the  diode.  Thus, 
for  any  value  of  z between  0 and  W 


Q(z)  dQ(z)/dz  -D(z)dQ(z)/dz  + /i(z)E(z)Q(z)  Iz(z) 

Q'(z)  “ dQ'(z)/dz  -D(z)dQ'  (z)/dz  + /x(z)E(z)Q'  (z)  = J'  (z) 


= Constant. 


(2.31) 


Therefore , 


Q(z) 

Q'(z) 


pW 

I Q(z) dz 



pW 

Q' (z)dz 

J0 


Qt  iz(°) 

Q't  J'(0)’ 


(2.32) 


where  and  Q^'  is  the  total  hole  number  and  total  hole  number  per 
unit  area  in  the  three-dimensional  and  one -dimensional  diode, 
respectively.  Equations  (2.31)  and  (2.32)  hold  for  any  values  of  S. 

S=0 , then  Iz(0)=Iq,  where  Iq  is  the  base  recombination  current  of 
the  three-dimensional  diode.  Let  J'  be  the  base  recombination  current 
density  for  the  corresponding  one -dimensional  diode.  Since  J'  = J(0), 
then  from  (2.32) 


Qt ' /J ' **  Tt'  ~ rt  = Qt^D-  (2.33) 

Therefore,  the  one -dimensional  transit  time,  rt' , is  the  same  as 
the  transit  time  of  the  shallow  junction  diode.  Since  no  assumptions 
were  made  regarding  the  shape  of  the  junction,  (2.33)  holds  for  any 
junction  morphology  if  the  thickness  is  very  small  compared  to  the 
diffusion  length  and  base  thickness  W.  Thus,  if  the  base  is  uniform, 
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closed  form  expressions  are  available  for  rt.  For  example,  if  S2  - «o 
and  W « (Dr)-*-/^,  then  rt  - W /2D.  In  the  general  case,  the  ordinary 
ferential  equation  (2.28)  would  have  to  be  solved  by  some  numerical 
methods.  This  is  much  simpler  than  solving  the  three-dimensional 
boundary  value  problem. 

The  derivation  of  (2.33)  was  based  on  the  assumption  of  zero 
recombination  current  at  the  front  surface.  Allowing  for  a finite  S]_, 
the  hole  recombination  current  becomes 

XD  = -D(°)  + M0)E  (0)Q(0)  + [Q(0)  - A(expg  - 1)  P0]S1,(2.34) 

where  A is  the  junction  area,  PQ  is  the  equilibrium  hole  density  at 
z=0  and  V the  applied  voltage.  The  term  in  the  square  bracket  in 

(2.34)  is  the  total  hole  number,  per  unit  distance,  on  the  uncovered 
front  surface.  Since  (2.31)  is  still  valid  if  the  charges  and  the 
currents  involved  are  replaced  by  their  saturation  values,  we  can 
write 

3Q(0)  Qq 

~D(0)  dz  +M(0)Ez(0)Q(0)  = — J0' , (2.35) 

ro 

where  JQ'  is  the  saturation  current  density  of  the  corresponding 
one -dimensional  diode,  and  Q0  is  the  saturation  value  (pre- exponential 
factor)  of  the  number  of  holes,  per  unit  length,  at  z=0.  Thus,  from 

(2.35)  and  (2.34)  we  obtain 

Qo  = do  + SiAP0)/(J'0/P0  + Sx).  (2.36) 


By  combining  (2.36)  with  (2.32)  and  (2.33)  we  obtain 
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rtpo  Qto  Qo  I0/po  + SqA 
VV  " Qto'  " “ Jo'/po  + Si* 

rt  1+SiAPo/I,, 

*V  1+Slpo/Jo'‘ 


Equation  (2.37)  provides  a way  of  measuring  the  product  S^Pq  if  the 
ratio  rt/rt'  is  known. 

Going  back  to  the  case  of  the  passivated  front  surface,  it  is 
possible  to  prove  that  (2.33)  holds  also  for  the  sinusoidal  steady 
state  excitation.  Here  rt'  - Qt'/J'  and  rt  - Qt/Id  where  Qt,  Jd,  Qt> 
Id  are  the  phasor  amplitudes  of  the  corresponding  quantities , of  the 
one  and  three-dimensional  devices,  oscillating  around  the  equilibrium 
point.  For  example,  ij)(t)  — Iq  + Id  exp(iwt).  It  is  assumed  that  the 
amplitude  of  the  a.c.  component  of  the  terminal  voltage  is  very  small 
compared  to  the  thermal  voltage.  To  prove  this,  it  is  enough  to 
rewrite  (2.38)  in  the  sinusoidal  steady  state 

2 

°(z)  + A(z)  — + [B(z)  - iw]  Q = 0 . (2.38) 

d z dz 


The  similarity  between  (2.28)  and  (2.38)  makes  evident  that  the  one 
and  three-dimensional  transit  times  are  the  same  in  the  sinusoidal 
steady  state,  also. 

We  may  think  of  the  device  in  Fig.  2.2  as  a transistor,  where 
emitter  and  base  retain  the  same  position  while  the  collector  extends 
under  the  base.  Now,  in  the  forward  active  mode  the  boundary  condition 
at  z = W becomes  S2  = 00  or  Q(W)  =0.  If  we  define  the  transit  time  in 
the  base  of  this  transistor  as 
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_ Q Q 

rtw  “ IZ(W)  ” -D(W)  dQ(W)/dz  ’ (2.39) 

then  based  on  the  discussion  preceding  (2.33),  one  can  prove  that 

r tw  ” r tw ' > (2.40) 

where  r^-w'is  the  transit  time  of  the  corresponding  one  dimensional- 
transistor.  Equation  (2.40)  is  valid  for  any  position  dependent 
recombination  velocity  in  the  front  surface  and  any  shape  of  the 
emitter. 

If  we  remove  the  assumption  of  the  infinitesimally  thin  emitter, 
we  can  no  longer  derive  closed  form  expressions  for  the  transit  time 
in  the  general  case.  It  is  possible,  however,  to  place  the  transit 
time  of  the  three-dimensional  base  between  upper  and  lower  limits.  In 
Fig.  2.3  a transistor  is  shown  having  an  emitter  of  finite  thickness. 
The  transit  time  can  be  written  as 

Qt  Qt  Qt(w_d)  Qt(W-d) 
rtw  “ IZ(W)  = Qt(W-d)  IZ(W)  > IZ(W)  rtw' (W‘d) > (2-41) 

where  Q^(W-d)  is  the  number  of  excess  holes  in  the  region  d < z < W 
and  (W-d)  is  the  transit  time  of  the  one -dimensional  transistor 
its  base  extending  from  z=d  to  z=W.  We  can  now  imagine  the 
junction  to  shrink  towards  plane  z=0  in  such  a way  that  its  original 
position  remains  an  equipotential  surface  (potential  in  the  sense  of 
the  minority-carrier  quasi-Fermi  potential) . This  can  be  practically 
done  by  extrapolating  the  minority-carrier  flow  lines  in  the  emitter 
and  drawing  the  equipotential  surfaces  that  are  orthogonal  to  them. 
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Fig.  2.3 


A planar  transistor  with  an  emitter  of  finite  thickness. 
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In  the  limit,  when  the  emitter  will  have  collapsed  on  the  front 
surface , the  transit  time  of  the  new  base  will  be  , the  same, 

according  to  (2.40),  as  in  the  one -dimensional  transistor  with  its 
base  extending  from  z=0  to  z=W.  For  a given  collector  current,  the 
charge  stored  in  the  base  of  this  hypothetical  transistor  will  be  more 
than  that  stored  in  the  actual  base  shown  in  Fig.  2.3  by  the  amount  of 
the  charge  stored  in  the  removed  emitter.  Thus  rtw  < rtw' (W) . From 
this  inequality  and  (2.41)  the  upper  and  lower  bounds  for  rtw  become 

'tw'(W-d)  < rt w <rtw'(W).  (2.42) 

Going  back  to  the  transistor  with  the  shallow  junction,  it  can  be 
seen  from  (2.31)  that  the  base  transport  factor,  BT,  of  this 
transistor  is  equal  to  the  transport  factor,  B’p'  , of  the  corresponding 
one -dimensional  transistor,  if  Sj-0.  If  Sx  is  finite,  then,  by  the 
same  procedure  that  led  to  (2.37),  we  obtain 


B'p 


l+AP0S1/I0 

l+P0Si/J0  ’ 


(2.43) 


where  IQ  and  JQ  are  the  collector  saturation  current  and  current 
density  of  the  three  and  one -dimensional  transistors  respectively, 
Iz(w)“Io[exP(eVAT)-l]  . The  assumption  of  a uniform  recombination 
velocity  on  the  base  surface  is  not  unrealistic  if  this  surface  is 
covered  by  a thin  polysilicon  layer  that  provides  isolation  of  the 
active  base  from  the  base  ohmic  contact.  Equation  (2.43)  can  in  fact 
be  used  to  determine  the  product  PQS^  characterizing  this  polysilicon 
layer  contact. 
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2-3.4  Base  Recombination  Currents  in  Three-Dimensional  Photovoltaic 
Devices 

In  Fig.  2.4  is  shown  a solar  cell  having  a textured  emitter  to 
improve  its  collection  efficiency.  The  saturation  value  of  the  hole 
particle  current  in  the  base  can  be  determined  without  solving  the 
three-dimensional  problem,  if  this  assumption  can  be  made:  a 

quasi -equilibrium  condition  exists  in  the  sense  that  the  minority  hole 
quasi-Fermi  level  is  position  independent  in  the  emitter  pyramids 
bounded  by  the  dashed  lines  from  above  and  by  the  dashed- dotted  line 
from  below.  This  assumption  is  good  when  (a)  the  dimensions  of  the 
pyramids  are  small  compared  with  an  average  base  diffusion  length  and 
(b)  either  the  base  thickness  is  much  greater  than  the  diffusion 
length,  or  the  back  surface  of  the  base  is  passivated  for 
minority-carrier  recombination.  In  many  silicon  solar  cells  having 
texturized  front  surfaces , these  conditions  will  hold  because  the 
dimensions  of  the  pyramids  are  a few  microns  whereas  the  diffusion 
length  and  cell  thickness  are  many  tens  of  microns . 

This  assumption  and  the  same  arguments  as  in  the  previous  section 
imply  a simple  determination  of  the  saturation  value  of  the  hole 
recombination  rate,  I0,  in  the  region  0 < z < W.  This  saturation 
current  of  the  structure  of  Fig.  2.4  equals  the  saturation  hole 
current,  Iq  , of  a related  one -dimensional  cell  having  a planar 
emitter,  the  same  area  and  a quasi-neutral  base  extending  from  z=0  to 
z=W  and  having  the  same  volume  parameters  and  back  surface 
recombination  velocity  S2.  For  completeness,  we  can  add  to  I0  the 
saturation  value  of  the  recombination  rate  in  the  quasi 
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Fig.  2.4 


A photovoltaic  device  with  a texturized  emitter. 
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equilibrium  pyramids:  /Pyramids  = [P0(r)/r (r) ] d3r  , where  PQ(r)  and 
r(r)  are  the  equilibrium  hole  density  and  lifetime  at  point  r.  Thus 
the  total  hole  saturation  particle  current  in  the  entire  base  can  be 
expressed  as 


Io  = V + /pyramids  (po/r)  d3r-  (2.44) 

2-3.5  Open-Circuit  Voltaee  and  Reciprocity  Relations  in  Planar 
Photovoltaic  Devices 

In  Fig.  2.5  is  shown  a solar  cell,  or  a photodetector,  having  a 
very  shallow  emitter  which  covers  only  part  of  the  front  surface.  The 
uncovered  front  surface  of  the  base  is  well  passivated  (S^O)  . The 
emitter  may  have  any  shape  and  area.  It  can  also  be  a strip  or  a 
number  of  strips.  The  illumination  is  assumed  independent  of  x and  y 
and  the  front-grid  coverage  negligible,  so  that  the  generation  rate  G 
is  a function  only  of  z:  G=G(z).  In  the  absence  of  emitter 

recombination,  the  open-circuit  voltage  of  this  cell  in  low  injection 
will  be  the  same,  under  the  same  illumination,  as  that  of  the 
one  - dimens ional  cell  having  the  same  base  and  a uniform  emitter 
covering  the  entire  front  surface.  To  prove  this  we  will  use  the 
reciprocity  theorem  presented  in  section  2.2  to  express  the 
collection  probability  Pc(r)  of  a minority  carrier  generated  at  point 
r as  Pc(r)=  Ps(r)/PG(r),  where  Ps  is  the  saturation  value  of  the 
excess  minority-carrier  density  at  r in  the  dark  and  P0  is  the 
equilibrium  minority-carrier  density.  Applying  this  relation  to  the 
cell  shown  in  Fig.  2 . 5 we  can  write  the  short-circuit  current  as 
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A planar  solar  cell  with  a shallow  emitter.  The 
illumination  is  uniform  in  the  x-y  plane. 
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base  J0 


Ps(x,y,z' ) 


W G(z')Qs(z') 


0 PQ(z') 


(2.45) 


where  Qs(z')  is  the  saturation  value  of  total  number,  per  unit 
distance,  of  excess  holes  in  the  z-z'  plane  in  the  dark.  If  Q's  is 
equivalent  quantity  for  the  corresponding  one -dimensional  cell,  then 
using  (2.32)  we  can  write 


where  IQ  and  JQ  are  the  saturation  current  and  saturation  current 
density  of  the  cell  shown  in  Fig.  2.5  and  the  one -dimensional  cell, 
respectively.  From  (2.45)  and  (2.46), 


Since  the  open-circuit  voltage  is  given  by  VQC  = (kT/e)Ln(Isc/IQ)  from 
(2.47)  it  follows  that 

Voc  - (kT/q)Ln(Isc/I0)  = (kT/e)Ln(Jsc'/J0' ) = VQC'  (2.48) 

If  the  front  surface  recombination  velocity  has  a finite  value  Sp, 
then  by  analogy  to  (2.37),  we  can  prove  that 


Qs(z')  - Qs'(z')I0/J0\ 


(2.46) 


where  Jsc'  is  the  short-circuit  current  density  one -dimensional  cell. 


ISC  JSC'  1 + APosl/Io 


1 + P0S1/J0' 


(2.49) 
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so  that  the  voltage  becomes 

V0C  = V'0C  ’ (kT/q)Ln[(l  + AP0S1/I0)/(1  + P0Si/J0')],  (2.50) 

where  PQ  is  the  equilibrium  density  of  holes  at  the  front  base 
surface.  Equations  (2.49)  and  (2.50)  provide  an  alternative  way  of 
measuring  P^.  Equations  (2.37),  (2.43)  and  (2.49)  show  that 
normalized  quantities  like  the  transit  time,  the  transport  factor  or 
the  ratio  Isc/^o  differ  from  the  analogous  values  of  the  corresponding 
one  - dimens ional  diode  or  transistor  by  a common  factor: 
(1+AP0S^/I0)/(1+S^P0/J0) . Of  course,  if  S^-0  these  values  become  the 
same.  In  the  particular  case  where  light  generates  carriers  only  at 
the  surface,  (2.47)  implies  a result  already  reached  by  Holloway, 
Brailsford  and  Sander  [9,10]  through  another  method.  This  result 
states  that  the  fractional  increase  in  the  collecting  area  of  a 
photodiode  due  to  lateral  diffusion  of  photogenerated  carriers  at  the 
surface  is  equal  to  the  fractional  increase  of  the  diode  saturation 
current  compared  with  that  of  the  one -dimens ional  diode.  Equations 
(2.48)  and  (2.49)  hold  also  for  the  interdigitated  [11,12]  or  point 
contact  [13]  cell  provided  the  back  diffused  contact  regions  are 
shallow.  In  such  cells,  however,  the  assumption  of  low  injection  could 
become  unrealistic  since  these  cells  usually  have  lightly  doped  bases 
and  are  subjected  to  high  concentration  ratios. 

The  fact  that  for  passivated  front  surfaces  the  ratio  Isc/I0  is 
emitter  geometry  independent  has  design  implications:  in  a solar  cell 
with  partitioned  emitter  for  a given  front  surface  coverage,  or  back 
surface  in  an  interdigitated  cell,  the  short-circuit  current  increases 
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when  finer  geometries  for  the  emitter  lines  or  point  contacts  are 
used.  This  results  because  IQ  increases  due  to  the  more  pronounced 
three-dimensional  effects. 

Reciprocity  can  also  be  used  to  relate  the  recombination  current 
of  the  two  different  devices  shown  in  Fig.  2.6.  Figure  2.6(a) 
illustrates  part  of  the  emitter  of  a planar  solar  cell.  This  emitter 
is  bounded  by  the  front  surface  which  is  assumed  passivated  with  S^=0, 
by  the  strip  ohmic  contact  and  by  the  injection  boundary.  The 
reciprocal  diode  is  shown  in  Fig.  2.6(b).  Here  a shallow  strip 
junction  replaces  the  ohmic  contact  and  the  boundary  at  z=-W  becomes  an 
infinite  ohmic  contact.  In  other  words,  the  injecting  boundary  and  the 
ohmic  contact  interchange  positions  from  Fig.  2.6(a)  to  Fig. 2. 6(b). 

The  bulk  properties  of  the  principal  regions  are  the  same  as  well  as 
sl-  Since  the  diode  shown  in  Fig.  2.6(b)  has  been,  to  some  extent, 
analyzed  [14]  it  is  of  interest  to  derive  reciprocity  relations  that 
would  express  the  recombination  current  of  the  device  in  Fig. 2. 6 (a)  in 
terms  of  that  of  the  diode  in  Fig.  2.6(b),  thus  eliminating  the  need 
to  solve  the  new  boundary  value  problem. 

We  introduce  again  the  one -dimensional  diode  that  corresponds  to 
the  device  shown  in  Fig.  2.6(b)  and  has  the  entire  z = 0 plane 
occupied  by  the  shallow  junction.  Let  Iao  and  Ibo  be  the  saturation 
value  of  the  total  particle  current  in  the  negative  and  positive 
direction  of  the  devices  in  Fig. 2. 6(a)  and  2.6(b),  respectively.  By 
J'bo  we  denote  the  current  density  of  the  corresponding  one- 
dimensional  diode.  Using  the  conventional  reciprocity  theorem 
relating  voltages  and  currents  in  the  base  of  a transistor,  we  derive 


38 


2-0 


/ I t«  ioN 

/ ♦ \ 

/ I V 


^Ohmlc  Contact 

jr  Si  -0 


2-W 


♦ 


Emitter 


fl  (W)  1 


L 


^Emitter-Base  Junction 

(a) 


z-o 


■7- 


Emitter 

S" 


S,  “0 


2-W 


il'.'o'N 

! \ 

I l 

f |»btW) 


Base 


(b) 


Ohmic  Contact 


Fig.  2.6  Reciprocal  bipolar  devices,  (a)  Part  of  the  emitter  of  a 
planar  solar  cell  showing  the  minority-carrier  flow  lines 
and  the  strip  contact,  (b)  A shallow  strip  p/n  junction  on 
a planar  base  having  an  ohmic  contact. 
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the  relation 


w°)  “ iboOO- 


(2.51) 


From  (2.31)  follows 


IboW/W0)  - Jbo'(W)/Jbo'(°)- 


(2.52) 


From  (2.51),  (2.52)  we  have 


Iao<°)  “ W°)  [Jbo'  (W)/Jbo'(0)]. 


(2.53) 


The  last  equation  expresses  the  recombination  at  the  ohmic 
contact  of  the  emitter  of  Fig.  2.6(a)  in  terms  of  the  recombination 
current  of  the  diode  in  Fig.  2.6(b)  and  the  recombination  currents  of 
the  one  dimensional  diode.  Analytical  expressions  similar  to  (2.53), 
though  more  complicated,  can  be  also  obtained  for  finite  Sb.  From 
(2.53),  from  relations  derived  in  reference  [14],  and  from  the 
reciprocity  theorem  we  can  write  an  approximate  formula  for  the 
recombination  rate  under  the  emitter  strip.  The  relative  increase,  R, 
of  this  rate  due  to  two-dimensional  minority-carrier  diffussion 
spreading  under  the  strip,  is  given  by 


where  L'  is  the  width  of  the  strip,  Lc  is  a characteristic  length, 
approximately  equal  to  the  minimum  between  an  average  minority- carrier 
diffusion  length  in  the  emitter  and  W.  In  conventional  solar  cells 
this  ratio  is  usually  negligibly  small.  In  some  of  today's  high 
efficiency  cells,  however,  especially  those  made  for  concentrator 
modules,  R could  have  non-negligible  values. 


R = (Lc/L')(Jbo'(W)/Jbo'(0)), 


(2.53) 
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2.3.6  Arbitrary  Boundaries 

So  far  we  have  discussed  planar  structures  where  the  only 
deviation  from  one-dimensionality  was  the  boundary  condition  at  the 
planar  surface.  Equation  (2.15),  however,  may  have  more  general 
applications.  Figure  2.7  shows  a quasi -neutral  region  limited  by 
arbitrary  surfaces  which  can  be  its  external  surfaces  or  the  injecting 
boundaries . The  device  is  either  infinite  in  the  x-y  direction  or 
terminates  at  a lateral  surface  orthogonal  to  the  x-y  plane.  If  the 
quasineutral  region  has  space -independent  transport  parameters,  then, 
in  the  steady  state  and  for  a generation  profile  of  the  form  G(z)=g 
exp(-az),  (2.15)  becomes 

2 

D - -2-  _ g exp(-az).  (2.53) 

dz  r 

Equation  (2.53)  has  a general  solution 

Q(z)  - Aexp( -x/L)  + Bexp(x/L)  + [g/(Da2 -1/r) ]exp(-az) , (2.54) 

where  A,  B are  constants  and  L=(Dr)-*-/^.  In  addition 
Iz(z)  - DAexp ( - x/L ) /L  - DBexp(x/L)/L  + Da [ g/(Da2 - 1/r ) ] exp( -az) . (2.55) 

To  determine  A and  B the  boundary  value  problem  would  have  to  be 
solved.  However,  even  with  A and  B undetermined,  (2.55)  provides  an 
analytical  expression  of  how  the  charge  and  the  current  in  the  z 
direction  changes  with  position. 
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A quasi-neutral  region  bounded  by  arbitrary  surfaces. 
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2 . 4 Conclusions 

In  the  first  part  of  this  chapter  we  presented  a reciprocity 
theorem  that  equates  the  collection  probability  in  the  base  of  a 
photovoltaic  device  to  the  normalized  excess  carrier  density  under 
electrical  bias  in  the  dark.  In  the  second  part  we  showed  that  planar 
devices  have  the  same  excess  carrier  decay  rate,  transit  time, 
transport  factor  and  open-circuit  voltage  as  the  corresponding  one- 
dimensional counterparts,  when  the  front  surface  recombination 
velocity  is  zero.  For  finite  values  of  the  recombination  velocity, 
analytical  formulas  were  derived  that  can  be  used  when  characterizing 
the  front  surface.  The  base  recombination  rate  per  unit  area  of  a 
texturized  emitter  solar  cell  is  the  same  as  the  recombination  rate  of 
the  same  cell  with  planar  emitter.  Partitioning  the  emitter  in  a base 
limited  cell,  to  reduce  the  area  and  the  associated  recombination, 
leaves  the  open-circuit  voltage  unchanged.  Finer  emitter  geometries 
result  in  higher  short-circuit  currents.  General  analytical  formulas 
were  also  derived  for  the  charge  and  current  of  non-planar  devices. 


CHAPTER  THREE 

SYSTEMATIC  DESIGN  THEORY  FOR  SILICON  SOLAR  CELLS  USING 
OPTIMIZATION  TECHNIQUES 

3 . 1 Introduction 

Upper  bounds  on  the  open-circuit  voltage  and  efficiency  of 
silicon  solar  cells  can  be  derived  by  considering  the  cell  performance 
as  constrained  only  by  detailed  balance  limitations  [15-18] . These 
bounds  are  general , in  the  sense  that  they  are  independent  of  the 
particular  device  design  and  geometry.  Such  derivations,  however,  have 
limited  design  implications.  This  work  will  use  a different  starting 
point  relating  to  bounds  on  solar-cell  performance  to  set  up  several 
practical  design  principles.  The  analysis  to  be  presented  is  not 
based  on  the  thermodynamics  of  interaction  between  radiation  and 
matter,  and  on  detailed  balance,  but  rather  is  focused  on  the 
transport  of  minority  carriers  in  low  injection.  It  is,  therefore, 
less  general  than  treatments  based  on  detailed  balance,  and  it  is 
restricted  to  one  spatial  dimension.  Our  analysis  applies,  however, 
for  arbitrary  doping  profiles,  drift  fields,  high- low  junctions  and 
includes  heavy  doping  effects;  moreover,  it  provides  explicit  rules 
for  design. 

Design  principles  follow  from  three  theorems.  Two  of  these, 
derived  here,  emphasize  recombination  in  the  non- illuminated  solar 
cell.  The  third  is  the  reciprocity  theorem  presented  in  the  previous 
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chapter.  The  theorems  exploit  linearity.  Thus  the  results  hold  only 
under  low  injection  conditions  in  base  and  emitter  regions.  For  the 
most  part,  the  work  concentrates  on  one -dimensional  models,  but 
extensions  are  made  to  accommodate  gridded  structures  on  the  back  and 
front  surface.  In  section  3.2,  lower  bounds  on  the  recombination 
current  are  derived  as  functions  of  surface  and  bulk  properties.  In 
section  3.3,  a theory  is  presented  that  allows  minimization  of  the 
recombination  current  for  any  given  thickness  and  (constant)  surface 
recombination  velocity.  In  section  3.4  the  model  developed  is 
compared  with  the  transparent  and  quasi -equilibrium  models,  and  bounds 
on  the  short-circuit  current  are  derived.  The  results  of  the  first 
three  sections  lead  to  optimum  design  principles,  presented  in  section 
3.5. 

3 . 2 Lower  Bounds  of  Dark  Recombination  Current  in  Quasi -Neutral 

Regions 

In  this  section  we  derive  upper  bounds  for  the  minority- carrier 
dark  recombination  current  of  regions  in  low  injection.  These  bounds 
affect  the  upper  limits  of  the  open- circuit  voltage  and,  consequently, 
of  the  conversion  efficiency.  The  following  analysis  focuses  on  the 
transport  parameters  of  minority  carriers  in  quasi -neutral  regions  and 
derives  the  absolute  bounds  on  the  recombination  current  by  letting 
the  lifetime  be  limited  only  by  the  fundamental  recombination 
mechanisms:  radiative  and  band-band  Auger  processes.  Low  injection  is 
assumed  throughout  this  work.  This  assumption  greatly  simplifies  the 
problem  since  the  differential  equations  describing  the  system  become 


linear. 
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The  analysis  presented  is  valid  for  any  doping  profile  in  one 
dimension.  It  accounts  for  the  presence  of  high-low  junctions  and  of 
heavy  doping  effects.  Extensions  to  three  dimensions  are  made  when 
possible.  Except  where  noted,  for  concreteness  we  treat  n-type 
silicon;  results  analogous  to  those  developed  below  hold  for  p-type 
silicon. 

The  two  basic  equations,  regarding  the  steady-state 
minority-carrier  (hole)  transport  and  recombination,  are 

Jp  =-eMP  dVp/dx  (3.1) 

dJp/dx  - -e(P  - PQ)/r.  (3.2) 

In  (3.1)  Vp  is  the  minority-carrier  quasi-Fermi  potential  referenced 
to  the  majority  carrier  (electron)  quasi-Fermi  potential  (which  is 
assumed  independent  of  position  [4]).  Also,  Jp , P,  PQ  and  r are  the 
minority  current  density,  carrier  density,  equilibrium  carrier  density 
and  lifetime,  respectively. 

The  densities  P and  PQ  are  related  to  Vp  by  P = P0exp(eVp/kT) . 
Using  this  equation  and  setting  u =>  exp(eVp/kT)  - 1,  then  from  (3.1) 
and  (3.2)  we  can  write 

d2u  d(DPQ)/dx  du  u 

+ " - 0,  (3.3) 

dx  DPq  dx  L 

where  L - JVr  is  the  minority  carrier  diffusion  length.  Equation 
(3.3)  holds  also  in  the  space-charge  region  (SCR)  of  a high-low  (H/L) 
junction  since  the  only  assumption  used  in  deriving  (3.3)  was  low 
injection.  This  is  expected,  because,  under  low  injection  the 
electric  field  there  is  independent  of  the  voltage  across  the 
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terminals.  Note  that  u is  the  fractional  excess  minority  carrier 
density:  (P-P0)/P0. 

In  a quasi -neutral  region,  where  the  majority  electron 
concentration,  N,  is  equal  to  the  ionized  donor  concentration,  we 

characterize  PQ  in  the  usual  manner: 

PG  - n^/Neff  (3.4) 

in  which 

Neff  “ nDD  exp(-6EG/kT) . (3.5) 

In  (3.4)  n^  is  the  intrinsic  carrier  density,  while  in  (3.5)  SEG  is 
the  effective  band- gap  shrinkage  (including  the  degeneracy  of  majority 
carriers) . 

In  the  SCR  of  a H/L  junction,  N(x)  (electron  density)  and  NDq(x) 
are  not  equal  and  6Eq  is  a function  of  both  N(x)  and  Nqq(x) , which 
makes  PQ(x)  hard  to  characterize.  Additional  difficulties  enter  for 
compensated  material.  To  overcome  this  problem,  we  assume  that  the  SCR 
is  so  thin  that  Vp  is  nearly  flat  through  it  and  that  the  integrated 
recombination  rate  through  the  volume  is  negligible.  This  is  a 
realistic  assumption  since  the  thickness  of  the  SCR  is  of  the  order  of 
the  Debye  length  in  the  low  doped  region  and  is  usually  a fraction  of 
a micron.  Then,  we  can  compress  the  H/L  junction  into  a point  where 
the  coefficients  of  (3.3)  may  have  step  discontinuities.  The 
continuity  of  u and  DPQdu/dx  — -Jp/e  can  serve  as  boundary  conditions 
connecting  the  solutions  of  (3.3)  in  the  quasi-neutral  regions  on 
either  side  of  the  H/L  junction  (point).  Thus,  Po  is  given  by  (3.4) 
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except  for  a number  of  points  equal  to  the  number  of  H/L  junctions. 

For  mathematical  simplicity  and  for  other  reasons  to  be  discussed,  D 
and  r are  also  assumed  as  functions  of  NDD  only.  From  (3.3)  and  (3.4) 
we  obtain 


d u d(D/Ne££)/dx  du  u 

7 + 7 “ 0.  (3.6) 

dx  D/^eff  dx  L 

Equation  (3.6)  is  the  same  as  the  homogeneous  part  of  (13)  of  del 
Alamo  and  Swanson  [4].  Note  that  this  equation  holds  between  two 
successive  H/L  junctions. 

The  assumption  that,  apart  from  temperature  T,  D is  a function 
only  of  Ndd  is  realistic  because  the  ionized  impurity  scattering  is 
the  principal  mechanism  that  determines  D;  D is  a weak  function  of  the 
fabrication  processing  at  least  for  dilute  doping  concentrations.  On 
the  other  hand,  the  lifetime  depends  strongly  on  fabrication.  Since 
we  are  mainly  concerned  here  with  lower  bounds  on  the  recombination 
current,  we  focus  mainly  on  fabrication  processing  that  gives  the 
longest  lifetimes.  Such  a lifetime  will  be  a function  only  of  NDD. 

If  d0(x)  is  the  saturation  value  of  the  recombination  current 
per  unit  area  in  the  region  extending  from  the  surface,  located  at 
x=0,  to  a distance  x,  as  shown  in  Fig.  3.1,  then  JQ(x)  = enj:2R(x)  , 
where 


J0(x)  1 JD  D du/dx 

R(x) — - . (3.7) 

e^i  enj.  u Neff  u 

Here  Jp  is  positive  when  it  flows  in  the  positive  x direction.  In 
words,  R(x)  is  the  pre-exponential  factor  (saturation  value)  for  the 
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Surface 


JQ(x)-en*  R(x,Rs) 


I 

I Junction 


r 


R(6,  Rs) 


Fig.  3.1  Illustration  of  the  region  under  study  and  definition  of 
the  function  R(x,Rs). 
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recombination  particle  current  in  the  interval  from  0 to  x, 
normalized  with  respect  to  the  square  of  the  intrinsic  carrier 
concentration.  The  introduction  of  R(x)  implies  the  transport 
velocity  transformation  [19],  often  used  in  modeling  of  quasi-neutral 
regions  through  the  factor,  Jp/u  - JpP0/(P-P0),  where  Jp/(P-PQ)  is 
the  effective  velocity  of  excess  carriers. 

Using  (3.7)  and  some  manipulation,  we  transform  (3.6)  into 


^ 1 Neff  2 

dx  “ D/Neff  ’ D R ’ (3'8) 

which  is  a first-order  nonlinear  differential  equation  of  the  Riccati 
form.  The  boundary  condition  for  (3.8)  is 


R(0)  - Rs  


1 JP(0)  1 P(0)-po(0)  Po(0) 


en^  u(0) 


u(0) 


n; 


Neff(O) 


(3.9) 


where  the  subscript  s refers  to  any  surface,  internal  or  external,  of 
the  solar  cell,  whose  surface  recombination  velocity  is  S.  The  ratio 
Rs  ” [S/Neff ]x_o  completely  characterizes  the  surface  for 
minority-carrier  recombination  in  steady-state.  Equation  (3.8) 
strictly  holds  only  in  the  quasi-neutral  regions  bounded  by  H/L 
junctions.  However,  because  of  the  previous  assumption  about  the 
carrier  transport  in  the  SCR  of  H/L  junctions,  R(x)  is  a continuous 
function  at  the  points  that  define  the  H/L  junction  in  our  model.  The 
continuity  of  R(x)  is  important  for  the  analysis  that  follows. 

Our  task  is  to  find  lower  bounds  for  R(x)  and  use  these  bounds 


to  determine  the  optimum  doping  profile  that  minimizes  R as  a function 
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of  Rs  and  distance  from  the  surface.  To  do  so,  a theorem  will  be 
presented  that  derives  lower  bounds  for  R(x)  as  a function  of  Rs  and 
another  bulk  parameter  Rb , which  depends  only  on  Neff  and  is  defined 
as 


Rb(x) 


Neff(x) 


( D(x)/r(x))V2 


D(x) 


Neff (x)L(x) ’ 


(3.10) 


Theorem  I 

R(5)  > min  [Rb(x),Rs],  0<  x <9  (3.11) 

where  9 is  the  thickness  of  the  region  under  study,  Fig.  3.1,  and  the 
prefix  min  stands  for  the  greatest  lower  bound  of  the  set  consisting 
of  all  values  of  Rb(x),  for  0 < x < 9,  plus  the  value  of  Rs . 

Proof:  Let  xm  define  the  point  in  the  closed  interval  [0,5] 

where  R(x)  assumes  its  minimum  value.  If  xm=0 , then  min[R(x) ] - Rs 
and  the  theorem  is  proved.  If  xm=5 , then  using  (3.8)  and  R(x)  > 0, 
we  have 


dR(x) 

lim  —3 <0  o 

x-0  dx 


1 Neff(5)  2 

Neff(0)r(5)'  D (9)  R (6)  - ° 


0 R(5)  > Rb(5)  > min[Rb(x)]  > min[Rb(x),  Rs],  0 < x < 9. 
If  0 < xm  < 9 , then 


d®-(xm)  3.  ^eff(xm)  2 

° Neff(xm)r(xm)  " D(xm)  R (xm)  =0  ° R(xm)  “ Rb(xm> ■ 


Thus  R( 5 ) > R(xm)  > min[Rb(x)]  > min[Rb(x),  Rs],  0 < x < 9,  Q.E.D. 
In  a similar  way  one  could  derive  an  upper  bound  for  R (5): 


R(5)  < max [Rb (x) , Rs] , 


0 < x < 9 . 


(3.12) 
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In  the  proof  of  the  above  theorem  we  implicitly  assumed  the  absence  of 
H/L  junctions.  However,  the  theorem  can  be  extended  to  cover  any 
number  of  such  points  by  using  the  continuity  of  R(x)  and  applying 

(3.11)  and  (3.12)  to  every  quasi  neutral  region  between  successive  H/L 
junction  points. 

Theorem  I states  that  if  a region  terminates  (x=0)  at  an  ohmic 
contact,  or  if  Rs  > min[R^(x)],  then  the  minimum  value  of 
^b(x)  ” (1/Neff) (D/r)^/^  between  0 and  9 imposes  a lower  bound  for 
R(0).  If  Rs  < min[Rb(x)],  the  lower  bound  becomes  Rs . If  9 
approaches  infinity,  then  theorem  I simplifies  to  R(9)  > min[R^(x)] 
since  the  condition  at  the  surface  has  no  effect  on  the  recombination 
current.  Theorem  I holds  for  any  drift  field  and  H/L  junction 
combination  that  constitutes  the  quasi-neutral  region  of  a solar  cell. 
In  this  sense,  the  smallest  value  that  R^  can  take  is  a fundamental 
property  of  the  semiconductor. 

We  now  introduce  notation  for  the  extrema  of  several  quantities 
that  will  be  defined  below.  Let  Y be  a parameter  that  depends  only  on 
the  doping  concentration  NDD.  Although  our  main  focus  in  this  work 
will  be  on  the  doping  concentration,  we  note  that  N - NDD  in  the 
quasi-neutral  region  under  study,  and,  for  brevity,  we  use  N in  the 
notation  that  follows  to  mean  the  doping  concentration  as  well  as  the 
majority  carrier  concentration.  Let  Ym  and  YM  be  the  minimum  and 
maximum  values  of  Y as  a function  of  the  doping  density.  When  the 
letter  f is  added  to  the  subscript,  we  mean  that  the  parameter  is 
limited  only  by  fundamental  mechanisms  (e.g.  Auger).  Superscripts  n 
and  p will  denote  n-type  and  p-type  material.  The  value  of  N 
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corresponding  to  a particular  value,  y,  of  the  parameter  Y is  denoted 
by  N(y).  In  this  notation,  R^m  is  the  minimum  value  of  Rb  and  is 
larger  than  or  equal  to  Rbmf,  which  for  silicon  is  limited  only  by 
Auger  recombination: 

Rbfm  “ t(DCA)1/2N/Neff]m>  (1.13) 

where  CA  is  the  Auger  coefficient.  The  quantity  N/Ne££  is  an 
increasing  function  of  N,  while  D is  a decreasing  one.  Although  there 
is  some  disagreement  in  the  literature  concerning  the  actual  values  of 
Neff  and  D as  functions  of  N for  N > 1018  cm'3,  the  consensual 
agreement  among  different  models  leads  to  the  conclusion  that  Rbmf 
occurs  in  the  vicinity  of  N - 1017  cm'3  for  both  n and  p-type 
material.  For  n-type  material  with  CAn  = 2.8xl0'31  cm6s'1  [20],  Rbmfn 
~ 1-7x10  cm4  s Similarly,  for  p-type  material  it  turns  out  that 
RbmfP  “ Rbmfn  « 1.7x10' 15  cm4  s'l. 

In  today's  silicon  cells,  where  the  lifetime  for  N < 1018  cm'3 
is  several  times  smaller  than  the  Auger  limited  value , Rbm  can  occur 
in  the  vicinity  of  N - 1018  cm"3  and  have  values  Rbmn  « RbmP  « 

3x10 cm4  s ' ^ or  more. 

Applying  theorem  I to  a solar  cell  having  Rs  > Rbmfn  ~ RbmfP  f°r 
the  surfaces  of  both  the  emitter  and  the  base  we  derive  the  following 
relation  for  the  saturation  current 

JQ  5:  en£  (Rbmfn  + RbmfP)  ~ 0-7x10' 13  A cm'2,  (3.14) 

holding  under  any  drift  or  back  surface  field.  In  (3.14)  and  from  now 
on,  we  take  en^  = 20  A cm'6s  (24  °C)  to  avoid  the  uncertainty 
presently  existing  about  effective  masses  and  their  temperature 
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dependence.  To  further  reduce  JQ,  Rs  must  be  reduced  below  the 
limiting  value  of  Rbmf  » 1.7x10  ^ cnA/s,  corresponding  to  a surface 
recombination  velocity  S — 1700  cm/sec  at  Ne££  = 10^®  cm'^  or  S = 17 
cm/s  at  Neff  - lO^  cm*^. 

The  previous  analysis  on  the  lower  bounds  of  R(x)  can  be 
extended  to  three-dimensional  quasi-neutral  regions  under  low 
injection.  This  is  done  by  replacing  Jp  in  (3.1)  and  (3.2)  with  the 
integrated  current  density  over  an  equipotential  surface  (potential  in 
the  sense  of  minority  carrier  quasi  Fermi-potential),  at  potential  Vp . 
In  this  case  Rb  in  (3.11)  can  be  replaced  by  its  limiting  value 

A[(l/Neffr)m/(Neff/D)M]1/2  < Rb , 

where  the  extrema  are  taken  over  an  equipotential  surface  having  area 
A.  In  the  general  case  these  surfaces  are  not  known.  Thus,  in  order 
to  derive  practical  bounds  for  Rb , the  extrema  can  be  taken  over  the 
entire  volume  of  the  region  under  study.  The  minimum  value  of  A can 
be  inferred  either  from  the  morphology  of  the  injecting  boundaries  or 
the  external  surfaces. 

Before  ending  this  section,  notice  that  despite  the  discrepancies 
in  the  published  data  on  the  band  gap -narrowing  and  minority- carrier 
diffusivity  in  heavily  doped  silicon,  the  parameter  Rb(x),  appearing 
in  theorem  I,  can  be  determined  relatively  accurately  by  experiment, 
at  any  doping  level.  This  can  be  done  by  measuring  the  recombination 
current  density,  eniRb,  in  a uniformly  doped  region  of  thickness  much 
greater  than  the  diffusion  length.  Similar  experiments  can  also 
determine  the  values  of  the  coefficients  of  the  basic  equation  (3.8), 
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1/Neffr  and  Neff/D,  as  was  recently  pointed  out  by  del  Alamo,  Swirhun 
and  Swanson  [21].  Because  of  the  absence  of  spatial  derivatives  of 
these  quantities,  we  may  consider  (3.8)  to  be  more  useful  than  its 
equivalent,  (3.6),  which  requires  knowledge  of  the  derivative  of  the 
experimentally  determined  parameter  D/Neff. 

3 • 3 Optimization  of  R(x,  Rc  ) 

Returning  to  (3.8),  we  observe  that  the  right  hand  side  is  a 
single-valued  and  continuous  function  of  x and  R in  a quasi-neutral 
region.  Thus,  through  every  point  in  the  (R,x)  plane  passes  only  one 
integral  curve  [22].  Therefore,  by  changing  only  the  recombination 
conditions  at  the  surface  from  Rgl  to  Rs2,  leaving  the  bulk  unchanged, 
we  have 


R(x,Rsl)  / R(x,Rs2)  o R(x,Rs1)  - R(x,Rs2)  / 0,  0 < x < 

where  R(x,Rs)  is  the  function  R(x)  corresponding  to  surface  value  of  R 
equal  to  Rs . Thus, 

sgn(Rsl  - Rs2)  - sgn[R(x,Rsl)  - R(x,Rs2)]  0 < x < (3.15) 

where  the  function  sgn  is  defined  by  the  equation  sgn(x)  = x/|x|  with 

x a non- zero  real  number. 

A direct  result  of  Equation  (3.15)  is  the  following  Lemma. 

Lemma  1:  If  R(x,Rs^)  and  R(x,Rs2)  are  the  normalized 

recombination  currents  of  the  same  region  having  different  boundary  R 
values,  Rsl  > Rs2,  then 


R(x,Rsl)  > R(x,Rs2)  , 


0 < x < R. 
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Using  the  same  arguments  as  in  theorem  I,  we  can  extend  the  above 
Lemma  to  a region  having  any  number  of  H/L  junctions.  This  Lemma 
leads  to  the  following  optimization  theorem. 

Theorem  II:  a)  Given  Rs  and  9,  the  doping  profile,  N0p(x,Rs), 

that  minimizes  R(0)  also  minimizes  R(x)  for  0 < x < 9.  b)  Given  Rs , 
the  optimum  doping  profile,  Nop(x,Rs),  is  the  one  that  at  any  x 
minimizes  dR/dx  or,  equivalently,  1/Neffr  - R2Neff/D. 

From  the  discussion  following  theorem  I we  know  that 
limx_KXJR0p(x,Rs)  — Rbm  a*^d  ■^nhi-KoN0p(x,Rs)  = ^(Rbm)  for  any  Rs , where 
R0p(x,Rs)  is  the  minimum  value  of  the  normalized  recombination  current 
as  a function  of  x and  Rs  and  N(Rbm)  is  the  doping  concentration 
corresponding  to  Rbm.  Also,  from  theorem  II  and  Equation  (3.8) 
written  in  the  form  dR/dx  - (1/Nefft)(l  - R2/Rb),it  follows  that  when 
dRop/^"*0-  then  Rop  "*  Rbm  and  vice  vetsa.  Therefore,  R (x,Rs)  is  a 
decreasing  function  of  x if  Rs  > Rbm,  an  increasing  function  of  x if 
Rs  < Rbm’  and>  for  large  x,  approaches  the  bulk-limited  value  Rbm 
(Rl§-  3.2)  . If  Rg  — Rbm>  then  ®-op(x,Rbm)  =*  Rbnj  = constant  and 
Nop(x,Rbm)  - N(Rbm).  Knowledge  of  Rop(x,Rs=0)  and  Rop (x , Rs=®) , 
combined  with  the  first  part  of  theorem  II  yields  Rop(x,Rs)  for  any 
Rs  ■ If  Rs  > Rbm>  the  function  Rop(x,Rs)  is  the  same  as  Rop(x,Rs=oo) 
shifted  to  the  left  until  it  intersects  the  R axis  at  Rs  (Fig.  3.2). 

In  more  analytical  terms  Rop(x,Rs)  - Rop (x+xs , Rs=®)  where  Rop(xs ,RS=») 
= Rs.  If  Rs  < Rbm,  the  shifting  applies  to  Rop(x,Rs=0).  To  derive  the 
optimum  doping  profiles  for  finite  values  of  Rs , we  can  apply  the  same 
translational  procedure  to  the  optimum  doping  profiles  N0p(x,Rs=«)  and 
^op ’ Rs  “0) , as  shown  in  Fig.  3.3. 
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Minimum  normalized  recombination  currents  for  R = n 
"s  = « and  their  relation  to  R (x,Rs)  for  any 


R 


Fig.  3.2 


paw 
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Fig.  3.3 


Optimum  doping  profiles  for  Rs  - 0 and  Rs  = ® and  their 
relation  to  Nop(x,Rs)  for  any  Rs . 
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The  second  part  of  theorem  II  provides  a way  to  obtain  Rop(x,Rs) 
and  Nop(x,Rs) . To  each  value  of  R corresponds  a doping  density, 
Nop(R),  which  minimizes  1/Neffr  - R2Neff/D  = f(N,R) . The  minimum  of 
f(N,R)  over  all  values  of  N at  a given  R will  be  a function  only  of  R: 
f(Nop(R),  R)  = R(R) • Thus,  the  differential  equation  for  Rop  is 

^op 

“ F(Rop}-  (3.16) 

Solution  of  (3.16)  provides  Rop(x),  while  Nop(x)  can  be  obtained  from 
the  function  Nop(R).  A practical  way  of  obtaining  Rop  is  to  use  a 
finite  element  method  and  divide  the  region  by  a sequence  of  points , 
xn,  into  short  segments,  each  of  which  can  be  assumed  uniformly  doped 
(Fig.  3.4).  We  obtain  the  optimum  doping,  Nop(Rs),  of  the  first 
segment  by  minimizing  f(N,Rs).  The  boundary  condition  for  the  second 
segment  is  Rb  = Rop(x^,Rs)  which  can  be  obtained  from  analytical 
expressions  available  for  uniform  regions.  Again  the  optimum  doping 
segment  2,  Nop(R1),  is  obtained  by  minimizing  (l/rNeff)  - RL2(Neff/D). 
Repeating  this  process  yields  Nop(Rn)  as  the  optimum  doping  of  the 
(n+1) th  segment.  From  now  on  we  will  implicitly  use  the  relation, 

N[(l/Neffr)m]  < N(Rbm)  < N[(Neff)/D)M]. 

This  inequality  is  realistic  for  silicon.  All  published  data  on 
band-gap  narrowing  and  mobility  indicate  that  Neff/D  is  an  increasing 
function  of  N,  at  least  for  N as  high  as  1019  cm-3  [4].  On  the  other 
hand,  several  different  models  for  the  lifetime  imply  that  1/Neffr 
decreases  as  N decreases  at  least  for  N as  low  as  1017  cm'3  [23,24], 
Thus  N[ (1/Neffr )m]  < [ (Neff )/D)M] . From  this  and  from  the  equality 
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1 Nop(Rs)  Nop(R.j)  ^ 

N0p(Rn) 

Rs  •*!  X,  ^'Xj  l 

*1  Ra 

"*]  'xn  * xn+1 

"n 

Fig.  3.4  Finite  element  method  for  obtaining  Nop(x,Rs)  and  Rop(x,Rs). 
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Rb  = [ (VNeffr)/(Neff/D)  ]1/2>  iC  follows  that  the  value  of  N that 
minimizes  Rb  lies  between  the  value  that  minimizes  1/Neffr  and  the  one 
that  maximizes  Neff/D. 

The  discussion  following  theorem  II  implies  that  it  is  important 

to  know  R0p(x,Rs-«) , ROp(x,Rs-0)  and  the  corresponding  functions 

NopCx.Rg-00) , Nop(x,Rs-0).  Dealing  first  with  the  case  of  an  ohmic 

contact,  we  see  that  if  Rs-«>,  then  minimizing  dR/dx  = 1/N effr  - 
2 

R Neff/°  is  equivalent  to  maximizing  Neff/D.  There  is  experimental 
evidence  that  Neff/D  is  an  increasing  function  of  N,  at  least  for 
n-type  silicon,  up  to  a doping  density  of  1020  cm'3  [4].  Therefore, 
in  order  to  minimize  the  recombination  rates  of  a region  terminated  at 
an  ohmic  contact,  the  doping  of  the  layer  adjacent  to  the  contact 
should  be  high,  approaching  the  solid  solubility  limit.  This  could 
account  for  the  passivation  effects  of  heavily  doped  polysilicon 
layers  between  single  crystal  silicon  and  a metal  contact.  The 
optimum  profile  N0p(x,Rs=«o)  is  a continuous  and  decreasing  function  of 
x since,  when  N drops,  1/Neffr  decreases  much  faster  than  does  Neff/D, 
provided  N > N(Rbm)  . As  explained  before,  Rop(x,Rs=«>)  is  also  a 
decreasing  function  of  x.  The  same  is  true  from  any  R(x,Rs=«) 
resulting  from  any  doping  profile  provided  the  doping  is  a decreasing 
function  of  x and  greater  than  N(Rbm).  This  is  a consequence  of 
theorem  I and  of  the  fact  that  Rb  is  an  increasing  function  of  N for  N 
> N(Rbm).  Thus,  for  a doping  profile  obeying  the  conditions  stated, 


dR 
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eff 
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This  observation  leads  to  the  conclusion  that,  to  minimize 
recombination  currents,  a diffusion  underneath  an  ohmic  contact  should 
be  made  deep  and  should  have  a surface  concentration  higher  than 

NOW  • 

Consider  now  a completely  passivated  surface.  Because  Rs  - 0 , 
(3.8)  implies  that  the  optimum  doping  density  at  the  surface  is 
apparently  the  density  that  minimizes  1/Neffr : Nop(Rs=0)  - 
N[ (l/NeffT )m] . Note,  however,  that  this  doping  concentration  will  not 
in  general  produce  Rs  - 0.  Rather,  this  will  be  some  other  doping 
concentration  Ns . Thus  the  optimum  doping  profile  near  the  surface 
will  consist  of  a thin  layer  of  concentration  Ns  followed  abruptly  by 
the  optimum  doping  concentration  stated  above.  There  is  a disagreement 
in  the  literature  on  the  value  of  Nop(Rs-0).  Rohatgi  and  Rai-Choudhury 
[23]  locate  this  doping  in  the  range,  10^  . cm*3.  por  sucb  an 

Nop(Rs=0)  the  optimum  doping  and  Rop(x,0)  are  derived  as  previously 
described.  As  shown  in  Fig.  3.3,  the  doping  density  increases  from 
N[ (1/Neffr )m]  to  N(Rbm).  This  optimum  doping  profile  as  an  increasing 
function  of  x creates  a field  that  pushes  minority  carriers  towards 
the  surface,  which  is  reasonable  because  near  the  passivated  surface 
the  carriers  recombine  at  a lower  rate  than  that  in  the  bulk.  The 
direction  of  the  field  remains  unchanged  as  long  as  Rs  < Rb  This 

should  be  contrasted  with  the  conventional  drift  field  that  pushes 
minority  carriers  away  from  a surface  having  a high  recombination 
rate,  such  as  an  ohmic  contact.  According  to  other  sources  [24], 

1/Neffr  is  a increasing  function  of  N,  at  least  for  doping  densities 
as  low  as  10^  cm'^.  Thus  Nop(Rs-0)  will  have  a small  value  yielding 
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a material  which  will  be  driven  into  high  injection.  In  this  sense, 
the  above  analysis,  about  NOp(Rs=0),  becomes  unrealistic  although  it 
is  still  valid  from  a mathematical  viewpoint.  Even  in  this  case, 
however,  the  doping  near  the  surface  should  be  kept  low  to  provide 
good  matching  between  the  passivated  surface  and  the  semiconductor 
bulk  characterized  by  Rbm.  Note  that  if  the  lifetime  is  limited  only 
by  radiative  and  Auger  processes,  then  1/Neffr  decreases  as  N 
decreases  and  saturates  at  [1/Neffr]mf  = Ra  for  N < 1016  cm'3.  Here  Ra 
is  the  radiative  recombination  constant  which  for  silicon  is  R„  = 

a 

2. 1x10' 15  cm3/s  [18]. 

In  Figs.  3.5  and  3.6  the  functions  Rop  and  Nop  are  shown, 
numerically  calculated  from  data  provided  by  Rohatgi  and  Rai-Choudhury 
for  n-type  silicon  [23].  From  these  figures  it  is  seen  that 
NOp(x,Rs-0)  stays  almost  constant  at  N[ (1/Neffr)m]  until  Rop(x,Rs=0) 
approaches  Rbm  and  then  it  sharply  increases  to  values  close  to 
N(Rbm).  The  same  behavior  of  Nop(x,Rs=0)  is  observed  for  many 
different  functional  dependencies  of  r on  N.  As  an  implication  of 
this,  a practical  way  of  obtaining  optimum  profiles,  when  Rs  < Rbm,  is 
to  keep  N constant  and  equal  to  N[ (l/Ne££T )m]  from  the  surface  to  a 
depth  where  R will  become  Rbm.  At  this  point  the  doping  density 
abruptly  changes  to  N(Rbm)  (H/L  junction)  guaranteeing  that  R will  not 
further  increase  beyond  Rbm.  In  the  case  of  an  ohmic  contact,  both 
functions  Rop  and  Nop  were  calculated  by  taking  Rs  as  the  ratio  of  the 
kinetic  limit  velocity,  S = 107  cm/s , and  Neff  = 1018  cm'3.  The 
values  of  Rop , Nop  corresponding  to  such  an  Rs  are  practically  the 
same  as  Rop(x,Rs=«>)  and  Nop(x,Rs=<»)  which  result  from  setting  Rs  = oo. 


(cm4 /sec) 
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Fig.  3.5  Numerical  results  for  the  minimum  value  of  R(x,Rs).  The 
inset  emphasizes  the  behavior  near  the  surface  for  Rs  = 


(cm-3) 
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Fig.  3.6  Numerical  results  for  the  optimum  doping  density.  The  inset 
emphasizes  the  behavior  near  the  surface  for  Rs  = <». 
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Note  that  in  Fig.  3.6  the  optimum  doping  adjacent  to  an  ohmic  contact, 

1 Q O 

which  maximizes  Neff/D,  is  approximately  1.5X10^^  cm  , according  to 
the  data  of  Ref.  23.  This  is  substantially  lower  than  the  value  of 
N0p(Rs-»)  ~ 10^®  cm'^,  indicated  by  del  Alamo  and  Swanson  [4]. 

It  is  of  interest  to  have  analytical  expressions  that  will 
constrain  both  ROp(x,Rs-0)  and  R0p(x,Rs=«>)  and  can  be  used  as 
closed-form  bounds  for  any  R(x,Rs).  Starting  with  Rs=0 , we  show  in 
Appendix  A that 

Rbm  tanh(x/lm)  < Rop(x,Rs-0)  < min[Rb (N]_)  tanhCx/LC^)  ) , 

Rbmtanh(x/L(N2))  ] - (3.18) 

where  lm  is  Rbm/(^-/^effT)m>  L(Nb)  is  the  diffusion  length  when  N = Nb 
- N[(l/Neffr)m]  and  N2  - N(Rbm) . 

The  lower  bound,  Rbmtanh(x/lm)  is  absolute  in  the  sense  that  the 
saturation  current,  in  low  injection,  of  an  one -dimensional  region 
having  thickness  x is  greater  than  or  equal  to  en^Rbmtanh(x/lm)  for 
any  Rs . The  length  lm  indicates  how  fast  ROp(x,Rs=0)  saturates  to  its 
bulk-limited  value  Rbm.  Using  data  by  Rohatgi  and  Rai-Choudhury  [23] 
for  n-type  silicon,  Rbm  - 3.6x10'^  cnA/s>  [1/Neffr]m  “ 5.1x10'^ 

O 

cm-ys  so  that  lm  = 70/xm.  A solar  cell  base  made  of  such  a material  and 
having  a passivated  back  surface  requires  a thickness  much  less  than 
70  fim  to  yield  an  optimized  saturation  current  less  than  en^Rbm  « 
0.7x10'^  A/cm. 

For  an  ohmic  contact  R0p(x,Rs=>»)  obeys  the  relation 

Rbm  coth(x/l'm)  < Rop(x,Rs=oo)  < min[Rb(N3)coth(x/L(N3) ) , 

Rbmcoth(x/L(N2)) ] > 


(3.19) 
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where  l'm  = [ <Neff/D)M  R^]’1  and  N3  = N[ (Neff/D)M] . The  length  l'm  is 
the  least  thickness  required  to  isolate  an  ohmic  contact  from  the 
bulk  and  limit  R0p(x,Rs=<»)  to  near  R^m-  Typical  values  of  lm'  vary 
from  3 to  10  pm.  The  bounds  derived  above  can  accommodate  any  Rs  by 
using  the  translational  relation  between  ROp(x,Rs=0),  R0p(x,Rs=<»)  and 
a given  Rop(x,Rs). 

3.4  Relation  of  the  Present  Model  with  the  Transparent  and 
Quasi -Equilibrium  Models  and  Optimization  of  the  Short-Circuit  Current 

Among  the  approximate  models  used  in  the  literature  to  derive 
analytical  expressions  for  JQ(x)  - en^R(x) , two  particular  models  are 
related  with  the  previous  analysis  and  can  themselves  yield  bounds  for 
the  recombination  current  and  internal  quantum  efficiency. 

The  transparent  model  [25],  based  on  the  assumption  that  the 
minority- current  density  is  independent  of  position,  can  be  derived 
from  the  basic  equation  (3.8),  dR/dx  - l/(Neffr)  - R2(Neff/D),  by 

O 

setting  R^Ngff/D  » 1/Neffr  (negligible  bulk  recombination) : 
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R(x)-  — . (3.20) 

Jo<Neff/D>  dx'+l/Rs 


The  quasi -equilibrium  model  [26],  based  on  the  assumption  that  the 
minority-carrier  quasi-Fermi  potential  Vp  is  constant,  can  be  derived 
from  (3.8)  by  setting  1/Neffr  » R2Neff/D: 
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(3.21) 


This  rederivation  of  both  models  demonstrates  that  the 
h^^ispar sncy  assumption  underestimates  R while  the  quasi-equilibrium 
assumption  overestimates  R.  Therefore,  the  expressions  for  R derived 
in  the  quasi - transparent  and  quasi-equilibrium  model  can  be  used  as 
lower  and  upper  bound  for  R,  respectively: 
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(3.22) 


From  the  discussion  in  the  previous  section  it  follows  that  the 
transparent  model  is  accurate  only  if  x « l'm  and  Rs  > Rbm.  This 
model  predicts  an  optimum  doping  Nop  = N[ (Neff/D)M) ] , which  maximizes 
Neff/D<  independently  of  Rs  and  x.  This  was  shown  to  be  true  for  thin 
regions  having  ohmic  contacts.  On  the  other  hand  the  quasi -equilibrium 
model  is  accurate  if  x « lm  and  Rs  < Rbm.  This  model  predicts  an 
optimum  doping  Nop  = N[ (1/Neffr )m] , which  was  shown  to  be  the  case  for 
a thin  region  having  a passivated  surface.  Therefore,  these  two 
models  converge  to  the  accurate  model  developed  in  the  previous 
section  for  limiting  cases. 

The  transparent  model  can  be  used  to  derive  upper  bounds  for  the 
internal  quantum  efficiency  of  a region  in  low  injection.  This 
derivation  is  based  on  the  reciprocity  theorem  of  Chapter  two.  Let 
ms (x)  he  the  saturation  value  of  the  normalized  excess  minority- 
carrier  density.  Because  the  transparent  model  underestimates  the 
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current,  (3.1)  implies  that  the  ms(x)  predicted  by  the  transparent 
model,  mstr(x) , is  higher  than  the  actual  one. 

Thus  us ing  (2.9) 

n8 

Isc  = e cp(x)  G(x)  dx  < e mstr(x)  G(x)  dx  < 

J0  J0 


JoNeff/D)dx'  + 1/Rs 

0 Jo(Neff/D)dx  + 1/Rs 


G(x)  dx, 


(3.23) 


where  G(x)  is  the  generation  rate.  Relation  (3.23)  in  the  form  of  an 
equation,  instead  of  an  inequality,  was  used  by  del  Alamo  and  Swanson 
[4]  in  deriving  the  internal  quantum  efficiency  of  a transparent 
emitter.  Here  we  have  shown  that  the  same  expression  serves  as  an 
upper  bound  for  the  internal  quantum  efficiency  of  any  region. 

Correspondingly,  a first-order  quasi-equilibrium  model  can  yield 
a lower  bound  for  the  short-circuit  current: 


Isc  = e Jocp(x)  dX 
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From  the  reciprocity  theorem,  the  zeroth  order  quasi  equilibrium 
model  predicts  100%  internal  quantum  efficiency,  which  is  the  same  as 
saying  that  in  a region  having  a flat  Vp  from  the  inj  ecting  boundary 
to  the  surface,  all  the  photogenerated  minority  carriers  are  collected 
by  the  junction.  Therefore,  the  flatness  of  Vp  (that  is,  dVp/dx  = 0) 
is  the  criterion  for  optimizing  a region  with  respect  to  current 
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collection,  under  low  injection  conditions.  Since  lower  recombination 
currents  help  keep  the  Vp  flat,  one  could  expect  that  a region  having 
minimized  recombination  losses,  R(x,Rs)  - Rop(x,Rs),  would  also  have  a 
high  internal  quantum  efficiency.  However,  optimization  of  the 
short-circuit  current  for  a given  excitation  is  a separate  problem 
which  leads  to  a solution  different  from  the  one  that  optimizes 
R(x,Rs) . 

The  optimum  doping  concentration  for  maximum  current  collection 
depends  on  the  optical  generation  profile.  If  the  generation  is 
confined  between  9 (injecting  boundary)  and  a point  x' , then,  based  on 
the  previous  discussion,  the  optimum  doping  between  the  surface  (x=0, 
Fig.  3.1)  and  x'  is  the  same  as  Nop(x,Rs),  0 < x < x' , the  doping 
minimizing  R(x',Rg).  For  the  case  where  the  generation  function  is  a 
delta  function  applied  at  x' , G(x)  — 5(x-x'),  the  optimum  doping 
between  x'  and  9 is  zero  (intrinsic  material)  if  both  Neff/D  and 
l/Ngffr  are  increasing  functions  of  N.  This  conclusion  is  based  on 
the  reciprocity  theorem  previously  mentioned.  It  becomes  unrealistic 
because  of  the  high  injection  condition  that  would  prevail  in  such  an 
intrinsic  material.  Even  in  high  injection,  however,  low  doping 
concentration  would  increase  the  collection  probability  at  x*  because 
of  long  diffusion  lengths.  This  is  all  the  more  true  when  the  field 
created  by  the  redistribution  of  the  electron-hole  plasma  helps  in 
separating  the  photogenerated  carriers.  Such  a low  doping 
concentration  could  be  expected  to  be  optimal  if  G(x)  is  uniform 
between  9 and  x.  In  a base  region  subjected  to  solar  illumination  the 
generation  near  the  injection  boundary  (x=0)  is  much  higher  than  in 
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the  bulk.  To  maximize  the  current,  a drift  field  towards  x = 9 could 
become  necessary.  Such  a field,  though,  could  conflict  with  the 
optimization  of  R(x,Rs),  as  presented  in  the  previous  section.  An 
optimized  combination  of  these  two  opposing  tendencies  will  be 
presented  in  the  last  section  of  this  chapter. 

In  the  literature  [18]  the  short-circuit  current  is  usually 
believed  to  be  an  increasing  function  of  the  base  thickness.  This  is 
true  for  a hypothetical  base  having  unity  collection  probability  at 
any  point.  It  is  also  true  for  any  actual  uniformly  doped  base  having 
Rs  - ^b-  However,  for  a base  with  a passivated  back  contact,  there  is 
an  optimum  thickness  that  maximizes  the  current  collection.  Let  a be 
the  absorption  coefficient  of  monochromatic  light,  L the  diffusion 
length  of  the  uniformly  doped  base  and  Rs-0 ; then  for  one  optical  pass 
this  optimum  thickness,  0op,  is  given,  within  a good  accuracy,  by  the 
approximate  expression 

^op’  1 2 

- Ln  (2.25) 

L 1 - La  La  + (La)^ 

This  optimum  thickness  also  exists  for  typical  AM  spectral  conditions 
and  is  function  only  of  L [27].  Figure  3.7  shows  the  ratio  0op/L  for 
an  AM2  spectrum  and  a uniform  base  having  Rs  - 0.  For  L in  a range  of 
practical  interest,  0op/L  varies  from  1 to  0.85  for  one  optical  pass, 
and  from  0.6  to  0.8  for  two  optical  passes.  In  the  same  figure  the 
maximum  short  circuit  current,  (ISC)M,  at  the  optimum  thickness  is 
shown,  normalized  with  respect  to  the  short-circuit  current,  (Isc^un* 
of  a uniform  base  having  infinite  thickness  and  the  same  doping.  As 
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uV*i)/w(s*i) 


Fig.  3.7  The  ratio  0Op'/L  and  the  optimized  optimum  current  as  functions  of 
the  diffusion  lenght  L for  AM2  light  and  for  R_  - 0. 
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Fig.  3.7  illustrates,  use  of  the  optimum  thicknesses  improves  the 
short  circuit  current  from  4 to  2%  for  two  optical  passes,  and  from  1 
to  0.4%  for  one  optical  pass,  this  improvement  being  a slowly 
decreasing  function  of  the  diffusion  length.  These  are  marginal 
improvements  considering  the  fact  that  they  correspond  to  completely 
passivated  surfaces  (Rs  - 0) . Texturizing  the  front  surface  is  a more 
effective  way  of  improving  Isc  as  has  been  pointed  out  by  the  recent 
developments  in  high  efficiency  Si  solar  cells  [28,29],  The  design 
implications  of  the  analysis  presented  in  the  last  two  sections  are 
discussed  in  the  next  section. 

3 . 5 Optimum  Design 

As  noted  previously,  if  on  both  the  front  and  back  surface  of 
the  cell  Rs  - [S/NDeff  or  S/NAeff]  > Rbmn,  RbmP,  then  the  saturation 
current  of  the  cell  will  be  more  than  JQ  = en^  (Rbmfn  + RbmfP)  = 

1 O o 

0.7x10  A cm  , where  this  value  of  JQ  and  a short-circuit  current 
of  36-40  mA  correspond  to  an  open-circuit  voltage  Voc  w 700  mV.  In 
order  for  the  efficiency  to  approach  theoretical  limits  [30],  Voc 
should  be  no  less  than  750  mV.  Such  an  increase  in  Voc  would  require 
values  of  Rs  much  less  than  Rbmf  and  thicknesses  substantially  less 
than  the  diffusion  lengths.  This  can  be  done  in  a shallow  diffused 
emitter  by  growing  quality  thermal  oxides  on  the  uncovered  regions. 

For  such  an  emitter,  the  base  generation  and  recombination  determine 
both  Isc  and  Voc. 

For  a base  region,  theorem  I implies  that  the  conventional  back- 
surface-  field  system  (terminating  in  an  ohmic  contact)  cannot  reduce 
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the  base  recombination  saturation  current,  eni2R(5),  below  eni2Rbm. 
This  condition  holds  also  for  a thick  base  region,  by  which  we  mean 
thickness  0 > lm,  independently  of  the  doping  profile  and  of  the  back 
surface  recombination  velocity.  Thus,  to  reduce  JQ  below  its 
bulk-limited  value,  J0  « eni2Rbmf  * 3.4xl0"14  A/cm2,  Rs  should  be  less 
than  Rbmf  « 1.7x10" 15  cm4/s  and  the  thickness  0 should  be  less  than  1 
= Rbm/(1/Neffr)m-  1116  commonly  accepted  principle  that  thinner  bases 
result  in  lower  recombination  rates  applies  to  an  optimized  cell  only 
if  Rs  < Rbm.  Otherwise,  the  optimum  R(x,Rs)  is  a decreasing  function 
of  x,  which  implies  that  the  base  should  be  made  thick  to  increase 
both  Voc  and  Isc-. 

As  noted  previously,  for  bases  with  well  passivated  back 
surfaces,  optimization  of  the  saturation  current  requires  a drift 
field  pushing  minority  carriers  to  the  surface.  This  field,  however, 
could  deteriorate  the  internal  quantum  efficiency.  Therefore,  it  is 
not  recommended  in  the  emitter  region  having  high  generation  rates,  or 
near  the  junction  edge  of  the  base.  A good  choice  for  the  doping 
concentration  of  a base  region  having  0 < lm  and  a well  passivated 
back  surface  (Rs  « Rbm)  follows  from  our  previous  considerations. 

This  doping  concentration,  N[ (1/Neffr )m] , is  the  one  that  minimizes 
(1/Neffr).  Such  a doping  concentration  would  minimize  the 
recombination  current  and  would  probably  be  low  enough  to  result  in  a 
long  diffusion  length  and  thus  high  internal  quantum  efficiency. 

For  regions  on  which  an  electrical  contact  should  be  made,  three 
solutions  are  applicable:  ohmic  contact,  polysilicon  (or  SIPOS) 
contact,  or  a thin  tunneling  oxide  between  metal  and  silicon.  If  an 
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ohmic  contact  is  chosen,  the  diffusion  under  the  metal  should  be  made 
deep  and  the  surface  concentration  high  to  maximize  the  ratio  Neff/D. 
As  mentioned  earlier,  the  junction  depth  needed  to  reduce  R(x,Rs=«) 
near  Rbm  is  approximately  [ <Neff/D)MRbm]  ' 1 - 3 - 10  /jm.  This  is  in 
agreement  with  numerical  results  of  del  Alamo  and  Swanson  on 
Gaussian-profile  emitters  [4].  The  deep  diffusion  under  the  ohmic 
contact  provides  the  added  benefit  of  increasing  the  reliability  by 
decreasing  the  possibility  of  metal  spikes  penetrating  into  the  base 
[31].  Moreover,  the  high  doping  density  near  the  silicon-metal 
interface  would  minimize  contact  resistance.  In  contrast  to  the 
doping  profile  under  the  metal  contact,  the  doping  in  the  uncovered 
and  passivated  emitter  regions  should  be  lighter  and  the  junction 
shallow  to  reduce  recombination  losses  and  maximize  current 
collection.  Series  resistance  constraints  place  a lower  limit  on 
emitter  doping,  so  that  the  optimum  doping  of  the  uncovered  emitter 
lies  in  the  range  10^  - 10^  cm  3.  Such  an  optimized  emitter  having 
step  discontinuities  in  doping  density  and  junction  thickness  is  shown 
in  Fig.  3.8. 

Polysilicon  contacts  have  shown  small  values  of  Rs (=1 . 5xl0'15 
cm4/s)  [32].  This  value  is  almost  the  same  as  Rbmf  and  smaller  than 
its  present  values.  Thus,  as  explained  earlier,  this  contact  should 
be  placed  as  close  to  the  junction  as  possible.  The  other  option, 
thin  tunneling  oxides,  has  also  been  used  to  decrease  surface 
recombination  [32].  Since  there  are  no  data  on  the  ratio  S/Neff 
characterizing  this  contact,  its  use  will  not  be  further  discussed 
here.  A grid  contact  and  a passivating  oxide  can  also  be  incorporated 
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~5  pm 


Fig.  3.8  Step  geometry  of  an  optimized  emitter  having  doping  density 
discontinuities  at  the  boundaries  between  the  metal-covered 
and  uncovered  regions  and  deeper  junctions  under  the  ohmic 
contacts . 
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at  the  back  contact  of  a thin  base  [34]  to  reduce  the  recombination 
below  the  bulk  limited  values.  Again,  a deep  diffusion  or  a 
polysilicon  contact  should  isolate  the  metal  at  the  contact  regions 
from  the  bulk  in  order  to  minimize  the  overall  effective 
recombination  velocity  at  the  back  surface. 

High  collection  efficiency  demands  that  care  be  taken  to  assure 
that  the  minority-carrier  quasi-Fermi  potential  Vp(x)  is  flat  in 
regions  having  high  generation  rates.  By  flat  we  mean  that  the 
variation  across  the  region  must  be  small  relative  to  kT/e.  This 
condition  can  be  satisfied  for  shallow  emitter  regions  having  doping 
concentrations  in  the  range,  l-5xl0^^/cm^ . As  an  example,  we  consider 
n-type  emitters  having  thicknesses  —0.2  pm,  constant  doping 
concentration  NDD  « 5xl019  cm'3  and  Rs=3xl0'15  cm4/s.  Taking  Neff/D  - 


the  quantum  efficiency  lower  bounds  derived  in  the  previous  section, 
we  find  that  the  collection  probability  at  the  surface  is  higher  than 
90%,  yielding  an  overall  internal  efficiency  for  the  emitter  very 
close  to  unity.  If  the  doping  concentration  were  103®  cm'3,  the 
collection  probability  at  the  surface  would  drop  to  nearly  50%.  In 
an  optimized  emitter  the  quasi -equilibrium  condition  (dVp/dx  » 0) 
should  be  combined  with  the  requirement  that  its  recombination 
current,  Joe,  is  much  less  than  the  base  recombination  current,  J^: 


5x10-*-^  cm'3  s, 


t — (2.8x10  1 s,  Ne££  — 1018  cm'3  and  applying 


where  Xj  is  the  thickness  of  the  quasi -neutral  emitter  not  under 
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metal.  This  relation  implies  that,  in  order  to  reduce  the  total  cell 
recombination  current  below  eni^K-bmf>  an  emitter  in  quasi -equilibrium 
having  a thickness  of  0.2  nm  and  Neff  = 1018/cm3  should  have  N < 1019 
cm"3  and  Rs  < 10 "•'■3  cm^/s . 

The  requirement  for  low  emitter  doping  is  not  necessarily- 
satisfied  by  assuring  that  the  average  value  of  or  NDD  is  low.  In 
a quasi -equilibrium  emitter  having  an  average  value  of  NDD  or  of 
10  cm'J  and  a surface  concentration  in  excess  of  1020  cm'3,  the 
recombination  at  the  thin  highly  doped  surface  layer  can  be  much 
higher  than  the  recombination  rate  in  the  rest  of  the  emitter  because 
of  the  quadratic  dependence  of  Auger  lifetime  on  majority  carrier 
concentration.  A practical  solution  for  a diffused  or  implanted 
emitter  would  be  a rather  square  profile  with  surface  concentration  no 
more  than  10^9  cm"3. 

This  doping  concentration  of  the  uncovered  emitter  is  consistent 
with  the  step  geometry  of  the  deep  diffusion  under  the  metal  shown  in 
Fig.  3.8.  If  the  conventional  planar  emitter  geometry  were 
incorporated,  an  emitter  doping  concentration  of  the  order  of  lcA9 
cm"  could  result  in  a rectifying  metal  contact.  On  the  other  hand,  a 
surface  concentration  in  the  uncovered  regions  in  excess  of  102®  cm"3 
could  make  it  hard  to  suppress  Rs  below  10 cm^/s  even  if  quality 
thermal  oxides  were  grown. 

An  indication  of  whether  the  emitter  satisfies  the  quasi- 
equilibrium condition  is  offered  by  the  limit  of  the  internal  quantum 
efficiency  at  low  wavelengths.  This  is  a result  of  the  reciprocity 
theorem  presented  in  the  previous  section.  Recent  data  on  high 
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efficiency  cells  indicate  that  this  limit  is  around  60-70%  even  with 
passivating  oxides  present  [35],  The  losses  thus  implied  could  result 
from  the  high  surface  impurity  concentration  characterizing  the 
emitters  of  these  cells.  This  observation  points  to  the  direction  of 
low  doping  in  the  uncovered  emitter  regions. 

As  was  mentioned  in  the  previous  section,  for  a base  with  a 
passivated  back  surface  there  is  an  optimum  thickness,  0O£,  that 
maximizes  the  collected  current.  Let  0op  be  the  thickness  that 
maximizes  the  product  VocIsc.  If  the  emitter  recombination  exceeds 
that  of  the  base,  then  0Op  — 0Op.  This  is  the  case  because  when 
5x10  ^ < JQ  < 5x10"^  A/cm^  and  Isc  = 36-40  mA/cm2,  then  the  product 
vocxIsc  is  about  27  times  more  sensitive  to  Isc  than  it  is  to  J0  = JQe 
+ ^ob  ’ anb  tbe  sensitivity  of  Job  to  a change  89  in  the  uniform  base 
thickness  is  at  most  89/6%.  However,  in  a cell  with  an  optimized 
emitter,  where  the  base  determines  both  JQ  and  Isc,  0op  could  be 
significantly  less  than  0op' . Figure  3.9  presents  the  ratio  0op/L  for 
a uniform  base  having  Rs  =■  0 and  two  different  values  of  Rb,  at  AM2 
spectral  conditions.  As  shown,  this  ratio  is  a decreasing  function  of 
L,  the  number  of  optical  passes  and  it  is  relatively  insensitive  to 
Rb-  As  expected  0Op  < 8Q p.  Figure  3.10  illustrates  the  maximum 
value,  (VocIsc)M,  that  the  product  VocIsc  assumes  at  0Qp.  This 
maximum  is  shown  normalized  with  respect  to  the  product,  (Vn„I  ) 
of  a uniform  and  infinitely  thick  base  having  the  same  doping.  As  can 
be  seen,  for  L > 70  /im  the  normalized  (VocIsc)j^  is  almost  independent 
of  L and  Rb ; in  other  words,  it  is  almost  independent  of  the  doping 
density.  This  observation  leads  to  the  conclusion  that  the  problem  of 
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Fig.  3.9  The  ratio  0op  as  a function  of  L for  two  values  of  Rh  and  for 
Rs  =0  (AM2  light). 
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The  normalized  optimum  value  of  the  product  VQC  Isc  as  a function  of 
L for  two  values  of  Rg  and  for  Rs  = 0 and  AM2  light. 
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finding  the  optimum  doping  and  thickness  which  will  maximize  VocIsc, 
for  a uniform  base,  consists  of  two  autonomous  problems.  First, 
optimum  doping  is  determined  by  requiring  that  V0CISC  be  maximized  for 
a uniform  and  infinite  base,  and  then  the  optimum  thickness  can  be 
inferred  from  Fig.  3.9.  Figure  3.10  also  shows  that  the  presence  of  a 
back  reflector  and  a passivated  back  contact  will  improve  the 
efficiency  by  about  6%  (or  at  least  1%  absolute)  if  the  right 
thickness,  0Op,  is  chosen.  The  data  by  Rohatgi  and  Rai-Choudury  [23] 
for  n-type  silicon  imply  an  optimum  doping,  2.3x10^  cm'^  and  an 
optimum  thickness  0Op  » 150  jum  for  one  optical  pass,  or  9 ~ 100  /an  for 
two  optical  passes. 

So  far  in  this  analysis  the  derivation  of  the  limiting 
recombination  losses  was  restricted  to  one -dimensional  regions.  In  the 
previous  chapter,  however,  we  pointed  out  that  in  shallow  junction 
structures,  especially  in  those  having  passivated  surface,  the  base 
recombination  current  per  unit  area  is  about  the  same  as  the 
recombination  current  density  in  the  one -dimensional  counterpart. 
Therefore,  the  conclusions  regarding  minimization  of  recombination 
losses  also  hold  for  three-dimensional  solar  cells,  like  the  ones  with 
texturized  front  surface  to  enhance  light  absorption  [28,  29]. 

3 . 6 Concluding  Remarks 

The  optimization  theory  presented  was  based  on  three  parameters 
that  determine  the  ultimate  performance  of  the  cell.  These  parameters 
are  Rs , 6,  and  1/Neffr , the  first  two  of  which  are  characterized  by 
two  limiting  values,  Rbm  and  lm.  It  was  shown  that  in  order  to  reduce 
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the  saturation  current  below  en^^R^m,  Rs  should  be  much  less  than  R^ 
and  9 less  than  lm.  If  these  conditions  are  satisfied,  the  optimum 
doping  lies  in  the  vicinity  of  N[ (1/Neffr )m] . The  fundamentally 
minimum  value  of  the  parameter  1/Neffr  sets  the  upper  bound  on 
efficiency.  The  dependence  of  the  fundamental  limit  minority- carrier 
lifetime  on  the  doping  concentration  is  still  debated  in  the 
literature.  For  N > 10^  cm'\  this  limit  is  probably  determined  by 
Auger  recombination.  For  lighter  doping,  however,  unavoidable  defects 
induced  by  processing  and  the  addition  of  impurities  can  lower  the 
lifetime  well  below  the  value  limited  only  by  Auger  and  radiative 
mechanisms  [36],  Of  the  same  importance  is  the  lowest  possible  value 
of  the  parameter  Rs  - S/Neff  characterizing  the  surface.  Further 
research  is  needed  to  investigate  the  fundamental  mechanisms,  if  any, 
that  limit  the  value  of  this  parameter. 


CHAPTER  FOUR 

LIFETIME  AND  DIFFUSIVITY  DETERMINATION  FROM  FREQUENCY -DOMAIN 

TRANSIENT  ANALYSIS 

4 . 1 Introduction 

In  this  chapter  we  propose  a new  method  to  measure  the  minority- 
carrier  lifetime  and  diffusion  coefficient  in  the  base  of  a solar 
cell.  We  also  point  out  that  several  time-domain  experiments  used  to 
extract  device  parameters  like  lifetime,  diffusion  coefficient  and 
recombination  velocity,  can  be  viewed  from  a different  than  the 
conventional  perspective.  In  methods  such  as  the  electrical  short- 
circuit  current  decay  [37]  and  photocurrent  decay  [38]  the  lifetime  is 
determined  by  waiting  until  the  response  drops  well  below  its  initial 
value  and  then  observing  the  almost  exponential  decay  of  the  signal. 

In  the  Haynes -Shockley  experiment  the  carrier  lifetime  and  diffusivity 
are  determined  by  monitoring  the  time  dependence  of  the  peak  value  and 
the  width  of  the  minority-carrier  pocket  injected  by  light.  Here  we 
demonstrate  that  in  these  experiments  we  can  exploit  the  full  response 
instead  of  focusing  on  a particular  part  or  point  of  the  curve.  This 
is  done  by  taking  the  Fourier  transform  of  the  pulse  response  of  the 
device  and  then  comparing  it  to  closed  form  expressions  (transfer 
functions)  relating  input  and  output  in  the  frequency  domain.  By 
fitting  the  entire  curve  for  the  amplitude  and  the  phase  of  the 
Fourier  transform  with  available  analytical  expressions  we  increase 
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the  accuracy  and  reliability  of  the  method.  The  next  two  sections 
demonstrate  the  technique  applied  on  the  pulsed  photocurrent  and  a 
variation  of  the  Haynes -Shockley  experiment.  Throughout  this  work  we 
assume  low- inj ection  conditions  so  that  input  and  output  are  linearly 
related  and  Fourier  transform  can  be  applied  to  derive  the  frequency 
response  of  the  device. 


4.2  Analytical  Formulation 

Consider  a one -dimensional  diode  where  the  quasi-neutral  base  is 
much  longer  compared  with  emitter  and  space-charge  region  combined  and 
is  excited  by  uniformly  absorbed  light  pulses.  As  shown  in  Appendix  B, 
the  short-circuit  current  in  the  frequency  domain  is  Jsc(w)  = e G(w) 
F(w) , where  G(u)  is  the  spatially  uniform  generation  rate  in  the 
frequency  domain,  and  F(w)  is  given  by 


F(w) 


[sinh(W/L')+(SL,/D)cosh(W/L,)-SL,/D] 
cosh(W/L' )+(SL'/D)sinh(W/L' ) 


(4.1) 


In  (4.1)  L'  is  the  complex  diffusion  length  L'=L/(l+jcor)^/^ , L is  the 
static  diffusion  length  L=(D r)^/^,  D is  the  diffusion  coefficient,  r 
is  the  lifetime,  S is  the  back  surface  recombination  velocity,  W is 
the  base  thickness  and  e is  the  elementary  charge.  Therefore,  when  a 
very  narrow  pulse  of  uniformly  absorbed  light  excites  the  cell 
(G(t)=S(t)  or  G(w)=constant) , the  Fourier  transform  of  the  short- 
circuit  pulse  response  current  Jsc(t)  is  the  same  function  of  the 
frequency  as  F(a>)  in  (4.1).  If  G(t)  is  not  a delta  function,  then  to 
obtain  F(w)  we  have  to  normalize  the  Fourier  transform  of  the  current 
by  dividing  Jsc(w)  by  the  Fourier  transform  of  the  input  G(w).  For  a 


85 


long  base  diode,  W/L  > 4,  Jsc(w)  is  proportional  to  (1+j ur)'^  and 
the  plot  of  F(w)  in  a log-log  scale  is  frequency- independent  for  <j  < 
1/r , a straight  line  with  a slope  of  -1/2  for  w » 1/r  and  the 
intercept  of  the  two  lines  is  the  corner  frequency  which  depends  only 
on  the  lifetime,  uc-l/r.  If  the  base  is  not  long  enough,  additional 
measurements  are  needed  to  extract  r,  D and  S.  However,  for  W/L  >1.5 
the  intercept  depends  primarily  on  the  lifetime  and  can  still  be  used 
for  a first  approximation  of  r. 

An  alternative  experiment  is  to  excite  the  device  using  strongly 
absorbed  light  pulses  striking  from  the  back  (after  opening  a window 
in  the  back  contact  ) and  then  monitor  the  current  collected  by  the 
shallow  front  emitter  junction.  This  is  a variation  of  the  well  known 
Haynes -Shockley  experiment  the  difference  being  that  the  minority 
carriers  move  vertically  and  in  one  dimension,  instead  of  laterally. 

If  G'  is  the  generation  rate,  per  unit  area  at  the  back  surface,  then, 
as  shown  in  Appendix  B,  the  solution  of  the  continuity  equation  is 
Jsc(w)  - e G'(w)  F' (w) , where 


F'(«)  - 


cosh(W/L' ) + (SL'/D)  sinh(W/L' ) 


(4.2) 


From  (4.2)  we  see  that  when  exp(-2W/L)  « 1 and  S < D/L  then 
F' (w)  is  proportional  to  exp(-W/L').  Thus,  if  we  plot  the  logarithm 
amplitude  of  F' (w)  versus  for  w » 2/r  we  obtain  a straight  line 

with  a slope  equal  to  -W(2D)*^/^  and  an  intercept  with  the  low- 
frequency  limit  at  (u>£  )^/^  = (2/r)-*-/^  . The  same  is  true  also  for 
the  slope  of  the  phase  plotted  versus  . Note  that  D and  r are 

independently  calculated  by  observing  different  portions  of  the  plot 


86 


of  F' (w)  vs  «.  For  higher  values  of  S the  slope  and  the  intercept  will 
slightly  change,  but  they  still  depend  primarily  on  D and  r and  can 
still  be  used  for  a first  estimate  of  D and  r. 

For  a more  accurate  determination  of  r,  D and  S one  can  combine 
the  uniformly  absorbed  light  and  strongly  absorbed  light  measurements. 
First,  from  the  corner  frequency  of  the  amplitude  of  Jsc(w),  for  the 
uniformly  absorbed  light,  we  derive  an  approximate  value  for  r.  Then, 
from  the  slope  of  the  amplitude  and  the  phase  of  Jsc(w),  for  the 
strongly  absorbed  light,  we  derive  a first  approximation  for  D.  To 
derive  S we  can  measure  the  quantum  efficiency  of  the  device  for 
strongly  absorbed  light  shining  from  the  back  surface.  This  quantum 
efficiency  is  a sensitive  function  of  S and  is  given  by  F' (0) . After 
having  made  a first  approximation  for  r,  D,  and  S,  a more  accurate 
determination  can  be  obtained  by  fitting  Jsc(w)  for  the  uniformly  and 
the  strongly  absorbed  light  with  the  analytical  expressions  for  F(w) 
and  F' (w)  as  given  by  (4.1)  and  (4.2). 

4 . 3 Demonstration  of  the  Method  and  Discussion 

A p-type  base  Si  cell  was  mounted,  with  the  emitter  facing  down, 
on  a TO- 5 header  after  its  back  contact  was  removed  except  for  a small 
dot  near  the  edge  to  bond  a wire.  The  cell  was  cut  from  a 2 x 2 cm^ 
cell  with  base  resistivity  of  0.16  ohm-cm.  The  base  doping,  as 
measured  by  the  capacitance  method,  was  1.22  ± 0. 04x10-'- ^ cm"^.  The 
base  thickness  was  309  ± 2 urn.  The  short-circuit  current  in  the  time 
domain  was  measured  by  monitoring  the  voltage  across  a 50-ohm  resistor 
connected  in  series  with  the  cell,  through  the  Tektronix  7D20 
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digitizer.  The  Fourier  transform  was  performed  on  the  digital  output 
of  the  digitizer.  As  a light  source  the  GR  1531 -AB  stroboscope  was 
used.  Long-pass  (A  > 1.1  /xm)  and  low-pass  (A  < 0.6  /an)  optical  filters 
were  used  to  create  the  uniformly  absorbed  and  strongly  absorbed  light 
pulses.  The  light  pulses  were  also  measured  with  a p-i-n  Si  detector 
for  the  strongly  absorbed  light  and  with  a Ge  detector  for  the 
uniformly  absorbed  light.  To  normalize  the  response,  the  Fourier 
transform  of  the  short-circuit  current  was  divided  by  the  Fourier 
transform  of  the  light  pulses.  The  amplitude  and  the  phase  of  the 
normalized  frequency  domain  response  are  shown  in  Figs.  4.1  and  4.2 
for  the  uniformly  absorbed  light,  and  in  Figs.  4.3  and  4.4  for  the 
strongly  absorbed  light.  The  quantum  efficiency  for  the  strongly 
absorbed  light  shining  from  the  back  contact,  as  measured  using  a He- 
Ne  laser,  was  about  0.7  %.  Using  the  procedure  described  in  the 
previous  section  the  best  fit  in  Figs.  4.1  and  4.2  was  found  for  rn  = 
12  ± 1 jxs , the  best  fit  in  Figs.  4.3  and  4.4  was  found  for  Dn  =>  17  ± 1 
cm^/s , and  from  the  quantum  efficiency  S-3 . 5x10^  ± 5x10^  cm/s.  The 
exact  value  of  S is  not  important  for  determining  rn  and  Dn  since  a 
change  in  S from  3x10^  to  4x10^  cm/s  has  a negligible  effect  on  the 
ratios  F(w)/F(0)  and  F' (w)/F' (0) . The  effect  of  the  lifetime  varying 
from  11  to  12  /xs  is  also  negligible  in  the  determination  of  Dn. 

We  believe  that  the  value  of  the  minority-carrier  diffusion 
coefficient  derived  above  is  a reliable  one.  This  is  supported  by  the 
good  agreement  between  the  experimental  results  and  theoretical 
expressions  for  the  strongly  absorbed  light  and  by  the  fact  that  the 
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UNIFORMLY  ABSORBED  LIGHT 


Fig.  4.1  Normalized  amplitude  of  the  frequency- domain  response  for 
uniformly  absorbed  light.  The  solid  line  is  the 
amplitude  of  the  measured  current  in  the  frequency  domain. 
The  broken  line  is  the  amplitude  of  the  ratio  F(«)/F(0) 
for  rn  = 12ns,  Dn  = 17  cnr/s  and  S - 3.5x10^  cm/s. 


PHASE  (DEGREES) 


89 


UNIFORMLY  ABSORBED  LIGHT 


Fig.  4.2  Phase  of  the  normalized  frequency- domain  response  for 

uniformly  absorbed  light.  The  solid  line  is  the  phase  of 
the  measured  current  in  the  frequency  domain.  The  broken 
line  is  the  phase  of  the  ratio  F(u)/F(0)  for  r = 12/xs, 
Dn  = 17  cm2/s  and  S = 3 . 5x10^  cm/s . 


NORMALIZED  AMPLITUDE 
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SQUARE  ROOT  OF  FREQUENCY  (Hz)**  1/2 


STRONGLY  ABSORBED  LIGHT 


Fig.  4.3  Normalized  amplitude  of  the  frequency- domain  response  for 
strongly  absorbed  light  shining  from  the  back.  The  solid 
line  is  the  amplitude  of  the  measured  current  in  the 
frequency  domain.  The  broken  lines  are  the  amplitude  of  the 
ratio  F'(w)/F'(0)  for  Dn  =17  ( ), 

Dn  “ 18  ( ) and  Dn  = 19  cm^/s  ( — — — ) . 

Here  rn  was  set  at  12  /*s  and  S at  3.5x10^  cm/s. 
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SQUARE  ROOT  OF  FREQUENCY  (Hz)**  1/2 
STRONGLY  ABSORBED  LIGHT 


Fig.  4.4  Phase  of  the  normalized  frequency- domain  response  for 

strongly  absorbed  light  shining  from  the  back.  The  solid 
line  is  the  phase  of  the  measured  current  in  the 
frequency  domain.  The  broken  lines  are  the  phase  of  the 

ratio  F'(w)/F'(0)  for  Dn  - 16  ( ), 

= 17  ( ) and  Dn  ■ 18  cm^/s  ( — — — ) . 

Here  rn  was  set  at  12  (is  and  S at  3.5xl04  cm/s. 
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phase  and  the  amplitude  yield  values  for  Dn  less  than  10  % apart.  Also 
the  corner  frequency  in  Fig.  4.3  corresponds  to  a lifetime  rn  = 

2/u'c  which  is  about  the  same  as  the  lifetime  derived  from  Fig.  4.1. 

The  value  Dn  - 17  ± 1 cm^/s  or  /*n  = 666  ± 39  cm^/Vs  at  23  °C  for 
minority  electrons  and  a doping  density  NA  = 1.22  x 10^  cm'^  , is 
within  experimental  error  equal  to  the  majority  electron  mobility  of 
696  cur- /Vs  for  the  same  donor  doping  density  as  reported  by  Thurber  et 
al  [39].  This  verifies  the  results  reported  by  Dziewior  and  Silber 
[40]  and  by  Neugroschel  [41] . 

The  frequency -domain  transient  analysis  method  can  also  be  used 
to  characterize  bulk  wafers  by  depositing  a Schottky  contact  on  the 
front  surface  and  a gridded  ohmic  contact  on  the  back  surface  and  then 
shining  both  the  strongly  and  the  uniformly  absorbed  pulsed  light  from 
the  back  surface.  It  can  also  be  used  in  dark  devices  by  using 
electrical  pulses,  as  in  the  electrical  short-circuit  current  decay. 

In  Si  devices,  however,  the  turn-on  part  of  the  response  is  determined 
by  space  charge  region,  emitter,  and  series  resistance  parasitics. 
Therefore,  only  the  turn-off  of  the  response  can  be  considered,  which 
corresponds  to  exciting  the  device  by  a step  input.  The  amplitude  of 
the  Fourier  transform  of  the  step  function  drops  as  1/w  and, 
consequently,  the  input  might  not  have  enough  power  to  excite  higher 
modes  in  the  device,  thus,  reducing  the  sensitivity  of  the  method  at 
high  frequencies.  This  also  points  out  the  advantage  of  using  light 
pulses  for  excitation  since  there  are  no,  at  least  for  some  frequency 
range , space-charge  region  and  resistance  limitations,  and  narrow 
pulsed  light  sources  are  readily  available. 
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The  lifetime  and  the  diffusion  coefficient  could  also  be 
extracted  from  the  Haynes -Shockley  experiment  by  fitting  the  time- 
domain  response  with  theoretical  expressions.  The  effect  of  r and  D, 
however,  spreads  over  the  entire  time-domain  pulse  response  and  the 
fit  would  be  a two  parameter  fit.  In  the  method  demonstrated  here  r 
and  D affect  different  parts  of  the  curve  and  the  two  parameters  can 
be  separated.  An  additional  drawback  of  the  time-domain  fitting  is 
mathematical  complexity:  the  pulse  response  of  the  system  would 
involve  an  infinite  series  of  analytical  functions  which  would  have  to 
be  convoluted  with  the  input  pulse. 

The  Fourier  transform  method  introduced  here  can  be  applied  to 
several  other  transient  techniques  such  as  photoluminescence  decay  or 
light  induced  open-circuit  voltage  decay.  In  the  latter,  care  should 
be  taken  that  the  terminal  voltage  is  low  enough  to  assure  linearity 
between  input  and  output,  and  that  the  capacitance  and  recombination 
in  the  space  charge  region  is  taken  into  consideration. 

So  far  in  Chapters  two,  three  and  four  we  implicitly  assumed  that 
when  a p/n  junction  is  short-circuited  and  illuminated  the  minority- 
carrier  density  at  the  space-charge  region  edge  of  the  base  is  zero. 
The  validity  of  this  assumption  is  discussed  in  Appendix  C. 

4 . 4 Conclusions 

In  conclusion  we  have  proposed  and  demonstrated  a new  method, 
frequency- domain  transient  analysis  (FDTA) , for  measuring  lifetimes 
and  diffusion  coefficients.  The  method  utilizes  the  whole  time-domain 
transient  response  and  through  Fourier  transform  gives  the  frequency 
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response  of  the  device.  We  showed  that  FDTA  is  more  plausible  and 
accurate  than  other  time-domain  techniques.  The  diffusion  coefficient 
and  lifetime  in  a p-type  0.16  ohm-cm  base  were  measured  and  found  to 
be  17  ± 1 cm^/s  and  12  jus , respectively.  The  method  can  be  applied  to 
characterize  solar  cells,  photodetectors  or  bulk  wafers.  Extension  to 
devices  operating  in  the  dark  is  possible  although  problems  arising 
from  emitter  or  space-charge-region  parasitics  can  reduce  the 
sensitivity  of  the  method  at  high  frequencies. 


CHAPTER  5 

NUMERICAL  SOLUTIONS,  PHYSICAL  INTERPRETATIONS,  AND  ANALYTICAL 
MODELS  FOR  AMORPHOUS  SILICON  P/I/N  SOLAR  CELLS 

5 . 1 Introduction 

This  chapter  presents  contact- to-contact  computer  solutions  of 
the  a-Si:H  p/i/n  solar  cell  and  treats  these  as  computer  experiments 
designed  to  reveal  the  approximations  and  insight  needed  for  the 
development  of  analytical  models.  A FORTRAN  77  program  was  written 
based  on  the  Poisson  equation  and  the  hole  and  electron  continuity  and 
current  equations.  We  ran  this  program  on  a Harris  800  computer, 
studying  many  aspects  of  internal  variables  such  as  electric  field, 
carrier  densities,  current  densities,  and  recombination  rates  versus 
position,  terminal  voltage  and  photon  flux  density.  Strongly  absorbed, 
weakly  absorbed  and  AMI . 5 light  were  studied  to  gain  the  insight 
needed  for  the  analytical  modeling.  Based  on  our  physical 
interpretation  of  the  computer  experiments,  we  propose  analytical  and 
equivalent-circuit  models  which  support  each  other  and  which  help 
explain  the  physical  origin  of  interdependencies  among  such  variables 
as  quantum  efficiency,  electric -field  and  recombination- rate 
nonuniformities  and  their  relations  to  the  current-voltage 
characteristic.  We  mathematically  treat  the  concept  of  the  limiting 
carrier  by  separating  the  current  into  photocollected  and 
back- injection  components  and  conclude  that  the  limiting  carrier  is 
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the  carrier  with  the  smaller  photocollected  current.  This  concept  is 
also  treated  in  our  critical  review  of  the  literature  where  we  employ 
reasoning  in  the  time  domain. 

This  Chapter  consists  of  five  sections.  Section  2 presents  the 
parameters  and  assumptions  used  for  the  numerical  simulations.  In 
section  3 the  numerical  solutions  are  interpreted,  the  concept  of  the 
limiting  carrier  is  emphasized  and  a new  equivalent -circuit  model  is 
developed.  Section  4 proposes  analytical  models  that  express 
current , voltages  and  recombination  rates  as  functions  of  material 
properties  such  as  lifetime,  mobility  and  absorption  coefficients.  In 
the  last  section,  we  review  selected  papers  concerning  the  modeling  of 
amorphous  p/i/n  devices. 

All  numerical  and  analytical  modeling  is  done  for  the  steady 
state . 


5 . 2 Numerical  Modeling 

Following  Hack  and  Shur  [42]  we  assume  an  exponential 
distribution  for  acceptor,  gA(E) , and  donor,  gD(E) , type  states  in  the 
mobility  gap,  as  shown  in  Fig.  5.1: 

gA(E)  =■  gAexp  [ ( E - Emc ) /Ea ] (5.1(a)) 

gD(£)  “ §DexP[-(E-Emc)/ED] • (5.1(b)) 

In  (5.1),  Ea  and  Ep  are  the  characteristic  energies  of  the 
acceptor  and  donor  type  traps.  In  amorphous  silicon  the  dominant 
recombination  mechanism  is  Shockley-Read-Hall  (SRH)  recombination 
through  the  gap  states  [43] . Assuming  a and  zero  temperature 
distribution  function  for  the  trapped  carriers  (this  assumption  is 
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Fig-  5.1 


Defect  states  spectrum  of  an  a-Si:H  solar  cell. 
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justified  since  EA  and  Eq  are  several  times  the  thermal  voltage  kT)  we 
derive  the  following  closed  form  expression  for  the  recombination 
rate : 

R - (Pn-ni2)  C Uth  an  {gA  EA  [Nc-1exp((Ec-Emc)/EA)((n+Cp)/Nc)kT/EA-1- 

exp((Ev-Emc)/EA)  CkT/EA+  gDED  [Nv'1exp( (EmcEv)/ED) ( (n+Cp)/Nv)kT/E^~ 1 

exp((-Ec+Emc)/ED)  C kT/ED  Nc\ (p+nC)/Nc) } . (5.2) 

Here,  Ut^  is  the  thermal  velocity  of  an  electron,  an  is  the  capture 
cross  section  of  a neutral  trap  and  C is  the  ratio  crc/an  where,  ac  is 
the  cross  section  of  a charged  trap.  Further,  n and  p are  the  carrier 
densities  and  Nc  and  Nv  the  effective  density  of  states  for  electrons 
and  holes  respectively.  The  intrinsic  carrier  density  is  given  by  n^  = 
[Nc  Nv  exp( -Eg/kT)  ] -*-/2  where  Eg  is  the  band  gap.  Under  the  same 
assumptions  regarding  density  of  states  and  occupation  statistics,  the 
expression  for  the  trapped  electron  density  becomes 

nt  = §A  EA  [exp((Ec-Emc)/EA)(n/Nc)((n+Cp)/Nc)kT/EA-1  + 

+ exp(-(Emc-Ev)/EA)(Cp/(n+Cp))((p+n/C)/Nv)'kT/EA  ],  (5.3(a)) 

while  for  the  trapped  holes  becomes 

Pt  - §D  eD  [(p/Nv)exp((Emc-Ev)/ED)((n+Cp)/Nvr/ED'1+ 

exp((-Ec+Emc)/ED)(nC/(nC+p))((n+p/C)/Nc)'kT/ED] } . (5.3(b)) 

The  expression  for  the  total  charge  density  is 


Q(n,p)  - e(pt  * nt  + ND  - NA) , 


(5.4) 
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where  e is  the  elementary  charge,  and  Np  - NA  is  the  net  ionized 
doping  density. 

The  cell,  Fig.  5.2,  is  assumed  to  have  an  overall  thickness  of 
0.6  /im.  The  doped  layers  have  a thickness  of  100  A and  the  light 
enters  through  the  P+  layer.  The  boundary  condition  at  the  ohmic 
contact  on  the  P+  layer  is  that  the  majority  carrier  (hole)  density 
remains  constant  at  p(0)  = 4.6xl0-*-^  cm"^  while  the  minority  carriers 
(electrons)  see  the  contact  as  a surface  characterized  by  a 
recombination  velocity  equal  to  Sn  - 5x10^  cm/s  (kinetic  limit 
velocity  for  electron  in  single-crystal  silicon).  Similarly  the 
boundary  condition  at  the  N+  layer  is  that  the  electron  density  stays 

i r o 

constant  at  N(W)  — 10-*-°  cm'~>  while  holes  see  a recombination  velocity 
of  Sp  - 2x10^  cm/s . To  correlate  the  doping  density,  the  trap  density 
and  the  majority  carrier  density  at  the  ohmic  contacts,  we  have 
assumed  that  in  equilibrium  the  net  charge  density  at  the  ohmic 
contact  is  zero.  Taking  the  density  of  ionized  dopants  in  the  heavily 
doped  layers  to  be  10-*-^  cm'-*-  forces  gp  in  the  P+  layer  to  be  g+  = 
2.5x10-^  cm" 3 eV*l.  In  the  N+  layer  gA  = g+  = 4. 2x10-*- ^ cm'^  eV'^.  In 
the  intrinsic  layer  the  values  of  gA  and  gp  were 
SA  ” §D  ” 0.2x10-*-^  cm eV* ^ . 

Two  cases  of  optical  monochromatic  excitation  are  considered.  In 
the  first,  the  absorption  coefficient  a - a(A)  is  high  (a  =2xl05  cm'-*-) 
and  the  light  is  strongly  absorbed  near  the  front  surface.  In  the 
second  case  the  absorption  coefficient  is  small  (a  - 2x10^  cm'-*-),  the 
light  is  reflected  once  on  the  back  surface.  The  reflection 
coefficient  is  0.8,  and  the  generation  rate  is  almost  uniform 
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Fig.  5.2  Band  diagram  of  an  a-Si:H  solar  cell. 
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throughout  the  cell.  In  addition,  the  effect  of  AMI . 5 light 
shining  through  the  P+  layer  is  studied  and  related  to  the 
monochromatic  excitations.  The  equations  used  in  the  numerical 
simulation  were: 


Jn  = ^n  n (dFn/dx) 

(5.5) 

Jp  “ bp  P (dFp/dx) 

(5.6) 

dJn/dx  = - e(G  - R) 

(5.7) 

dJp/dx  - e(G  - R) 

(5.8) 

dE/dx  - Q/e 

(5.9) 

d(E^/e)/dx  = E. 

(5.10) 

Equations  (5.5)  and  (5.6)  express  the  electron  and  the  hole  current  in 
terms  of  the  gradient  of  the  respective  quasi-Fermi  levels.  Equations 
(5.7)  and  (5.8)  are  the  continuity  equations  for  electrons  and  holes, 
(5.9)  is  the  Poisson's  equation  and  (5.10)  expresses  the  electric 
field  as  the  gradient  of  the  electrostatic  potential  on  the  electron 
energy  scale.  The  standard  relations,  n = n^  exp((Fn  - E^)/kT)  and  p 
= n^  exp((E^  - Fp)/kT) , connect  Fn  and  Fp  with  the  rest  of  the 
equations . 

The  system  of  equations  (5. 5) -(5. 10)  was  solved  using  Newton's 
method.  The  results  appear  in  the  next  section.  We  conclude  this 
section  by  presenting  a list  of  the  parameters  used  in  the  numerical 


simulation. 
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Table  1. 

List  of  Parameters  Used  in  Numerical  Modeling 


= 10  cm2/(Vs) 
A*p  - 1 cm2/ (Vs) 
gA  = 0.2x10^  cm‘2eV"^ 


Electron  mobility 
Hole  mobility 

Acceptor  states  preexponential  factor 
Donor  states  preexponential  factor 
Acceptor  states  characteristic  energy 
Donor  states  characteristic  energies 
Energy  of  gap  states  minimum 
Mobility  gap 

Effective  density  of  states 
Capture  cross  section  x thermal  velocity  (neutral  trap)  = 10'H  cm2s‘^ 


= 0.2x10 

16  cm"2 

eV-1 

ea 

= 0.053 

eV 

ed 

= 0.088 

eV 

Emc 

- 0.65 

eV 

Eg 

= 1.72 

eV 

Nc  - Nv  ■ 

= 1019  i 

cm" 2 

Capture  cross  section  x thermal  velocity  (charged  trap)  = 10" 9 cm2s"^ 


5 . 3 Interpretation  of  the  Numerical  Results 
Figure  5.3(a,b)  presents  the  numerical  results  on  the  electric 
field  magnitude  profiles  at  short-circuit  conditions  as  a function  of 
the  photon  flux  density  of  strongly  and  weakly  absorbed  light.  We  can 
see  that  the  field  is  much  more  disturbed  from  its  equilibrium  value 
for  uniform  absorption  than  it  is  for  strong  absorption.  For  weakly 
absorbed  light  the  field  has  a deep  minimum  near  the  N+  layer,  while 
for  strongly  absorbed  light  the  field  is  relatively  uniform  with  its 
minimum  shifted  towards  the  P+  contact.  Figure  5.4(a,b)  shows  the 
effect  an  applied  forward  bias  has  on  the  field  profile.  For  strong 
absorption  the  applied  voltage  reduces  the  field  mainly  in  the 
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(b) 


Fig.  5.3  Electric  field  profiles  at  short  circuit  as  functions  of 
the  photon  flux  density,  (a)  Strongly  absorbed  light, 
curves  1,2,3  correspond  to  F = 101',  1016-  0 cm^s’1, 
respectively,  (b)  Uniformly  absorbed  light,  curves  1,  2,  3 
correspond  to  F — 10^®,  10  0 cm'^s"^,  respectively. 
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(a) 


(b) 


Fig.  5.4  Electric  field  profiles  as  functions  of  forward  bias  at 
short-circuit,  1,  maximum-power,  2,  and  open-circuit,  3. 

(a)  Strongly  absorbed  light,  F = lO17  cm'^s'1,  V » 0, 

0.725,  0.927  Volts  (1  2,  3,  respectively.),  (b)  Uniformly 

absorbed  light,  F = 1018  cm‘2s'1,  V = 0,  0.760,  0.918  Volts 
( 1,2,3  respectively). 
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intrinsic  layer;  large  voltages  reverse  the  direction  of  the  field  in 
the  bulk  of  the  i-layer.  For  uniform  absorption  the  field  drops  mainly 
at  the  p/i  interface  and  its  minimum  value  is  only  slightly  affected. 
In  both  cases  the  position  of  the  minimum  is  relatively  bias 
independent . 

To  understand  the  effect  the  generation  profile  has  on  the 
electric  field,  we  have  to  take  into  account  the  difference  between 
the  electron  and  the  hole  collection  length  (lc  = E^r) . In  this  model 
electrons  have  higher  mobility  (nn  = 10  cm^/Vs)  than  holes  (^p  = 1 
cm^/Vs).  In  addition,  the  electron  lifetime  is  higher  than  the  hole 
lifetime  because  there  are  more  hole  traps  than  electron  traps  (EA  < 
Eq) . Because  of  the  longer  electron  collection  length  and  under 
uniform  excitation,  more  photogenerated  electrons  arrive  at  the 
electron  ohmic  contact  than  holes  at  the  hole  ohmic  contact.  Since  in 
one  dimension  the  total  current  must  be  position  independent,  the 
electric  field  is  distorted  to  forward  bias  the  i/n  interface  and  to 
cause  the  subsequent  electron  back- inj ection  current  compensating  for 
the  difference  between  the  two  photocurrents. 

This  is  shown  in  the  equivalent -circuit  model  of  Fig.  5.5.  The 
current  sources  1^  and  I2  are  the  electron  and  hole  photocollected 
current  (defined  in  the  next  section) . Diodes  and  D2  represent  the 
p/i  and  i/n  interfaces,  respectively,  while  D3  accounts  for  the 
carrier  injection  caused  by  the  terminal  voltage.  The  limiting  carrier 
is  the  carrier  with  the  least  photocollected  current.  For  uniform 
excitation  l2»I]_,  holes  are  the  limiting  carriers  and  is  reverse 
biased  while  D2  is  forward  biased  under  short-circuit  conditions.  This 
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Fig.  5.5 


Equivalent-circuit  model  for  an  a-Si:H  solar  cell. 
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explains  the  high  field  density  at  the  p/i  interface  and  the  deep 
minimum  near  the  i/n  interface  shown  in  Fig.  5.3(b).  The  equivalent 
circuit  of  Fig.  5.5  also  explains  why  the  strongly  absorbed  light  only 
slightly  perturbs  the  electric  field, as  shown  in  Fig.  5.3(a):  almost 
all  carriers  are  generated  near  the  P+  layer  and  this  compensates  for 
the  shorter  hole  collection  length  and  makes  electrons  the  limiting 
carriers.  Thus  1^  is  only  slightly  larger  than  I2  and  the  bias  across 
and  D2  is  low. 

When  a terminal  voltage  is  applied,  under  uniform  excitation  this 
voltage  drops  almost  entirely  across  since  a very  small  voltage 
drop  across  the  forward  biased  D2  can  accommodate  any  difference  in 
the  terminal  current.  This  accounts  for  the  sharp  drop  of  the  field 
with  increased  bias  at  the  p/i  interface  and  the  minor  dependence  of 
the  minimum  value  of  the  field  on  terminal  bias  ( Fig.  5.4(b)).  In  the 
case  of  strong  absorption  the  applied  voltage  is  nearly  evenly  divided 
between  and  D2.  Field  changes  in  the  i-layer  sustain  the  applied 
voltage  and  the  field  minimum  becomes  sensitive  to  the  external  bias 
(Fig.  5.4(a)). 

The  previous  analysis  about  the  dependence  of  the  voltage  across 
D]_  and  D2  on  excitation  can  explain  the  recombination  profiles  shown 
in  Fig.  5.6(a,b).  For  uniform  excitation  the  large  difference  between 
currents  1^  and  I2  circulates  through  D2  and  creates  the  high 
recombination  rate  near  the  i/n  interface  (Fig.  5.6(b)).  At  this 
interface  the  recombination  rate  is  almost  terminal  bias  independent 
because  the  voltage  across  D2  is  nearly  terminal  bias  independent.  For 
strong  absorption  the  recombination  rate  at  short  circuit  is  higher 
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(a) 


(b) 


Fig.  5.6  Recombination  rate  profiles  as  functions  of  the  forward 

terminal  bias,  (a)  Strongly  absorbed  light,  (b)  Uniformly 
absorbed  light.  Curves  1,  2,  3 correspond  to  short-circuit, 
maximum-power  and  open-circuit  conditions,  respectively. 
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near  the  p/i  interface  than  it  is  in  the  rest  of  the  device  because  of 
the  slight  forward  bias  across  Dj_.  For  higher  terminal  voltages  the  pn 
product  depends  more  on  the  terminal  bias  and  less  on  the  generation 
profile  resulting  in  flatter  recombination  profiles  as  shown  in  Fig. 
5.6(a). 

Now  we  shift  our  attention  to  AMI. 5 light.  The  recombination  rate 
(Fig.  5.7)  is  higher  near  the  N*  layer  (especially  at  low  applied 
voltages)  since  electrons  have  higher  photocollected  current  and  the 
recombination  current  through  D2  is  higher  than  through  D^_  (Fig.  5.5). 
Figure  5.7  also  shows  the  generation  rate  profile  which  is  shown  to  be 
higher  than  the  recombination  rate  in  the  front  part  of  the  device 
even  at  open  circuit.  The  difference  between  bulk  generation  and 
recombination  is  compensated  by  the  electron  surface  recombination  at 
the  P+  ohmic  contact.  This  surface  recombination,  caused  primarily  by 
the  low  hole  density  in  the  P+  layer  and  the  subsequent  high  electron 
injection  in  this  layer,  is  the  major  loss  mechanism  for  the  device 
under  study.  Another  aspect  of  Fig.  5.7  is  that  AMI. 5 light  results  in 
a recombination  profile  having  the  same  qualitative  features  as  the 
recombination  profile  of  the  strongly  absorbed  light. 

In  Fig.  5.8  the  electric  field  profile  is  shown  at  short-circuit 
maximum-power  and  open-circuit  conditions.  We  see  that  it  combines 
features  from  both  the  strongly  absorbed  and  the  uniformly  absorbed 
light  field  profiles  (Fig.  5.4(a,b)).  This  indicates  that  the  concept 
of  the  limiting  carrier  is  not  so  well  defined  under  AMI. 5 light  as  it 
is  for  strongly  and/or  uniformly  absorbed  light.  Fig.  5.9  shows  the 
electron  and  hole  current  profiles.  The  sharp  increase  of  the  electron 
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Fig-  5.7  Recombination  profiles  as  functions  of  the  applied  forward 
bias  for  AMI. 5 light.  Curves  1,  2,  3 correspond  to  V = 0 
(short-circuit),  0.740  (maximum -power)  and  0.922  Volts 
(open-circuit),  respectively.  Curve  4 is  the  generation 
profile. 
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Fig.  5.8  Electric  field  profiles  as  functions  of  the  applied  forward 
bias  for  AMI. 5 light.  Curves  1,  2,  3 correspond  to  V = 0 
(short-circuit),  0.740  (maximum -power)  and  0.922  Volts 
(open-circuit),  respectively. 
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Fig.  5.9  Electron  and  hole  current  profiles  as  functions  of  the 
applied  forward  bias  for  AMI . 5 light.  Curves  3,2,1 
(hole  current)  and  4,  5,  6 (electron  current)  correspond  to 
V - 0 (short-circuit),  0.740  (maximum -power)  and  0.922 
Volts  (open-circuit) , respectively. 
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current  near  the  P+  contact  is  due  to  the  high  generation  rate  at  this 
contact.  For  voltages  near  the  open-circuit  voltage  the  electron 
current  drastically  decreases  with  bias  and  reverses  sign  at  the  front 
part  of  the  cell  because  of  the  back- inj ection  into  the  P+  layer. 
Figures  5.10  (a,b)  present  the  carrier  density  profiles.  The  carrier 
density  near  the  p/i  interface  is  more  sensitive  to  forward  voltage 
because  the  applied  voltage  mainly  drops  in  the  reverse  biased  D]_ 
diode  of  the  equivalent -circuit  model  of  Fig.  5.5.  From  Fig.  5.8,  5.9 
and  5.10  one  concludes  that  the  order  of  magnitude  of  the  electron 
density  in  the  i-layer  can  be  expressed  by  the  drift  only 
approximation  as  n = Jn/(q  Mn  E)  . A similar  relation  holds  for  holes. 
The  approximation  is  better  for  low  forward  voltages  where  the  field 
is  still  high  in  the  i-layer. 

The  carrier  lifetimes , defined  as  the  ratio  of  the  carrier 
densities  over  the  recombination  rate,  are  shown  in  Fig.  5.11  (a,b) . 
Near  the  N+  layer  the  carrier  lifetimes  are  basically  bias  independent 
because  the  applied  voltage  drops  almost  entirely  on  D^.  As  a result 
the  hole  lifetime  decreases  with  forward  voltage  near  the  p/i 
interface.  At  the  same  interface  the  electron  lifetime  increases 
because  the  electron  density  increases  faster  than  the  recombination 
rate  with  forward  bias.  In  the  heavily  doped  layers  the  carrier 
lifetime  is  bias  independent  because,  there,  the  minority-carrier 
density  and  recombination  rate  are  linearly  proportional  to  each 
other.  This  linear  relation  can  be  derived  from  (5.2)  when  p » n or 


n » p . 
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(a) 


(b) 


Fig.  5.10  Free  carrier  density  profiles  as  functions  of  the  applied 
forward  bias  for  AMI . 5 light,  (a)  Electron  density,  (b) 
Hole  density.  Curves  1,  2,  3 correspond  to  V = 0 (short- 
circuit),  0.740  (maximum- power)  and  0.922  Volts  (open- 
circuit),  respectively. 
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(a) 


(b) 


Fig.  5.11  Free  carrier  lifetime  profiles  as  functions  of  the  applied 
forward  bias  for  AMI . 5 light,  (a)  Electron  lifetime,  (b) 
Hole  lifetime.  Curves  1,  2,  3 correspond  to  V = 0 (short- 
circuit),  0.740  (maximum -power)  and  0.922  Volts  (open- 
circuit)  , respectively. 
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5.4  Analytical  Modeling 


5.4.1  Definition  of  Photocollected  Currents 

To  define  mathematically  what  1^  and  I2  in  Fig.  5.5  are,  we  use 
our  numerical  results  for  the  electric  field,  E(x)  and  the  carrier 
recombination  lifetimes,  rn  and  r p,  to  separate  the  currents  into 
photocollected  and  back- inj ection  components.  For  a given  generation 
profile  the  electron  (hole)  photocollected  current  results  from  the 
photogenerated  electrons  (holes)  that  cross  the  cell,  under  the 
influence  of  the  numerically  calculated  electric  field  E(x)  and 
lifetime  rn  (r p) , and  become  collected  by  the  electron  (hole)  contact. 
On  the  other  hand,  the  back- inj ection  electron  (hole)  current  is  the 
current  that  would  flow  in  the  dark  if  the  electric  field  were  E(x) 
and  the  lifetime  rn  (rp) • In  calculating  the  photocollected  and  back- 
injection  currents  the  recombination  velocities  at  the  ohmic  contacts 
retain  the  values  assigned  in  the  introduction,  Sn  and  Sp . 

In  more  analytical  terms  the  electron  photocollected  current 
is  the  solution  of  the  system  of  the  equations : 


The  boundary  condition  for  (5.11)  at  the  P+  ohmic  contact  is 
represented  by  the  recombination  velocity  for  electrons  Sn,  JA(0)  = 
eSnn' (0)  , while  at  the  N+  ohmic  contact  is  n'(W)=0.  On  the  other  hand, 
the  equations 


JA  - e nn  n'  E + e Dn  dn'/clx 


(5.11(a)) 


dJ^/dx  ~ - e(G-n'/rn) • 


(5.11(b)) 


Jbin  = e Mn  n"  E + e Dn  dn"/dx 


(5.12(a)) 
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^bin/ dx  “ e n"/rn>  (5.12(b)) 

define  the  electron  back  injection  component,  Jb:Ln.  The  boundary 
conditions  for  (5.12)  at  the  P+  ohmic  contact  are  represented  by  the 
same  recombination  velocity  for  electrons  Sn,  while  at  the  IT1"  ohmic 
n"(W)  — N(W)  - 1016  cm'^,  where  N(W)  is  the  equilibrium  electron 
density  at  the  electron  ohmic  contact.  In  (5.11)  and  (5.12),  E and  rn 
are  the  numerically  calculated  electric  field  and  electron  lifetime. 
Equations  similar  to  (5.11)  and  (5.12)  hold  for  holes  also. 

Based  on  the  apparent  linearity  of  the  carrier  transport 
equations,  we  can  see  that  under  optical  and  terminal  excitation  the 
electron  (hole)  current  is  the  sum  of  the  photocollected  and  back- 
injection  components: 


J = Jph  + Jbi-  (5-13) 

We  emphasize  that  this  separation  of  currents  does  not  imply  that  the 
superposition  principle  holds  for  a-Si:H  p/i/n  solar  cells.  The 
dependence  of  the  electric  field  and  carrier  lifetimes  on  bias 
invalidates  this  principle.  The  separation  of  currents  attempted  here 
only  seeks  to  clarify  the  concept  of  the  limiting  carrier. 

According  to  the  model  shown  in  Fig.  5.5  the  carrier  with  the 
least  photocollected  current  is  the  limiting  carrier  and  determines 
the  cell  photocurrent.  The  photocurrent  is  defined  as  the  algebraic 
difference  between  the  terminal  current  under  illumination  and  the 
terminal  current  in  the  dark  at  the  same  voltage.  To  support  the 
predictions  of  this  model  we  present  in  Fig.  5.12(a,b)  numerical 
results  on  the  cell  terminal,  photocollected  and  back- inj ection 
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Fig.  5.12  Current-Voltage  characteristics,  (a)  Strongly  absorbed 

light,  (b)  Uniformly  absorbed  light.  Curves  1,  2,  3,  4,  5, 
6 are  the  terminal  current,  electron  photocollected 
current,  hole  photocollected  current,  electron  back- 
injection  current,  hole  back- inj ection  current  and  cell 
photocurrent,  respectively. 


119 


current  densities  for  strongly  and  weakly  absorbed  light.  As  discussed 
in  the  previous  section,  for  strongly  absorbed  light  the  carriers 
(electrons)  associated  with  the  nonilluminated  contact  are  the 
limiting  carriers.  This  is  illustrated  in  Fig.  5.12(a),  where  the 
electron  photocollected  current  is  shown  to  be  very  close  to  the 
terminal  current  and  even  closer  to  the  terminal  photocurrent.  The 
higher  hole  photocollected  current  forward  biases  the  p/i  interface 
and  creates  the  hole  back- inj ection  current.  The  electron  back- 
injection  current  is  negligible  at  low  terminal  voltage,  where  the  i/n 
interface  is  reverse  biased,  and  increases  at  voltages  near  Voc,  where 
both  interfaces  are  forward  biased.  The  hole  back- inj ection  current 
and  the  back  diffusion  of  photogenerated  electrons  compensate  for  the 
difference  between  the  two  photocurrents.  This  difference  is  small 
resulting  in  high  internal  quantum  efficiency  85%)  for  the  strongly 
absorbed  light.  For  uniformly  absorbed  light  (Fig.  5.12(b))  the  large 
difference  between  the  electron  and  the  hole  photocollected  currents 
reduces  drastically  the  cell  photocurrent  and  the  internal  quantum 
efficiency  (~  45%).  At  low  voltages  the  small  difference  (==.5  mA) 
between  the  hole  photocollected  current  and  the  cell  photocurrent  is 
due  to  the  back  diffusion  of  photogenerated  electrons  in  the  P+  layer. 
At  voltages  near  Voc  this  difference  becomes  larger  due  to  the 
disturbance  of  the  electric  field  in  the  P+  layer.  The  high  electron 
back- inj ection  current  implies  the  strong  forward  bias  of  the  i/n 
interface . 

Figure  5.13  shows  the  same  currents  for  AMI . 5 light.  Holes  are 
here  the  limiting  carriers  because  their  photocollected  current  is 


(mA/cm 
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Fig.  5.13  Current- Voltage  characteristics  for  AMI . 5 light.  Curves  1, 
2,  3,  4,  5,  6 are  the  terminal  current,  electron 
photocollected  current,  hole  photocollected  current, 
electron  back- inj ection  current,  hole  back- inj ection 
current  and  cell  photocurrent,  respectively. 
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less  than  the  electron  photocollected  current.  The  difference  between 
the  two  currents,  however  is  small  since  the  high  generation  rate  near 
the  P+  layer  compensates  for  the  smaller  hole  collection  length.  For 
the  device  under  study  the  AMI. 5 short-circuit  current,  open-circuit 
voltage,  fill  factor  and  conversion  efficiency  were  13.93  mA,  0.922  V, 
0.614  and  8.2%,  respectively. 

5.4.2  Analytical  Solutions  for  Strongly  Absorbed  Light 

Based  on  the  above  discussion,  we  can  write  the  cell  photocurrent  as 

Jph  “ JL  + Jbd>  (5-14) 

where  JL  is  the  limiting  carrier  photocollected  current,  calculated  at 
the  respective  ohmic  contact,  and  is  a correction  term  to  include 
the  back  diffusion  of  photogenerated  minority  carriers  near  the  same 
contact.  To  be  able  to  derive  closed  form  expressions,  the  electric 
field  and  carrier  lifetimes  are  assumed  spatially  independent.  The 
limitations  of  this  assumption  are  discussed  later  in  this  and  the 
last  section.  Under  these  assumptions  and  for  strongly  absorbed  light 
shining  through  the  P+  layer  the  expression  for  the  cell  photocurrent 
becomes,  as  derived  in  Appendix  B, 

2a[l/Lg  + 1/L  exp(W/L2) 

Jph  - Jn  - e F — — ,(5.15) 

a/  + 2a/LE  - 1/1/  1 - exp[W(l/L2  - 1/LX) ] 

where  Lg=2 (kT/e | E | ) , L is  the  diffusion  length,  L=(Dr)l/2,  F is  the 
photon  flux  density  and  and  L2  are  defined  as 


L]_  - [1/Le  + (I/Le+I/L2)1/2]!1  L2  = [1/Le  - (1/Le+1/L2)1/2]:1  (5.16) 
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In  (5.15)  the  cell  photocurrent  was  set  equal  to  the  electron 
photocollected  current  since  the  hole  back- diffusion  current  near  the 
N'1'  contact  is  negligible.  If  the  electron  collection  length  lc 
satisfies  the  relation  lc  » 4kT/(e|E|),  which  is  normally  satisfied 
for  voltages  below  Voc,  then  (5.15)  can  be  simplified  to 


Jph 


e | E 

/kT 

a + e 

E | /kT 

exp(-W/lc) . 


(5.17) 


In  the  drift  only  approximation  exp(-W/lc)  is  the  probability 
that  an  electron  generated  at  x = 0 is  collected  by  the  electron  ohmic 
contact  located  at  x - W.  Therefore,  if  there  were  no  diffusion,  the 
photocurrent  would  be  Jp^  — e F exp(-W/lc)  since  almost  all  carriers 
are  generated  near  x = 0 because  of  the  strong  absorption.  The  term 
e E /(a  + e E /kT)  accounts  for  the  back  diffusion  of  electrons 
generated  near  the  hole  contact  [43].  To  make  (5.15)  and  (5.17)  better 
correspond  to  the  actual  field  distribution,  we  introduce  in  these 
equations  the  field  magnitude  at  the  p/i  interface,  Ep , and  the  field 
magnitude  in  the  i-layer  bulk,  E^.  Equation  (5.15)  now  becomes 


2«[1/Lep  + 1/L2]V2  exp[W(l/L2)] 

Jph  * e F . (5-18) 

a2  + 2a/LE  - 1/L2  exp[W(l/L2  - l/I^)  ] 


where  LEp  - 2kT/(eEp)  , LE^  = 2kT/(eE^)  , and  L^_ , L2  are  given  by  (5.16) 
where  now  LE^  replaces  LE.  Similarly,  (5.17)  becomes 


Jph 


e F 


e|Ep|/kT 


a + e 


[e^IAt 


exp(-W/lc) , 


(5.19) 


where  in  (5.19)  lc  = nnE^r . In  this  model  we  assume  that  only  Ep  and 
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E}.,  change  with  terminal  voltage.  The  value  of  the  field  at  the  p/i 
interface  is  of  the  order  of  10^  V/cm  and  exhibits  a weak  bias 
dependence  (Fig.  4(a)).  A suitable  expression  for  the  field  dependence 
on  the  terminal  voltage  V is 

Ep  = 0.8xl05  (Vp  - V)/Vp , (5.20) 

where  Vp  is  a parameter  to  be  determined  by  the  numerical  fitting.  The 
field  in  the  intrinsic  layer  is  expressed  as 

Eb  - <V0  - V)/W,  (5.21) 

where  VQ  is  a characteristic  voltage  defined  by  the  property  that 
V - VQ  makes  the  average  field  zero  in  the  intrinsic  layer.  This 
characteristic  voltage  normally  assumes  values  in  the  vicinity  of  Voc. 
The  best  fitting  was  obtained  for  Vp  - 1.2  V and  VQ  — 0.88  V.  Figure 
5.14  shows  the  numerical  and  analytical  photocurrent  calculations  for 
strongly  absorbed  light.  We  can  see  that  there  is  good  agreement 
between  numerical  and  analytical  results . The  parameter  to  which  the 
fitting  is  most  sensitive  is  the  limiting  carrier  fir  product.  The 
value  of  A‘nrn>  1-4-xlO'^  cm^/V,  used  in  Fig.  5.14  agrees  with  the 
numerical  results  for  /inrn.  Therefore,  by  fitting  experimental  I-V 
characteristics  for  strongly  absorbed  light  to  the  analytical 
expression  (5.19),  one  can  extract  the  limiting  carrier  fir  product. 

5.4.3  Analytical  Solutions  for  Uniformly  Absorbed  Light 

For  uniform  absorption  holes  are  the  limiting  carriers.  Under  the 
assumptions  of  uniform  field  and  lifetimes  the  hole  photocollected 


(A/cm2 ) 
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Fig.  5.14  Numerical  and  analytical  photocurrents  versus  forward 

bias,  for  strongly  absorbed  light.  Curves  (1)  and  (2)  are 
derived  from  (5.18)  and  (5.19),  respectively,  and  curve  3 
is  the  numerically  calculated  photocurrent. 
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current  density  can  be  expressed,  as  shown  in  Appendix  B,  as 


L2[(l/L2  - l/L1)exp(W/L2)  + l/I*  - exp(W/L2  - W/Lj.)/!*] 


Jp  = eG 


1 - exp(W/L2  - W/Lx) 


(5.22) 


where  L]_  and  L2  are  given  by  (5.16),  G is  the  generation  rate  per  unit 
volume,  and  L is  now  the  hole  diffusion  length.  If  lcu  is  the  limiting 
carrier  (holes)  collection  length,  lcu  = | E | /Xp r p , and  if  lcu  » 
4kT/e|E|,  a condition  satisfied  for  voltages  below  Voc,  then  (5.22) 
can  be  written  as 


Jp  - qG[lcu(l  - exp(-W/lcu))  - (kT/q|E|)(-l  + 2exp( -W/lcu) ) ] . (5.23) 


The  first  term  in  the  square  bracket  in  (5.23),  lcu(l  - exp ( -W/lcu) ) , 
dominates  and  is  the  same  as  the  expression  for  the  internal  quantum 
efficiency  in  the  Hecht  formula,  also  derived  by  Crandall  [44],  We 
emphasize  that  in  (5.23)  lcu  is  the  collection  length  of  the  carrier 
with  the  least  photocollected  current  and  not  the  sum  of  the  electron 
and  hole  collection  as  concluded  in  [44] . Since  the  Hecht  formula  is 
based  on  the  drift  approximation,  we  can  consider  the  second  term  in 
(5.23)  as  a correction  to  include  diffusion.  If  lc  is  the  electron 
collection  length,  then,  based  on  (5.23),  the  electron  back-diffusion 
current  near  the  P+  contact  can  be  written  as 


Jbd  = -qG(kT/q|E|)(-l  + 2exp ( -W/lc) ) . 


(5.24) 


From  (5.14)  the  cell  photocurrent  density  becomes 
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Jph  " eG[lcu(l  - exp ( -W/lcu)  - (kT/e | E | ) x 

( -2+2exp( -W/lcu)+2exp ( -W/lc) ) ] . (5.25) 

The  dominant  term  in  (5.25)  is  lcu(l  - exp(-W/lcu))  because  of  the 
high  value  of  |e|.  The  two  other  terms  represent  diffusion  losses.  As 
we  have  seen  in  the  numerical  results  section,  under  uniformly 
absorbed  light  a total  generation  rate  of  1.88x10^  cnT^  s ‘ ^ strongly 
perturbs  the  electric  field  the  minimum  of  which  near  the  i/n 
interface  is  slightly  affected  by  a small  forward  bias.  Such  a bias 
mainly  drops  at  the  p/i  interface.  Therefore,  to  account  for  the  field 
nonuniformities  in  distance  and  bias,  we  introduce  in  (5.25)  a 
quadratic  dependence  of  the  field  strength  on  the  bias: 

|E|  - [Vbi  - V(V/VN)]/W,  (5.26) 

where  is  the  built-in  barrier  height  and  VN  a normalization 

parameter  accounting  for  the  field  nonuniformities.  The  best  fit 
between  analytical  and  numerical  results  were  obtained  for  = Vb*  — 
1.16  V.  Figure  5.15  shows  numerical  and  analytical  calculations  for 
the  photocurrent  density  as  a function  of  the  voltage.  Curve  1 
corresponds  to  (5.25)  while  curve  2 is  the  plot  of  the  Hecht  formula 
Jph  =■  qGlcu[l  -exp(W/lcu) ] . These  analytical  formulas  are  in  good 
agreement  with  the  numerical  photocurrent  (curve  4)  for  voltages  up  to 
VQC . The  value  of  the  hole  pr  product  used  to  generate  curves  1 and  2 

Q O 

is  1.5x10" y ctar/V,  a value  close  to  the  numerical  results.  At  voltages 
near  VQC  the  strong  perturbation  of  the  field  in  the  P+  layer  causes  a 
strong  electron  current  towards  the  P+  contact.  This  is  not  predicted 
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Fig.  5.15  Numerical  and  analytical  photocurrents  versus  forward 
bias,  for  uniformly  absorbed  light.  Curves  (1)  and  (2) 
are  derived  from  (5.22)  and  the  Hecht  formula, 
respectively,  curve  3 is  the  hole  photocollected  current 
and  curve  4 is  the  numerically  calculated  photocurrent. 
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by  (5.25).  We  think  this  is  the  reason  for  the  disagreement  between 
numerical  and  analytical  results  near  Voc.  This  view  is  supported  by 
the  good  agreement  among  curves  1 , 2 and  curve  3 which  is  the 
numerically  calculated  hole  photocollected  current  density.  Note  that 
if  the  total  generation  rate  were  much  smaller  (=10^  cm'^s'^)  the 
field  would  be  more  uniform  and  we  could  model  the  field  dependence  on 
the  bias  by  simply  setting  |e|  = (V^  - V)/W.  In  Appendix  B an 
analytical  formula  is  presented  for  the  photocurrent  under  AMI. 5 
illumination. 

5 . 5 Discussion  and  Critical  Review  of  the  Literature 

A basic  concept  used  in  the  modeling  of  a-Si:H  is  the  concept  of 
the  limiting  carrier  [42],  [44],  [45],  [46].  Determination  of  the 
limiting  carrier  provides  insight  about  mechanisms  that  constrain  the 
quantum  efficiency  and  open- circuit  voltage  as  well  as  about  the 
asymmetric  performance  of  the  cell  when  illuminated  from  different 
sides.  Crandall  [44]  and  Okamoto  [45]  conclude  that  the  cell 
performance  is  determined  by  the  sum  of  the  hole  and  the  electron 
collection  lengths,  respectively.  Hence,  they  conclude,  the  carrier 
with  the  longest  collection  length  is  the  limiting  carrier.  Similar 
results  have  been  reached  by  Sichanugrist  [46] . In  the  above  models 
electric  field  and  carrier  lifetimes  were  assumed  spatially 
independent.  To  understand  why  the  above  conclusions  are  incorrect  we 
have  to  reconsider  the  legitimacy  of  the  constant  field  approximation. 

To  do  so  we  introduce  a phenomenological  way  of  thinking  and 
follow  the  sequence  of  events  that  take  place  after  a sudden 
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application  of  an  optical  excitation  to  a short-circuited  cell.  We 
assume  a uniform  optical  generation  and  that  electrons  have  higher 
mobility  and  lifetime  than  holes.  Immediately  after  the  excitation 
photogenerated  electrons  and  holes  will  move  towards  their  respective 
contacts.  Because  electrons  have  higher  mobility  and  lifetime  than 
holes , the  electron  photogenerated  current  collected  at  the  N"*-  layer 
is  higher  than  the  hole  photogenerated  current  collected  at  the  P+ 
layer.  Since  in  one  dimension  the  total  current  is  spatially  uniform 
and  since  the  P+  layer  cannot  inject  electrons  in  the  i-layer,  there 
must  be  a mechanism  that  compensates  for  the  difference  in  the  two 
photogenerated  currents.  This  mechanism  is  the  distortion  of  the 
electric  field  which  forward  biases  the  i/n  interface  and  reverse 
biases  the  p/i  interface.  The  forward  biasing  of  the  i/n  interface 
causes  an  electron  back- inj ection  current  that  compensates  for  the 
difference  between  the  two  photocurrents  and  reduces  the  total 
electron  current  at  the  N+  layer  to  the  same  level  as  the  hole  current 
at  the  P+  contact.  The  reverse  biasing  of  the  p/i  causes  a negligible 
hole  back  injection  current. 

The  conclusion  is  that  the  assumption  of  the  constant  electric 
field  can  be  applied  only  to  the  carrier  with  the  least  photocollected 
current.  If  we  apply  the  constant  field  approximation  to  the  other 
carrier,  we  ignore  its  back- inj ection  current  which  could  be 
comparable  to  the  photocollected  current.  Field  nonuniformities  do 
affect  the  current  of  the  carrier  with  the  least  photocollected 
current  but  their  effect  is  much  more  severe  on  the  carrier  with  the 
higher  photocollected  current. 
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We  stress  here  the  difference  between  a bipolar  device,  like  a 
p/i/n  solar  cell,  and  a device  like  an  n/i/n  photoconductor.  In  an 
a-Si:H  p/i/n  solar  cell  the  hole  current  at  the  P+  contact  must  be 
about  the  same  as  the  electron  current  at  the  N-*"  contact.  This 
forces  the  carrier  with  the  smaller  collection  length  to  determine  the 
performance  of  the  device  under  illumination.  On  the  contrary,  the 
carrier  collection  of  a n/i/n  photoconductor,  where  both  ohmic 
contacts  inject  the  same  type  of  carrier,  is  based  on  different 
principles.  To  facilitate  the  discussion,  we  can  imagine  an  n/i/n 
photoconductor  under  terminal  bias  and  uniform  optical  excitation. 
Again  we  assume  that  the  hole  collection  length  is  less  than  the 
electron  collection  length.  Here,  as  opposed  to  the  p/i/n  case,  for 
every  photogenerated  hole  trapping  event  in  excess  of  photogenerated 
electron  trapping,  an  electron  will  be  injected  from  the  negative 
contact  to  compensate  the  charge  of  the  trapped  hole . Hence , the 
current  is  proportional  to  the  total  distance,  per  unit  time,  covered 
by  all  photogenerated  carriers  before  they  recombine.  Thus,  in  a 
n/i/n  photoconductor,  the  carrier  with  the  longer  collection  length 
indeed  determines  the  current.  On  the  contrary,  in  the  p/i/n 
structure  the  P+  contact  cannot  inject  electrons  and  the  continuity  of 
the  total  current  is  preserved  through  forward  biasing  the  i/n 
interface  and  subsequent  electron  back  injection. 

The  conclusion  derived  here,  that  the  limiting  carrier  in  the 
p/i/n  structure  is  the  carrier  with  the  smaller  collection  length,  has 
also  been  reached  by  other  workers  [47] , [48] , [49] . They  based 
their  conclusion  on  the  minority-carrier  concept  borrowed  from 
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single -crystal  transport  theory.  However,  this  concept  is  not 
particularly  applicable  to  a-Si:H  p/i/n  cells  since  low  injection 
prevails  only  in  the  heavily  doped  layers. 

Hack  and  Shur  [42]  arrived  also  at  a similar  conclusion  regarding 
the  limiting  carrier.  They  show,  by  using  numerical  solutions,  that 
the  short-circuit  current  decreases  when  the  limiting  carrier  mobility 
decreases.  In  general,  their  numerical  solutions  are  among  the  most 
complete  in  the  literature.  They  have  taken  into  account  the 
continuous  spectrum  of  trapped  states  and  have  presented  numerous  data 
on  the  dependence  of  internal  device  parameters  on  absorption 
coefficient  and  band  mobilities.  Their  numerical  results  agree  with 
ours.  The  new  element  in  our  work  is  the  physical  interpretation  of 
the  numerical  solutions,  based  mainly  on  the  equivalent  circuit  model 
of  Fig.  5.5,  and  the  subsequent  analytical  formulation  presented  in 
Section  5.4. 

Numerical  solutions  have  also  been  presented  by  Swartz  [50] , 
Manfredotti  [51] , Schwartz  [52]  , where  the  amphoteric  nature  of  the 
dangling  bonds  was  taken  into  consideration,  and  Chen  [53].  Their 
solutions  added  valuable  insight  into  device  understanding.  At  the 
same  time,  each  of  these  solutions  is  based  on  different  assumptions 
and  device  structures.  Thus  no  unifying  view  of  these  numerical 
results  is  possible,  which  prevents  informed  general  conclusions  for 
the  device  operation.  It  is,  therefore,  important  that  consistent 
numerical  solutions  are  developed  that  will  enable  computer 
experiments  designed  to  reveal  the  approximations  and  insight  needed 
for  the  development  of  analytical  models  and  other  related  tools  such 
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as  equivalent-circuit  models.  All  this,  interwoven  with  experimental 
assessment,  will  lead  to  informed  design. 

5 . 6 Conclusions 

There  are  two  major  results  in  this  Chapter:  the  equivalent- circuit 
model  of  Fig.  5.5,  which  clarifies  the  concept  of  the  limiting 
carrier,  and  the  derivation  of  analytical  expressions,  5.19  and  5.25, 
for  the  cell  photocurrent  density.  A central  point  in  our  analysis  was 
the  separation  of  currents  into  photocollected  and  back- inj ection 
components.  The  carrier  with  the  least  photocollected  current  is  the 
limiting  carrier  and  determines  the  performance  of  the  device.  For 
strongly  absorbed  light  the  limiting  carrier  is  the  carrier  associated 
with  the  non- illuminated  contact.  Equation  (5.19)  expresses  the 
photocurrent  density  in  terms  of  the  interfacial  field  and  the 
limiting  carrier  collection  length.  For  uniformly  absorbed  light  the 
limiting  carrier  is  the  carrier  with  the  shortest  collection  length. 
The  Hecht  formula  can  be  used  in  the  case  of  uniformly  absorbed  light, 
if  a proper  assessment  of  the  field  dependence  on  terminal  voltage  is 
made.  The  equivalent  circuit  of  Fig.  5.5  can  have  practical 
ramifications.  When  designing  the  cell  care  should  be  taken  that  1^ 
and  I2  are  not  only  large  but,  also,  about  equal.  Therefore,  the 
tailoring  of  the  thin  film  composition,  either  in  terms  of  alloying  to 
change  the  band  gap  or  in  terms  of  doping  to  change  the  field  profile, 
should  focus  on  increasing  the  limiting  carrier  collection 
probability. 


CHAPTER  SIX 

CRYSTALLINE  VERSUS  AMORPHOUS  SILICON:  PHYSICS  AND  MODELING 

6 . 1 Introduction 

Although  the  fundamental  chemical  structure  (tetrahedrally  bonded 
silicon)  is  the  same  in  both  the  crystalline  and  amorphous  silicon, 
the  electronic  and  optical  properties  of  the  two  materials  are 
entirely  different.  The  major  differences  include  intrinsic  and 
extrinsic  dark  conductivity  and  photoconductivity,  optical  absorption 
coefficients,  photo luminescence  spectrum,  Hall  effect  and  thermopower. 
All  of  these  different  properties  originate  in  the  nature  of  the 
electronic  states  in  the  two  materials  and  can  be  traced  back  to  the 
methods  used  to  fabricate  them. 

In  single  crystal  silicon,  either  crystallized  from  molten 
silicon  or  obtained  through  epitaxial  methods,  most  of  the  native 
atoms  are  relaxed  in  their  minimum  energy  position.  The  high 
temperature  during  fabrication,  above  1200  °C,  gives  the  silicon  atoms 
the  ability  to  minimize  the  free  energy  of  the  lattice.  Well  defined 
point  and  line  defects  or  metal  impurities  give  rise  to  discrete 
energy  levels  in  the  band  gap.  On  the  other  hand,  amorphous  silicon  is 
deposited  in  thin  films,  either  by  decomposition  of  silane  gases  or  by 
sputtering,  at  relatively  low  temperatures,  200  - 500  °C.  The  silicon 
atoms  are  "frozen"  in  metastable  equilibrium  sites  and,  although  the 
short  range  order  of  tetrahedral  bonding  is  still  observed,  the  bonds 
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between  nearest  neighbors  are  stretched.  The  new  lattice  has  no  long 
range  order,  and  a continuous  spectrum  of  defect  states  is  introduced 
in  the  mobility  gap.  The  disruption  of  lattice  periodicity  relaxes  the 
crystal  momentum  conservation  rules  for  optical  transitions.  This 
increases  the  a-Si:H  absorption  coefficient  by  about  an  order  of 
magnitude  above  the  absorption  coefficient  of  crystalline  silicon, 
which  is  known  to  be  an  indirect-gap  material,  and  enables  the  use  of 
a-Si:H  for  thin  film  solar  cells.  The  lack  of  long  range  order  makes 
questionable  for  amorphous  silicon  both  the  applicability  of  Bloch's 
theorem  and  the  nearly  free  electron  approximation,  which  is  the 
basis  of  the  charge  transport  theory  in  crystalline  semiconductors. 

This  Chapter  does  not  provide  answers  to  the  theoretical 
questions  raised  above  concerning  amorphous  silicon.  We  will  assume, 
as  we  implicitly  did  in  Chapter  5,  that  the  free  electron 
approximation  is  still  a valid  hypothesis  and  that  the  charge  is 
transferred  through  drift  and  diffusion  of  extended  carriers.  Section 
6 . 2 focuses  on  the  charge  transport  mechanisms  in  the  two  materials 
and  explains  the  differences  in  terms  of  device  dimensions  and 
influence  of  localized  states.  Section  6.3  compares  the  modeling  of 
devices  and  discusses  the  applicability  to  amorphous  silicon  of 
concepts  such  as  low  injection  and  minority  carrier  borrowed  from 
silicon  device  modeling. 

6 . 2 Charge  Transport  in  Crystalline  and  Amorphous  Silicon  Devices 

An  excellent  example  that  illustrates  the  difference  in  the 
charge  transport  mechanisms  in  the  two  devices  is  the  celebrated  time 
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of  flight  experiment.  In  crystalline  silicon  this  experiment  is 
usually  referred  to  as  Haynes -Shockley  experiment.  Haynes  and  Shockley 
[54]  through  this  experiment  proved  the  existence  of  excess  carriers 
in  semiconductors  and  measured  the  minority-carrier  mobility  and 
lifetime.  In  this  experiment  a light  pulse  excites  a silicon  bar  (the 
original  experiment  was  done  on  germanium)  from  one  end  . A blocking 
contact  exists  near  the  other  end  while  two  other  contacts  and  an 
external  bias  create  an  electric  field,  E,  across  the  device.  The 
injected  pocket  of  minority  carriers  drifts  under  the  influence  of  the 
field  and  at  the  same  time  spreads,  because  of  diffusion,  and  forms  a 
Gaussian  function  in  space.  The  voltage  at  the  blocking  contact  will 
rise  when  the  package  of  the  minority  carriers  arrives  at  this 
contact.  Closed  form  expressions  for  the  time  delay  between  input  and 
output  are  used  to  extract  the  minority-carrier  diffusion  coefficient 
and  lifetime.  That  the  carrier  Debye  length  is  much  less  than  the 
length  of  the  bar  and  that  the  carrier  dielectric  relaxation  time  is 
much  less  than  the  transit  time  is  critical  to  this  behavior,  and  is 
related  to  the  quasi  neutrality  condition. 

In  amorphous  silicon  the  time  of  flight  experiment  is  done  on 
thin  film  p-i-n  devices  similar  to  the  one  analyzed  in  Chapter  5.  The 
device  is  reverse  biased,  to  increase  the  field  in  the  i- region  and 
make  it  more  uniform,  and  a pulse  of  strongly  absorbed  light  excites 
the  cell,  usually  through  the  P+  layer.  The  output  that  is  monitored 
is  the  transient  current  through  the  cell.  When  the  light  pulse 
strikes,  a pocket  of  electrons  and  holes  is  generated  at  the  p/i 
interface.  The  holes  are  immediately  swept  in  the  P+  layer,  while  the 
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electrons  drift  toward  the  N+  layer.  One  important  point  that 
distinguishes  this  experiment  from  the  previous  Haynes -Shockley 
experiment  is  that  there  is  no  delay  between  input  and  output:  the 
output  current  raises  to  its  maximum  right  after  the  excitation,  stays 
practically  flat  or  decays  slowly  until  the  electrons  arrive  at  the  N+ 
layer  and  then  drops  to  zero  with  a short  time  constant.  This  results 
from  the  fact  that  in  the  i- layer  the  transit  time  is  less  than  the 
dielectric  relaxation  time,  and  the  displacement  current  is 
communicated  at  the  speed  of  light  because  of  the  absence  of  carrier 
screening.  The  large  band  gap  and  the  numerous  free  carrier  traps 
result  in  low  carrier  densities  and  relatively  long  dielectric 
relaxation  times.  The  excess  photogenerated  electrons  cross  the  i- 
layer  before  the  holes  can  neutralize  their  charge.  During  the 
transient  the  current  in  most  of  the  i-layer  is  displacement  current 
caused  by  the  space  charge  of  the  moving  electrons. 

In  crystalline  silicon  the  minority  carriers  are  transferred 
through  a quasi-neutral  region  the  transit  time  of  which  is  much 
longer  than  the  dielectric  relaxation  time  due  to  the  numerous 
majority  carriers.  The  majority  carriers  shield  the  minority-carrier 
charge  and  the  two  pockets  move  together  until  they  arrive  at  the 
junction  of  the  blokcing  contact  where  the  strong  field  separates  the 
carries  and  causes  a current  to  flow  in  the  external  circuit.  In  a 
thin-film  amorphous  silicon  cell  the  whole  i-layer  is  a space-charge 
region,  a strong  electric  field  ( >104  V/cm)  extends  throughout  the 
device  and  the  separation  of  photogenerated  carriers  begins 
immediately  following  excitation. 
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Returning  to  the  time  of  flight  experiment,  we  identify  the  time 
interval  between  excitation  and  the  fast  output  current  decay  with  the 
transit  time,  t^,  of  the  device  [55].  If  E is  the  average  field  and  W 
the  device  thickness , then  one  could  write  the  relation 

^e=-lV  (6-1} 

for  the  electron  mobility.  However,  this  not  the  free-carrier  mobility 
but,  rather.it  is  an  effective  drift  mobility.  To  extract  the  free 
carrier  mobility  we  would  have  to  take  into  account  the  temporary 
trapping  of  the  carriers  especially  in  the  tail  states.  This  trapping 
delays  the  carriers,  makes  ne  in  (6.1)  less  than  the  actual  extended 
mobility  and  can  lead  to  dispersive  transport  [55] . The  derivation  of 
the  mobility  and  the  lifetime  from  the  time  of  flight  experiment 
depends  on  the  particular  model  for  the  interaction  between  free 
electrons  and  localized  states  or  phonons  [56,57,58],  Elaboration  on 
these  models  is  beyond  the  scope  of  this  Chapter.  We  would  only  like 
to  stress  that  in  crystalline  silicon  the  transport  properties  of  the 
free  carriers  can  be  extracted  directly  from  transient  experiments,  as 
was  done  in  Chapter  four.  This  can  be  contrasted  with  the  model - 
dependent  interpretation  of  transient  experiments  in  a-Si:H  devices. 

6 . 3 Modeling  Aspects  of  Crystalline  and  Amorphous  Silicon  Solar  Cells 
From  a mathematical  viewpoint  analytical  modeling  of  a bipolar 
device  is  equivalent  to  decoupling  the  original  six  Shockley  equations 
into  a system  of  linear  differential  equations  which  are  independent 
and  related  only  through  boundary  conditions.  This  linearization  of 
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the  system  allows  derivation  of  analytical  formulas  that  enable 
interpretation  of  experimental  results,  allow  understanding  of  the 
charge- transport  mechanisms  and  enable  derivation  of  equivalent - 
circuit  models. 

In  single  crystal  silicon  almost  all  the  modeling  is  based  on 
Shockley's  concept  of  considering  a device  as  consisting  of  quasi- 
neutral and  space-charge  regions.  Under  low  injection  the  field  and 
the  minority-  carrier  lifetime  are  bias  independent  in  the  quasi - 
neutral  regions  and  the  minority-carrier  density  can  be  described  by  a 
linear  second-order  differential  equation  (3.15).  A similar  equation 
can  be  derived  in  high  injection  by  using  the  ambipolar  transport 
equation.  This  makes  possible  the  calculation  of  recombination 
currents  and  carrier  densities  in  the  quasi-neutral  regions.  The 
recombination  in  the  space-charge  region  is  treated  by  making  the 
usual  assumptions  of  the  depletion  layer  approximation,  for  reverse 
bias,  or  the  quasi  equilibrium  approximation  for  forward  voltages.  For 
low  or  high  injection  in  the  quasi-neutral  regions,  the  linear  nature 
of  carrier  transport  in  silicon  solar  cells  allows  use  of  the 
superposition  principle  to  separate  the  current  into  independent  dark 
recombination  and  photogenerated  components. 

The  situation  is  quite  different  in  a-Si:H  p-i-n  solar  cells. 

Here  the  carrier  lifetimes  and  electric  field  are  bias  (terminal  or 
optical)  dependent.  We  can  still  derive  (3.15)  but  now  the 
coefficients  will  depend  on  the  carrier  densities,  which  renders  the 
system  nonlinear  despite  the  apparent  linearity  of  (3.15).  It  is 
clear,  therefore,  that  the  superposition  principle  is  not  valid  in 
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a-Si:H  solar  cells  and  this  is  pictorially  illustrated  in  the 
equivalent-circuit  model  for  amorphous  silicon  solar  cells  shown  in 
Fig.  5.5.  Although  a separation  of  currents  into  photocollected  and 
back- injection  components  is  still  possible,  as  was  pointed  out  in 
Chapter  5,  the  two  currents  are  not  independent:  the  back-  injection 
component  depends  on  the  generation  profile  and  also  the  cell 
photocurrent  depends  on  the  applied  bias. 

The  nonlinear  nature  of  the  system  of  equations  that  describe 
carrier  transport  in  a-Si:H  solar  cells  makes  necessary  introduction 
of  simplifying  approximations,  like  the  constant  field  and  lifetime 
approximation.  Under  these  approximations  we  can  derive  closed  form 
expressions  for  the  cell  photocurrent.  Such  expressions  were  presented 
in  Chapter  5,  where  it  was  emphasized  that  these  assumptions  should  be 
applied  only  to  the  limiting  carrier  to  avoid  neglecting  the  back- 
injection  current  of  the  non- limiting  carrier. 

In  a device  like  an  a-Si:H  p/i/n  solar  cell  the  only  regions  that 
can  be  thought  as  being  in  low  injection  are  the  heavily  doped  P+  and 
N+  layers.  From  (5.2)  and  (5.3)  follows  that  an  n(p)-type  region  is  in 
low  injection  if  n(p)  » Cp(n).  This  inequality  guarantees  that  the 
electric  field  and  the  minority  carrier  lifetime  are  bias  independent. 
The  numerical  solutions  of  Chapter  5 suggest  that  this  is  the  case  for 
N+  layer  under  reasonable  bias.  This  is  a result  of  the  high  (~ 
10^/cm^)  equilibrium  electron  density  in  the  N+  layer.  Therefore,  the 
recombination  losses  in  the  N+  layer  can  be  modeled  as 
In(exp (Fnp/kT) - 1) , where  Fnp  is  the  separation  of  the  two  Fermi  levels 
at  the  i/n  interface  and  In  the  saturation  recombination  current  of 
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the  layer.  In  the  P+  layer,  because  of  the  low  (~  5x10-*-^) 
equilibrium  hole  density,  the  strong  generation  rate  under  AMI . 5 light 
and  the  relatively  high  electron  lifetime,  the  space  charge  can  change 
with  bias.  Therefore,  despite  the  fact  that  in  this  layer  the  holes 
outnumber  the  electrons  and  the  electron  lifetime  is  nearly  bias 
independent,  the  electric  field  can  change  with  illumination  or 
electrical  bias.  Thus,  the  P+  layer  is  not  really  in  low  injection 
unless  the  light  is  reduced  and  the  electrical  bias  is  kept  low. 

6 . 4 Conclusions 

The  transient  and  the  steady-state  experimental  and  theoretical 
analysis  of  crystalline  and  amorphous  silicon  photovoltaic  devices 
reveals  different  charge- transport  mechanisms.  The  interaction  between 
free  carriers  and  traps  or  tail  states  determines  the  experimental 
behavior  of  amorphous  silicon  devices  especially  during  transients.  In 
crystalline  silicon  devices  the  carriers  are  transferred  mainly 
through  quasi-neutral  regions  and  the  separation  of  charges  takes 
place  at  the  p/n  junction  space-charge  region.  In  amorphous  silicon 
thin  film  devices  the  transit  time  is  shorter  than  or  is  comparable  to 
the  dielectric  relaxation  time  and  the  concept  of  quasi-neutrality  is 


not  applicable. 


CHAPTER  SEVEN 
SUMMARY 


Analytical  models,  computer  solutions  and  optimal  design  for 
solar  cells  were  the  major  contributions  of  this  research.  Analytical 
models  and  optimum  design  principles  were  developed  for  crystalline 
silicon  solar  cells.  Computer  solutions  were  obtained  for  amorphous 
silicon  solar  cells  to  understand  the  physics  of  the  device  and  to 
enable  analytical  modeling. 

In  silicon  solar  cells  the  effects  of  certain  three-dimensional 
geometries  were  assessed.  Texturizing  the  emitter  to  increase  the  base 
absorbance  results  in  the  same  base  recombination  as  in  the  one- 
dimensional device  with  planar  emitter.  Also,  partitioning  the  emitter 
in  a base-limited  solar  cell  leaves  the  open-circuit  voltage 
unchanged.  The  analysis  about  three-dimensional  effects  was  extended 
to  bipolar  diodes  and  transistors.  It  was  shown  that  for  certain 
planar  devices  the  transit  time  and  transport  factor  can  be  expressed 
in  terms  of  the  same  parameters  of  corresponding  one -dimensional 
devices.  These  results  have  implications  in  device  characterization 
since  they  allow  measuring  the  recombination  velocity  of  the  front 
surface  of  a planar  device. 

Optimum  design  principles  for  silicon  solar  cells  were  derived  by 
focusing  on  the  minimization  of  recombination  losses.  Lower  bounds  for 
the  recombination  current  were  obtained  in  terms  of  minority-carrier 
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bulk  transport  and  surface  parameters.  An  optimized  emitter  structure 
was  proposed  that  minimizes  recombination  under  the  metal  contact  and 
maximizes  current  collection  in  the  uncovered  areas . 

Concluding  the  modeling  of  silicon  solar  cells,  a new  method  was 
presented  and  demonstrated  that  allows  determination  of  the  minority- 
carrier  lifetime  and  diffusivity  in  the  base  of  the  device.  The  method 
exploits  the  Fourier  transforms  of  transient  responces  caused  by 
pulsed  optical  excitation  to  derive  the  frequency  responce  of  the 
base . 

Thin  film  solar  cells  made  of  amorphous  silicon  were  also 
modeled.  Contact-to-contact  computer  solutions  of  the  a-Si:H  p/i/n 
solar  cell  were  presented  and  treated  as  computer  experiments  designed 
to  reveal  the  approximations  and  insight  needed  for  the  development  of 
analytical  models.  The  numerical  results  allowed  study  of  many  aspects 
of  internal  variables  such  as  electric  field,  carrier  densities, 
current  densities  and  recombination  rates  versus  position,  terminal 
voltage,  and  photon  flux  density.  Strongly  absorbed,  weakly  absorbed 
and  AMI . 5 light  were  studied  to  gain  insight  needed  for  the  analytical 
modeling.  The  resulting  analytical  and  equivalent-circuit  models 
support  each  other  and  explain  the  physical  origin  of 

interdependencies  among  such  variables  as  quantum  efficiency,  electric 
field  and  recombination  rate  profiles  and  their  relation  to  current- 
voltage  characteristics . The  concept  of  the  limiting  carrier  was 
mathematically  treated  by  separating  the  current  into  photocollected 
and  back- inj ection  components.  This  leads  to  analytical  expressions 
for  the  cell  photocurrent  in  terms  of  the  limiting-carrier  collection 
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length.  The  limiting  carrier  is  the  carrier  with  the  least 
photocollected  current. 

Finally,  the  physics  and  the  modeling  of  crystalline  and 
amorphous  silicon  photovoltaic  devices  were  compared.  The  differences 
in  the  transport  properties  of  the  two  materials  are  interpreted  in 
view  of  the  Haynes -Shockley  and  time  of  flight  experiments.  In  silicon 
solar  cells  the  transit  time  is  much  longer  than  the  dielectric 
relaxation  time,  the  device  is  mostly  quasi  neutral  and  a 
photogenerated  carrier  has  to  arrive  at  the  junction  to  contribute  to 
the  terminal  current.  In  amorphous  silicon  solar  cells  the  transit 
time  is  comparable  to  the  dielectric  relaxation  time,  the  space  charge 
region  extends  throughout  the  device  and  a photogenerated  carrier  can 
contribute  to  the  external  current  by  causing  displacement  current. 


APPENDIX  A 

PROOF  OF  (3.18)  AND  (3.19) 
From  (3.8)  and  Rs  - 0 we  have 


dR 


op 


dx 


1 Rop^ 

(1 -)  > ( 

Neffr  Rb2 


Rop^ 

) a ;). 

^effr  m Rbm 


(A.  1) 


since  if  Rs  < R^m  then  R0p  < R^m  < R^.  The  solution  of  the  following 
differential  equation 

^ (1  - ) (A.  2) 

dx  (^effT)m  Rbm 

subject  to  the  boundary  condition  of  y(0)=0,  is 

y - Rbm  tanh(-f-)  , (A.  3) 

where  lm  - RbmAVNeff  Om-  From  (A.l)  and  (A.  2) 

R (x,Rs=0)  > y - Rbm  tanh(-y— ) . (A. 4) 

■Lm 

Both  functions  R]-)(N^)  tanh(x/L(N]_) ) and  R^m  tanh(x/L(N2) ) are  greater 
than  ROp(x,Rs=0)  since  they  are  the  functions  for  the  normalized 
recombination  current  of  a uniform  region  having  Rs  = 0,  N = N]_  = 

N[ (1/Nef£r )m]  and  N = N2  = N(R^m) , respectively.  By  definition, 
ROp(x,Rs=0)  is  less  than  either  one  of  these  two  functions.  This 


144 


145 


proves  (3.18).  In  a similar  way  one  could  prove  (3.19)  by  using  the 
relation 


<-^Rop  ’ Rs“°°) 
dx 


1 

Neffr 


Neff 

- -5—  Rop2(x>Rs=co)  > 

Neff 

( 5“)m  (Rbm2  " Rop  (x,Rs=“)) • 


(A. 5) 


APPENDIX  B 

DERIVATION  OF  PHOTOCURRENT  RELATIONS 


Here  we  apply  the  reciprocity  theorem  of  section  2.2  to  derive 
the  expression  for  the  photocurrent  of  an  one -dimensional  uniform 
region.  Let  L be  the  minority- carrier  (electrons)  diffusion  length  and 
S the  back  surface  (x  = W)  recombination  velocity  S.  The  expression 
for  the  minority-carrier  density  is 


n(x)  - A exp(-x/L)  + Bexp(x/L).  (B.l) 

The  boundary  condition  at  x - W is , 

Jn(W)  - e D dn/dx  = - e S n(W) . (B.2) 

Coefficients  A and  B are  determined  from  (B.2)  and  from  the 
normalization  condition  which  requires  n(0)  — 1.  Thus, 


1 exp(W/L)  + (LS/D) exp (W/L) 

2 cosh(W/L)  + (LS/D) sinh(W/L) ’ 

1 exp ( -W/L)  - (LS/D) exp ( -W/L) 

2 cosh(W/L)  + (LS/D) exp (W/L) ' 

Using  (2.9),  the  collection  probability  at  point  x is 

Pc(x)  = n(x)  - 

cosh((W-x)/L)  + (LS/D) sinh( (W-x)/L) 
cosh(W/L)  + (LS/D)  sinh(W/L) 
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Equation  (B.4)  can  be  extended  to  the  sinusoidal  steady  state  by 
replacing  L with  L'  (w)  = L/(l  + j wr)^,  where  r is  the  lifetime.  For 
a uniform  generation  rate  G,  the  photocurrent  in  the  frequency  domain 
becomes 


f W 

Jsc  (u)  = e G Pc(x,w) 

J0 


fW 

dx  = e G n(x,w) 

J0 


dx  - 


sinh(W/L' ) + (SL'/D)  cosh(W/L')  - SL'/D 
cosh(W/L' ) + (SL'/D)  sinh(W/L' ) 


(B . 5) 


The  last  relation  proves  (4.1).  Equation  (4.2)  can  also  be  derived  by 
setting  F' (w)  equal  to  Pc(W,w)  =*  n(W,w) • 

Focusing  now  on  a-Si:H  p/i/n  devices  let  us  assume  that  the 
carriers  not  only  diffuse  but  also  drift  under  the  influence  a 
spatially  uniform  electric  field  E.  Then,  from  the  continuity  equation 
and  the  expression  of  the  current  in  terms  of  drift  and  diffusion,  we 
have  for  the  limiting  carrier 

D ^-2-  + n E — — = 0.  (B . 6) 

dx  dx  r 


The  general  solution  of  this  equation  can  be  written  as 

n(x)  = A'  exp(x/L^)  + B'  exp(x/L2) , (B.7) 

where 

L]_  - [1/Le  + (1/LE2+  1/L.2 ) 1/2 ] " 1 , L2  = [1/LE  - (1/LE2  + l/L2)1/2]:1 

Here  LE=-2kT/qE  and  E>0.  For  infinite  recombination  velocity  at  x * M, 
n(W)  = 0 and  for  n(0)  = 1 (normalization  condition),  A'  and  B'  become 
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- exp (W/L9 -W/Lt ) 
-exp(W/L2-W/L1)  + 1 ’ 


-exp(W/L2-W/L1)  +1  ‘ 


(B.8) 


If  strongly  absorbed  light,  exp(-aW)  « 1,  shines  from  the  back 
surface,  then  the  photocurrent  becomes 


r w 

Jph  = e F a expa(x  -W)  Pc(x)  dx  = 

J 0 


f W 

eFa  I { [A'expCx/L!)  + B ' exp(x/L2) ]/exp( -qEx/kT) } exp [a(x-W) ] dx , (B.9) 

0 


where  exp (-qEx/kT)  is  proportional  to  the  equilibrium  electron  density 
of  a region  having  spatially  uniform  electric  field  E.  From  (B.9), 
(B.8)  and  exp(-aW)  « 1 (5.15)  can  be  derived.  Equation  (5.22)  can  be 
derived  from  (B.8)  and  the  following  expression  for  the  photocurrent 
for  uniformly  absorbed  light: 

fW 

Jph  = e G Pc(x)  dx  - 

J0 

f W 

e G { [A'exp(x/L1)  + B ' exp ( x/L2 )] /exp (-qEx/kT ) } dx.  (B.10) 

J0 


For  an  arbitrary  generation  profile,  G(x) , like  the  one  created  by 
AMI. 5 light,  the  photocurrent  is  given  by 

f W 

Jph  = e G(x) { [A'exp(x/L1)  + B ' exp (x/L2) ] /exp ( -qEx/kT) } dx.(B.ll) 
J0 


APPENDIX  C 

MINORITY- CARRIER  ACCUMULATION  AT  THE  BASE  EDGE  OF  A JUNCTION 
SPACE- CHARGE  REGION  UNDER  SHORT-CIRCUIT  CONDITIONS 

The  purpose  of  this  appendix  is  to  set  forth  accurate  boundary 
conditions  for  a short-circuited  p/n  junction  solar  cell  under 
illumination.  These  boundary  conditions  pertain  to  the  base  edge  of 
the  junction  space -charge  region  (SCR) . The  usual  assumption  is  that 
the  excess  electron  and  hole  densities  vanish  at  x = W.  This 
corresponds  to  an  assumed  infinite  effective  surface  recombination 
velocity  S.  Chapters  1 through  4 implicitly  used  this  assumption. 
Another  approximation  consists  in  setting  S equal  to  the  kinetic  limit 
or  unilateral  thermal  velocity  Sjq^  = Ut^/(67r) , where  Uth  is  the 
root  mean  square  thermal  velocity  of  the  electrons  considered  as  a 
Boltzmann  gas  [59], 

Here  we  derive  and  demonstrate  more  accurate  boundary 
conditions.  For  low  injection,  we  define  S as  the  ratio  of  the 
minority-carrier  diffusivity  to  the  extrinsic  Debye  length.  For 
moderate  or  high  injection  we  derive  a nonlinear  expression  for  the 
excess  hole  and  electron  densities.  The  same  expression  for  S holds 
also  as  the  boundary  condition  at  the  collector  edge  of  a 
quasi-neutral  base  (QNB)  of  a transistor  in  the  forward- active  mode. 
The  new  value  of  S is  a factor  of  six  smaller  than  Sj^l  for  lightly 
doped  ( = 10^  cnT^)  silicon,  which  has  practical  consequences 
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particularly  for  transient  experiments  used  in  determining  carrier 
lifetimes  and  diffusion  coefficients.  For  higher  doping  densities  the 
derived  expression  for  S converges  to  SKL  and  the  effect  of  the 
accumulation  of  minority  carriers  at  the  SCR  edge  of  the  base  is  found 
negligible.  The  results,  therefore,  reached  in  the  first  four  Chapters 
are  still  valid  since  they  were  derived  for  regions  in  low  injection, 
which  implies  at  least  a moderate  doping  density  ( > 10^  Cm'^) . Thus, 
the  expressions  derived  here  become  important  only  for  lightly  doped 
regions  in  high  injection. 

Figure  C.l  shows  an  illuminated  n/p  junction  solar  cell  under 
short-circuit  conditions.  The  equation  for  the  electron  current  is 

Jn(x)  ” e Dn(x)dn(x)/dx  + e /xn(x) E(x)n(x)  , (C.l) 

a first-order  linear  differential  equation  for  n(x) . By  integrating 
this  equation  [60]  we  can  express  the  electron  density  in  the  SCR  in 
terms  of  the  electric  field  and  electron  current  as  follows: 

n(x)  - exp[  -J*Q(e/kT)E(x'  )dx'  ] (n(0)  + 

Jo(Jn(x')/eDn(x'))exp[jQ(e/kT)E(x")dx"]dx' ) , (C.2(a)) 

where  the  Einstein  relation  was  used. 

As  shown  in  Fig.  C.l,  the  point  x = 0 is  the  emitter  edge  of  the 
SCR  and  n(0)  is  the  majority  carrier  density  at  x = 0.  Taking  V(x)  to 
be  the  electrostatic  potential  in  the  electron  energy  scale  referenced 
to  V at  x = 0,  we  can  write  (C.2(a))  as 
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Fig.  C.l  An  illuminated  solar  cell  under  short-circuit  conditions. 

The  electron  accumulation  at  x = W is  emphasized  through 
the  separation  of  the  quasi-Fermi  levels,  Epn  - Epp. 
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n(x)  = exp( -eV(x)/kT) [n(0)  + 

J*(Jn(x')/eDn(x'))exp(eV(x')/kT)dx'].  (C.2(b)) 

If  W is  the  SCR  thickness  and  Vb  is  the  SCR  barrier  height,  then  from 
(C.2(b))  the  electron  density  at  the  edge  of  the  base  becomes 

n(W)  = exp( -eVb/kT)n(0)  + 

exp(-eVb/kT)  /Q(Jn(x)/eDn(x) )exp(eV(x)/kT)dx.  (C.3) 

We  now  examine  separately  the  low  injection  and  the  moderate  or  high 
injection  cases  in  the  base.  The  emitter  is  heavily  doped  and  is  in 
low  injection  for  any  level  of  excitation.  First  we  deal  with  low 
injection. 

If  low  injection  prevails  in  the  base  and  if  we  neglect  series 
resistance,  then  in  (C.3)  Vb  becomes  Vb^,  where  Vb^  is  the  SCR 
built-in  voltage  in  equilibrium.  Since  nQ(W)  = exp( -eVb^/kT)n(0)  is 
the  equilibrium  minority- carrier  density  at  the  base  edge  of  the  SCR, 
the  excess  minority-carrier  density  at  this  edge  becomes 

Sn(W)  - exp(-eVbi/kT)  (Jn(x)/eDn(x) )exp(eV(x)/kT)dx.  (C.4) 

Since  the  term  exp(eV(x)/kT)  increases  dramatically  with  x for 
0 < x < W while  Jn(X)  and  Dn(x)  are  much  weaker  functions  of  x,  the 
integrand  in  (C.4)  has  the  qualitative  shape  pictured  in  Fig.  C.2. 
Here,  Sx  is  a small  distance  in  the  vicinity  of  x = W where 
exp(eV(x)/kT)  stays  relatively  constant.  Therefore,  (C.4)  becomes 
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q D(W) 


e 


Fig.  C.2  The  quantity  (Jn/eD)exp(eV/kT)  as  a function  of  the 
distance  in  the  space  charge  region. 


Fig.  C.3  Band  diagram  near  the  base  edge  of  the  space  charge  region. 
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5 n(W)  « exp(-eVbi/kT)(Jn(W)/eDn(W))exp(eVbi/kT)5x, 


Sn(W)  « (Jn(W)Sx/eDn(W)) 


(C . 5 (a) ) 


and  the  effective  recombination  velocity  can  be  written  as 


S = Dn(W)/Sx. 


(C . 5 (b) ) 


We  now  prove  that  Sx  equals  the  extrinsic  Debye  length  in  base 
at  x = W.  We  define  the  SCR  edge  of  the  base  as  the  point  where  the 
space  charge  density  is  one-tenth  of  the  net  doping  density.  Let 
V' (x)  be  the  difference  between  the  electrostatic  potential  at  x = W 
and  that  at  any  point  in  the  bulk  of  the  base  (Fig.  C.3).  Since  at 
x =>  W the  density  of  the  majority  holes  is  nine -tenths  of  the  bulk 
value  (the  net  doping  density  is  assumed  to  vary  negligibly  within  an 
extrinsic  Debye  length  Lb  in  the  quasi-neutral  base)  then  V' (W)  will 
obey 


where  nondegenerate  statistics  have  been  used.  Since  V' (W)  is  very 
small  it  changes  exponentially  with  distance  for  a few  Debye  lengths. 
Hence  near  x - W 


0.9  - exp(-eV' (W)/kT)  - V' (W)  « 0.1  kT/e , 


(C.6) 


V'(x)/V'(W)  * exp [ - (x-W)/Ld] 


(C.7) 


and  the  potential  V(x)  as  used  in  (C.2(b))  becomes 


V (x)  = Vbi  - V'(x)  = Vbi  - V'(W)exp[-(x-W)/LD] , (C.8) 


in  the  vicinity  of  W.  From  (C.6)  and  (C.8),  we  write  exp(eV(x)/kT)  as 
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exp(eVbi/kT)[l  - exp( - (x-W)/LD) (eV' (W)/kT) ] » 


exp(eVbi/kT)[l  - 0 . lexp( - (x-W)/LD) ] 


(C.9) 


near  x - W.  Therefore,  exp (eV/(x)/kT)  remains  substantially  constant 
over  a Debye  length  at  x = W.  Consequently  Sx  ~ Lq  and  from  (C.5(b) 
the  boundary  condition  becomes 


If  the  base  was  n-type  having  a Debye  length  Lb  then  the  recombination 
velocity  for  holes  would  become  S - Dp/Lp.  Equation  (C.10)  holds  also 
under  reverse  bias  if  5n(W) , as  given  by  (C.5),  is  much  higher  than 
the  equilibrium  minority  carrier  density  at  x - W. 

We  now  switch  our  attention  to  high  injection.  Under  moderate  or 
high  injection  the  electric  field  in  the  base  changes  and  causes  the 
SCR  barrier  height  to  alter  by  an  amount  SVb£  = Vb^  - Vb . Using  the 
analysis  presented  in  the  low  injection  case,  we  can  write  (C.3)  as 

5 n(W)  - n(W)  - nQ(W)  = nQ(W) [exp(eSVbi/kT)  - 1] 


where  nQ(W)  is  the  equilibrium  electron  density  at  x = W.  In  (C. 11(a)) 
Lb  is  given,  for  Boltzmann  statistics,  by 


S * Dh/Lh. 


(C.10) 


+ Jn(W)LD/(eDn), 


(C. 11(a)) 


Ld  = J [kTe/(p(W)  + n(W) ) ] / e, 


(C. 11(b)) 


where  p(W)  and  n(W)  are  the  hole  and  electron  densities  at  the  SCR 
edge  of  the  base  and  e is  the  dielectric  constant.  Equation  (C. 11(b)) 


156 


serves  as  an  accurate  nonlinear  boundary  condition  for  an  illuminated 
short-circuited  p/n  junction  solar  cell  for  which  the  base  is  in 
moderate  or  high-level  injection.  It  relates  the  hole  and  electron 
densities  and  the  electron  current  density  at  the  space-charge-region 
edge,  x = W,  through  a variable,  SV^,  that  accounts  for  all  electric 
potential  variations  in  the  quasi-neutral  base  and  for  series 
resistance  effects  in  the  heavily  doped  emitter. 

Returning  now  to  low  injection,  equation  (C.10)  resembles  the 
expression  S = D/W  for  a uniformly  doped  short-base  region  of 
thickness  W,  defined  by  the  absence  of  bulk  recombination  and 
terminated  by  an  ideal  ohmic  contact.  For  both  (C.10)  and  for  S - 
D/W,  the  excess  minority-carrier  density  is  forced  to  vanish  within  a 
short  distance  characteristic  of  the  problem  under  study.  In  the  case 
of  the  short  base  region,  this  distance  W separates  the  injecting 
boundary  from  the  ideal  ohmic  contact.  In  the  case  of  (C.10),  this 
distance  Lp  separates  the  space -charge -region  edge  of  the 
quasi-neutral  base  from  the  high  field  part  of  the  space-charge 
region.  Equation  (C.10)  holds  only  if  Dn/Lp  < S^;  otherwise  S = S^. 

Here  we  emphasize  that,  provided  Dn/Lj}  < S^,  the  effective 
recombination  velocity  is  not  the  kinetic-limit  velocity.  Rather  it 
is  determined  by  the  rate  at  which  the  minority  carriers  diffuse 
toward  the  high- field  region  of  the  junction,  a rate  measured  by  the 
velocity  Dn/Lj). 

Equations  (C.10)  and  (C.ll)  serve  as  a boundary  condition  for  an 
illuminated  solar  cell.  They  also  apply  to  the  collector  edge  of 
quasi-neutral  base  of  a transistor  in  the  forward-active  mode, 
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provided  the  collector  is  in  the  low- current  regime  defined  by  J/e  « 
Nc  Vgat  which  is  the  product  of  the  collector  doping  and  the  saturated 
velocity. 

If  p is  the  majority-carrier  density  at  the  edge  of  the 
quasi-neutral  base  then,  for  Boltzmann  statistics  and  low  injection 
(C.10)  becomes 

S « eDny  [p/kTe],  (C.12) 

for  electrons.  A similar  equation  holds  for  holes  in  an  n-type  base. 
For  moderately  doped  silicon,  the  value  of  S in  (C.12)  is  about  the 

C. 

same  as  the  kinetic  limit  velocity  Sj^l  for  electrons,  Sj^l  = 5x10° 
cm/s,  or  holes,  SKL  = 3x10^  cm/s,  [61].  For  doping  densities  higher 
than  1017  cm'^  the  kinetic  limit  velocity  imposes  an  upper  bound  on  S. 
For  light  doping  densities,  NA  = 10^  cm"^,  (C.12)  yields  a value  six 
times  less  than  S^.  Also,  at  low  temperatures  the  ratio  S/S^  « 1 
since  Sj^l  an<^  D decrease  according  to  power  laws  while  deionization 
causes  the  majority  hole  density  to  decrease  exponentially  with 
temperature . 

We  now  investigate  the  effect  that  a finite  value  of  S has  on  the 
short-circuit  current  of  a solar  cell.  The  excess  minority-carrier 
density  at  the  space-charge  region  edge  of  the  quasi-neutral  base  is 
5n  =Igc/eS.  This  particle  accumulation  forces  a component  of  current 
to  flow  in  a direction  opposite  to  that  of  the  photocurrent, 

51  = IoBSn/n0  - IoBIsc/(eSn0),  (C.13) 

where  IoB  is  the  saturation  value  of  the  base  recombination  current. 
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Let  Igc*  t>e  the  short-circuit  current  of  the  same  cell  when  S = ®. 

Low  injection  implies  linearity  and  the  superposition  principle 
allows  Igc*  to  be  expressed  as 

ISC>  = ^SC  + 51  = Igc  [1  + IoB/^eSno)l-  (C.14) 

Therefore,  the  fractional  loss  in  short-circuit  current  arising 
from  the  finite  value  of  S is 

IoB/(eSn0)  - (IoB/eni2)N(W)/S,  (C.15) 

where  N(W)  is  the  doping  density  at  the  space-charge  region  edge  of 
the  quasi -neutral  base. 

O 

In  present  silicon  solar  cells  the  ratio  IoB/en£4  is  at  least 
10*^  cnA/s  • Thus,  even  if  N(W)  becomes  5x10^  cm'^  equation  (C.15) 
indicates  that  the  relative  current  loss  of  no  more  than  1%.  Similar 
results  hold  for  GaAs  solar  cells.  Consequently,  the  implicit 
assumption  of  infinite  S that  was  used  in  Chapters  1-4  is  a realistic 
one . 

Using  the  formula  <5n/n0  - IgQN(W)/(eSn^2)  and  realistic  values 
for  IgQ  and  N(W) , then  for  silicon  under  one-sun  illumination  and 
short-circuit  conditions  the  separation  of  the  quasi-Fermi  potentials 
at  the  space-charge-region  edge  of  the  base  is  about  0.44  V ~ 17  kT/e 
at  T - 300  K.  For  a GaAs  cell  this  separation  is  about  0.7  V = 

27  kT/e.  It  is  likely,  therefore,  that  the  recombination  in  the 
space -charge  region  will  be  more  important  than  the  back  injection 
into  the  base  especially  at  low  temperatures.  Thus,  the  effect 
discussed  could  appreciably  decrease  the  short-circuit  current  of 
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polycrystalline  silicon,  GaAs  or  other  solar  cells  having  a large 
density  of  recombination  centers  in  the  junction  region. 
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